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GAS TURBINE FACE SEAL THERMAL DEFORMATION AND COMPUTER PROGRAM 
FOR CALCULATION OF AXISYMMETRIC TEMPERATURE FIELD 


by Terrence E. Russell, Gordon P. Allen, Lawrence P. Ludwig, 

and Robert L. Johnson 

Lewis Research Center 
SUMMARY 

The effects of seal deformation on seal balance are pointed out. In particular, di- 
vergent sealing faces are shown to be undesirable. As a first step in determining seal- 
ing face deformation, a computer program is presented for calculating the temperature 
field in the arbitrary body of revolution with arbitrary thermal boundary conditions. The 
program, which is written in FORTRAN IV language for an IBM 7094 digital computer, 
is general and can be applied to a variety of problems. Application of the program to a 
specific gas turbine mainshaft seal is demonstrated in detail; both computer input and 
output data are given. The program is also used to solve an example problem which has 
an exact solution; thus a partial check on program accuracy is provided. 


INTRODUCTION 

Face seals in advanced gas turbine engines will be subjected to temperatures, pres- 
sures, and sliding surface speeds higher than that of current engine practice. Increases 
in temperature arise from continuing increases in flight speed and turbine inlet temper- 
ature (ref. 1), increases in pressure result from use of higher pressure ratio compres- 
sors, and increases in surface speeds are due in part to increases in shaft diameter 
needed to carry increases in torque. As an example, proposed advanced engines have 
seal sliding speeds to 500 feet per second (127 m/sec) (ref. 2) as compared to 350 feet 
per second (88. 9 m/sec) in engines in current use. In order to cope with these higher 
temperatures, pressures, and speeds, the contacting-type face seal may be replaced 
with noncontacting face seals such as the hydrostatic and hydrodynamic type (ref. 3). 
However, regardless of the type of face seal used, deformation of the sealing faces, due 
to thermal gradients, pressures, and centrifugal forces, will have significant effects on 
seal performance (ref. 4). In particular, deformations due to thermal gradients are a 



major problem area. The first step in determining thermal deformation is calculation 
of temperature distribution in the seal assembly. Such calculation by computer program 
is the subject of this report. The specific objectives of this study were to (1) present in 
detail a computer program for calculation of temperature distribution in axisymmetric 
bodies, ( 2) show how the computer program is applied to mainshaft seal thermal analy- 
sis, and (3) demonstrate program accuracy by comparison with a thermal problem having 
an exact solution. 

The computer program for calculation of temperature distribution is written in 
FORTRAN IV computer language for use on an IBM 7094 digital computer. The program 
is quite general and can be applied to a variety of axisymmetric body problems. The cal- 
culation procedure requires that these bodies be divided into an arbitrary finite number 
of axisymmetric volume elements or nodes. These nodes need not be equal in cross sec- 
tion. The program will take into account contact resistance at the interface between 
nodes and will also account for material properties that vary from node to node. Also 
provisions are made in the program to handle varying gas temperatures along the seal 
boundaries and internal viscous heat generation within the fluid at the boundaries. 


SYMBOLS 

A area, ft^ 

2 

A x cross section area, ft 

a outer radius of rotor or stator, ft 

C thermal conductance, Btu/(sec)(°F) 

Cj skin friction coefficient 

C M torque coefficient 

C specific heat, Btu/(lbm)(°F) 

Jr 

D hydraulic diameter, ft 

D l characteristic length for liquid film cooling, ft 

d differential 

F function of 

f view factor for radiation 

Gr Grashof number 

Gr D Grashof number based on diameter, D^g o (T^ - T g ^/i/j ’ T f 
Gr L Grashof number based on length, L 3 g Q (T b - T s )/^f T f 
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So 

H 

h 

K 

k 

L 

AZ 

N 

Nu 

P 

w 

Pr 

Pr f 

Q 

q 

q f 

R 

R c 

Re 

Re c 

Re y 

Re 

w 

Re x 

Re u 

r 

Ar 


gravitational acceleration, ft/sec 
convective conductance, Btu/(sec)(°F) 
heat transfer coefficient, Btu/(ft^)(sec)(°F) 

/3/co 

thermal conductivity, Btu/(ft)(sec)(°F) 
length of cylinder or surface, ft 
axial length, in. or ft 
radiative conductance, Btu/(sec)(°F) 

Nusselt number, hD^/k 
wetted perimeter, ft 
Prandtl number 

Prandtl number of boundary film, [v /k l 
heat flux, Btu/sec 
internal heat generation, Btu/sec 
fluid internal heat generation, Btu/sec 
thermal resistance, (sec)(°F)/Btu 
contact resistance, (sec)(°F)/Btu 
Reynolds number 

critical Reynolds number for transition 
Reynolds number based on total velocity, 
Reynolds number based on mass flow, W D/ /j. f 

a I 

Reynolds number based on flow length, V^X/ 
Reynolds number based on rotation, cor / 
radial coordinate, ft 
radial width, in. or ft 
average radius of element, in. or ft 
inner radius of annulus between cylinders, ft 
outer radius of annulus between cylinders, ft 
clearance (e.g. , seal gap), ft 
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temperature, °F 

average (or film) temperature, °F 
temperature entering fluid element, °F 
temperature leaving fluid element, °F 


bulk fluid temperature, °F 
Taylor number of boundary film, 


cors^^/ i 


axial velocity, ft/sec 


linear velocity of inner cylinder, ft/sec 


total fluid velocity, ft/sec 


free stream velocity, ft/sec 
mass flow rate, lbm/sec 

o 

mass flow per square foot of cross section, W/A x , lbm/(ft )(sec) 
over- relaxation factor 


X flow length, ft 

z axial coordinate, ft 

a coefficient of thermal expansion 

|3 angular velocity of fluid, rad/sec 

3 partial derivative 

A difference 


€ emissivity for radiation 

p absolute viscosity, lbm/(ft)(sec’ 1 

v kinematic viscosity, ft /sec 

O 

p density, lbm/ft 

2 

Tx tangential shear stress at wall, lbf/in. 
u) rotor angular velocity, rad/sec 
Subscripts: 

b at fluid bulk temperature 

c convective 


e environment 

ei between element and environment 
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es from environment to element 

f at boundary film temperature 

i element 

j neighbor element of i 

ji between neighbor and element 

k environment ’’area” affecting element i by radiation 

ki between environment "area" k and element 
Z environment ’’area" affecting element i by convection 

Zam laminar 

Zi between environment ’’area" Z and element 
r radiation 

s surface 

se from element to environment 

tr transition 

turb turbulent 


DISCUSSION OF THERMAL DEFORMATION IN SEALS 

Figure 1 is a schematic drawing of a gas turbine mainshaft seal for bearing sumps. 
The seal restricts leakage of surrounding high temperature, high pressure gas into the 



Figure 1. - Mainshaft seal for bearing sump in gas turbine engine. 



sump. The gas leakage, through the seal into the sump, acts as a sweep gas which helps 
prevent oil leakage. However, it is desirable to minimize this gas leakage, since the 
hot gas degrades the lubricant and adds to the cooling requirements of the lubricant sys- 
tem. 

Seal leakage and wear are grossly affected by seal force balance and it has been 
recognized that deformation of the sealing faces (which form the dam) can have a signifi- 
cant effect on the seal force balance (ref. 5). A discussion of seal force balancing can be 
found in reference 6. The manner in which the sealing face deformation affects the seal 
force balance and stability of a conventional face seal is discussed in reference 7. In 
general, thermal deformations cause the seal faces (dam) to form a nonparallel leakage 
path, and the pressure within these sealing faces is greatly dependent on the shape (con- 
vergent or divergent) of this leakage path; thus the force balance is affected by the ther- 
mal deformation. 

Thermal gradients usually produce sealing gaps divergent from the dam inside diam- 
eter to the outside diameter. Thus, it is usually advantageous to locate the higher pres- 
sure at the seal outside diameter since this will result in convergence with respect to the 
leakage direction. However, many face seals are designed with the higher pressure at 


Parallel faces Divergent faces Seal rub caused 


by divergence 



300 400 500 600 700 800 

Sealed air temperature, K 

Figure 2. - Effect of thermal deformation on air leakage of orifice 
compensated hydrostatic seal. Sliding speed, 200 feet per second 
(61.0 m/sec); Dressure differential, 100 pounds per square inch 
(6. 89X10 4 N/rA (Data from ref. 4. ) 
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the dam inside diameter and the gas turbine face seal is an example of this construction 
practice. The advantage of placing the higher pressure gas at the dam inside diameter 
puts the oil mist at the outside diameter; thus centrifugal force also acts to prevent oil 
leakage. 

Thermal deformation also affects operation of hydrostatic seals. Figure 2, from ref- 
erence 4, shows the effect of sealing face deformation on leakage through a hydrostatic 
seal operating at 200 feet per second (61. 0 m/sec) and 100 pounds per square inch 

(6. 89x10^ N/m^) pressure differential. With a sealed air temperature of 120° F (322 K), 

0 

the seal leakage is approximately 13 SCFM (0.0062 standard m /sec) and as the air tem- 



Figure 3*. - Overall wear pattern, surface profile traces, and photomicrographs of tungsten carbide seal seat. Sodium temperature, 1000° F (538° C); 
operating time, 4 hours; sliding velocity, 79 feet per second (24 m/sec); pressure, 50 pounds per square inch gage (34 N/cmZ gage). }ata from ref. 6.) 
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perature is increased, the leakage decreases. This leakage decrease is due to increas- 
ing angular deformation of the sealing gap. The result is that the seal runs with closer 
clearance. Eventually as the air temperature increases still further, the deformation 
causes the inside diameter of the nosepiece to rub against the seal seat and failure occurs. 

Seal thermal deformation effects have also been noticed in other applications. As an 
example, reference 6 reports that the main problem in sealing liquid sodium with a face 
contact seal was thermal deformation of the sealing faces. Some results of this study are 
repeated in figure 3 which shows the wear pattern made by a nosepiece sliding against a 
seal seat. Surface profile traces have been taken at intervals around the wear track. 
These profile traces show that the sealing faces were not parallel during operation, the 
greatest wear occurring at the inside diameter. The slope of the wear pattern suggests a 
divergence in the range of 0. 006 inch per inch (0. 006 m/m) of radial distance. 

Regardless of the seal type (face contact, hydrostatic, or hydrodynamic), the sealing 
gap is usually small (0.0001 to 0.0010 in. or 0.00025 to 0.0025 cm for gas film seals) 
and deformations can easily be of the same magnitude, thus affecting seal performance 
significantly. These sealing face deformations could arise from temperature, pressure, 
or centrifugal force effects. Careful design can minimize undesirable effects due to 
pressure and centrifugal force. Deformations due to temperature gradients are usually 
difficult to eliminate. Even the heat generated in shearing a thin film of gas within the 
interface can introduce a significant and undesirable thermal gradient in the seal rings. 

The preceding evidence indicates that thermal deformation can significantly affect 
seal performance. The first step in analysis of these thermal deformations is calculation 
of the temperature fields in the seal assembly, and the computer program for these calcu- 
lations is given in the following sections. 


ANALYSIS OF TEMPERATURE FIELDS IN AXISYMMETRIC BODIES 

Shaft seals are composed of basically axisymmetric bodies. Also, the circumferen- 
tial temperature gradient approaches zero for most applications. Therefore, the cylin- 
drical coordinate system is used as a basis for analysis. The following restrictions are 
placed on the thermal analysis: 

(1) A steady state must exist. 

(2) An axisymmetric temperature field must exist. 

Therefore, with internal heat generation, the heat conduction equation in cylindrical 
coordinates is (ref. 8) 


1 _d_ 

r Sr 




+ q(r, z) = 0 


( 1 ) 
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r 



Figure 4. - Seal axisymmetric nosepiece subdivided into 
finite number of three-dimensional axisymmetric 
volume elements. 



where r is the radial coordinate, k is thermal conductivity, z is the axial coordinate, 

T is temperature, and q is internal heat generation. 

Equation (1), which makes no restrictions of constant thermal conductivity k, can be 
approximated by a finite difference numerical approach in which thermal conductivity is 
held constant over small regions. In this finite difference approach, each of the axisym- 
metric bodies, which comprise the seal, is subdivided into a finite number of three- 
dimensional axisymmetric volume elements of rectangular cross section (see fig. 4). 
These elements need not be equal in cross-sectional area. In general, there will be P 
elements of which m are elements whose temperature is obtained by solving a heat bal- 
ance equation and P - m are elements in which the temperature is specified. Figure 4 

J.U 

indicates the position of a typical interior element (i Ln element). The numbers assigned 
to the neighboring elements will be determined in the construction of the element pattern, 
but for the purpose of developing the heat balance equation these neighboring elements are 
numbered 1, 2, 3, and 4 (see fig. 5). 

For the numerical solution, the application of the law of conservation of energy at 
each volume element allows good generality to be maintained. Referring to figure 5, a 
heat balance for steady- state conduction at element i results in the following difference 
equation: 
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Q li + Q 2i + Q 3i + Q 4i + q i “ 0 


where the 


terms are the heat fluxes from the neighboring elements or 


C li< T l - T i> + C 2i< T 2 ' T i> + C 3i< T 3 - T i> + C 4i< T 4 - T i> * «1 = 0 


(2a) 


in which the conductances are computed as follows with reference to figure 5: 


Az. 


Az. 


R li = 


'li 


27tRj Ar-jkj 


+ R. 


2ttR. Ar.k. ^li 


AR. 


AR, 


R 2i = 


'2i 


2 jrR mi Az i k i 


+ R„ 


2nR m2 Az 2 k 2 2i 




J 


(2b) 


where R.. is thermal resistance. As an example, for heat conduction from neighbor 


n 


element 1 to the i 
.th 


.th 


that of the i 


element, the resistance is made up of the resistance of one-half 
element plus one-half that of neighbor element 1 plus the contact resis- 


tance R . Thus the calculated temperature is that of the element center. However, 

ji 

the thermal conductivity assigned to the whole element volume is the value at the temper- 
ature of the element center. 

Equation (2a) can be stated concisely as 


2m - v + = ° 


(3) 


j=l 


where i is the element under consideration, the j's are its four neighboring elements, 
and Uism. At the outer surfaces of the axisymmetric bodies, which comprise the 

it. 

seal, the i a element will have less than four solid neighbor elements but will have en- 
vironmental radiation and convection effects on the external surface (or surfaces) (fig. 6). 
This means that some of the j neighbors are replaced by environmental regions which 
are identified as the k Ln area for radiation and the l U1 area for convection. It should 
be noted that an external surface could have both radiation and convection effects. The 

f Vi 

k m and l 1 areas would, therefore, be superimposed. 
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J-2 

neighbor 
node # 

Q q : 

j =3 

neighbor 
node i 

i node 


x ei 


j =4 

neighbor 
node # 


Environment 


z 


Figure 6. - Typical element with environment convection and/or radiation 
to one surface. 


th 

An i element could have radiation on more than one face. This radiation heat 
flux is 

4 

* Y N “ (Tk " Tl) 

k=l 

one face with convective effects 
4 

Y H « (T * ‘ T i ) 

1 = 1 

Combining equations (3), (4), and (5) results in the difference equation for the heat bal- 
ance of the i th element: 


An i^ 1 element could have more than 


Q Zi“ 


( 6 ) 


4 


4 


4 


X C i‘ <T i ’ T i ) + Z Nkl<Tk ' T i ) + Z Ha<T * ’ T ‘ 

dzl 




Conduction 


V 

Radiation 




2=1 


Convection 


+ -qj = 0 


7 

Internal 

heat 

genera- 

tion 


It should be noted that equation (6) establishes the convention that heat flowing from an 
element has a negative value and heat flowing into an element has a positive value. 

In addition to the boundary conditions of radiation and convection, two other boundary 
surface conditions arepossible. The heat flux may be a specified function 


r 


Q ■< 

^ei 



( 7 ) 


or 

C(r,z) 


where C(r, z) = 0 for a perfectly insulated surface. Or the heat transfer coefficient may 
be a specified value such as 


h ei = f(r, z) .(or H ei = f(r, z)) (8) 

Also the temperature of any element may be specified as 

T. = f(r, z) (9) 

The program has the special feature of calculating varying gas temperatures in flow 
over a surface or in small passages. In addition, heat generated within a fluid, such as 
that due to viscous shear in the sealing interfaces, can be accounted for. Referring to 
figure 7 reveals the heat balance for the center fluid element is 

WC p (T in - Tout) + h l A l< T sl - T i„> + h 2 A 2< T s2 ‘ T in> + % = 0 < 10 > 
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Solid body elements! 


Annular neighbor | 
(of flow element) . 


Fluid 

flow 


m f c P T fn 


I 

j Flow neighbor 
I (of flow element) 

I 

b -- 


h 2 A 2 (Ts2 - T in ) 



“in 




T out 

I m f c p T out 


■ n ‘^ Ts ‘ ; 

h l A l(T S l - T fn > 

Flow element t 
Solid body elements -|— 


Figure 7. - Fluid volume element for varying temperature forced 
convection heat transfer. 


in which q^ is the internal heat generated. Equation (10) enters into the element heat 
balance (eq. (6)) only through the determination of the fluid temperature. 

Returning to the element heat balance (eq. (6)) matrix notation can be used to express 
this system of linear equations as follows: 



[A] {T} {TL} 

(ID 


where [A] is the m x m conductance matrix, {T} is the m x 1 temperature vector, and 
{TL } is the m x 1 thermal load vector. 

In equations (6) and (11) it should be noted that, for interior elements, the N, . and 
terms vanish; also numbers are assigned to the j elements. 

The matrix [A] can be shown to be real, symmetric, and positive definite. The tem- 
peratures are obtained using the successive over- relaxation technique. By introducing 
the over- relaxation factor w, the successive over- relaxation algorithm is given for the 
temperature at the i th element at the (K + l) th iteration by: 
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where 


a. ..Tp 1 = a, .T? + W 


i-1 


m 


-Tt'Ui'f* 1 - Z a ‘.i 1 f + (TL) i- a M T i C i 


j=l 


3=1+1 


and 


1 < W < 2 and 1 «= i ^ m 


vivivz- 

j=l k=l 1=1 


diagonal term of i 
row in matrix [A] 


,th 


off-diagonal terms of 
a iJ P a j>i i^ 1 row in matrix [A] 


(TL)j = ^ ^ N^T^ + ^ ^li^l + ^i i^ 1 thermal load term 


k=l 


Z = 1 


This system of m equations is iteratively solved until a specified level of convergence 
is obtained. 

The boundary condition conductances and are calculated with the aid of the 

empirical relations listed in table I and discussed later. 


LIMITATIONS AND SCOPE 

This program has a present capacity of 400 volume elements, 200 boundary condition 
elements, 95 varying temperature boundary elements (see Computer Input in appendix A), 
and 15 different materials. Several minor limitations are discussed in the section on in- 
put where the parameters involved are introduced. The program is designed for an IBM 
7094 digital computer and is written in FORTRAN IV language. 

Provision is made for adding subroutines to calculate additional boundary condition 
coefficients without interrupting the sequence of statement numbers or altering present 
statements. Up to four forced convection, five free convection, and two varying tempera- 
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TABLE I. - BOUNDARY HEAT TRANSFER COEFFICIENTS 


Heat transfer type 
Forced convection 


Free convection 



Radiation 


Any surface 


ture conditions can be added by writing appropriate subroutines and by adding the required 
statements in subroutine CQN. The CQN statements required can be deduced from the 
portions labeled for present coefficients. Changes and some limitations are discussed in 
more detail in appendix B. 

The present coefficient equations are based on properties evaluated at the boundary 
film temperature (average of bulk fluid and surface temperatures). Unless the initial tem- 
perature estimates (required as a starting point) are very close to the final temperatures 
calculated, the input required will have to be changed and the problem rerun at least once 
to obtain the best results. 


CONCLUDING REMARKS 

The program, which is written in FORTRAN IV language for an IBM 7094 digital com- 
puter, is general and can be applied to a variety of axisymmetric problems. Various con- 
vection and radiation boundary conditions which can be used are given in the mathematical 
formulation. The program listing and flowcharts for steady- state thermal solution of an 
axisymmetric solid in cylindrical coordinates are given. Computer program application 
to a nosepiece assembly is given in detail. Additional demonstration of the program is 
provided in a comparison of the program solution to the exact solution for a finite cylin- 
der. This comparison shows a maximum error of 15 percent when a coarse mesh is 
used. Thus a partial check on program accuracy is provided. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 10, 1969, 

126-15. 
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APPENDIX A 


COMPUTER PROGRAM 
Preparation of Input Data 

The recommended first step is the preparation of the element pattern. For the seal 
analyzed in appendix C, a layout of the seal assembly five times actual size was used. 
There is no restriction on the numbering sequence except the use of consecutive numbers 
(i. e. , every number from 1 to P), but it is preferable to use an orderly system in order 
to avoid confusion and to expedite assembling the data required for input. The drawing 
also includes, for convenience, the radial dimension and mean radius of each element 
row and the axial dimension of each column. 

Adjacent elements should have the same common dimension. Slight variations will 
not cause appreciable errors if a contact resistance is added to compensate for the in- 
creased conduction path length of off-center nonmatching elements. 

Also, for convenience, the boundary conditions (pressure, temperature, fluid, and 
so forth, needed to determine the values of required properties) are included. The ma- 
terials of the various parts are included in a table (on the drawing if room is available). 

The items of data required are listed in the following section on computer input. A 
compact outline of most items can be obtained from the section Program Symbols and 
Definitions or from table II. The conduction element matrix (TD1) and boundary coeffi- 
cient lists (TD2) include input data. 


17 



TABLE H. - SUMMARY OF TD2 LISTS FOR SUBROUTINES TO CALCULATE COEFFICIENTS 


[Parameters in parentheses are calculated in subroutine. Except for TD2(20) of SN3 and VT2 and TD2(11) of FC1 and FC2, parameters after hA are 
merely to make writing FORTRAN expression for h easier. ] 


TD2 

Subroutine 



word 













or 

SN1 

SN2 

SN3 

SN4 

SN5 

FC1 

FC2 

FC3 

RA1 

1 SCI 

VTl 

VT2 

sub- 

(duct flow, 

(sides of 

(radial seal 

(concen- 

(flat 

(horizon- 

(vertical 

(liquid 

(radiaJ 

(specified 

(varying tempera- 

(varying tempera- 

script 

liquid) 

rotors) 

gap) 

trie cyl- 

plate) 

tal 

cylinder 

film 

tion) 

heat 

ture duct flow, 

ture radial seal gap) 





inders) 


cylinder) 

or plane) 

cooling, 


transfer 

liquid) 










gravity 


coeffi- 











flow) 


cient) 



1 

T b 

T b 

T b 

T b 

T b 

T b 

T b 

T b 

T e 

T b or T e 

T b 

T b 

2 

k t 

k f 

k f 

k f 



tf 

W/L 

f es 

h 

k f 


3 

c p,f 

c p,f 

C P,f 

c P,f 

c P,f 

c p,f 

S.f 

C P 

f se 

(hA) 

C„ f 

p, f 

c P,f 

4 

M f 

M f 


Mf 

H 

Mf 

M f 

d l 

e e 


M f 

pf 

5 

W a 

"f 

v f 

" f 

v f 

K i 

v f 

(h) 

e s 


W 

v i 

6 

D 

r av 

r av 

r o 

X 

D 

L 

(hA) 

(h) 


D 

r av 

7 

X 

s 

s 

r^ 

v_ 

(PrA 

(PrA 


(hA) 


X 

s 


10 


11 


12 




N /3 ) 


13 

(Re w ) 

(hf) 

14 

(Pr { ) 

(hA) 

15 

(h f ) 

<V r av> 

16 

(hA) 

(^a/' 

17 

(D/X) 

(K) 

18 

f (kf/D) 


19 

([(D/X)] 0 ' 7 ' 

> 
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(prl/3) 

(Re c ) 

(h f ) 

(hA) 

(Heat genera- 
tion factor) 


<Re v ) 

K /3 ) 

(hf) 

(hA) 

< V R> 

Kt) 

( V tot> 

(s) 


(Re x ) 


(Pr f ) 


(h) 


(Gr D ) 

(hf) 


(Gr L ) 

(hf) 

(hA) 


(hA) 

(hA) (Pr f Gr D ) (Pr f Gr D ) 
(AT) (AT) 


(1 - e e ) 



(T ) 

V 2LV' 


(T, 


(kj/D) (kj/L) 


Annular neighbor 
conduction number 
Flow neighbor con- 
duction number 

Combining flow 
indicator 
Flow initiation 
indicator 

A Y 


Combining flow 
indicator 
Flow initiation 
indicator 
Annular neighbor 



conduction number 

(Re w ) 

Flow neighbor con- 
duction number 

(Pr f ) 

W 

(h f ) 

< Re J 

(hA) 

(Pr f 1/3 ) 

(DA) 

(Re c ) 

(kj/D) 

(hf) 

([DA] 0 ' 7 ) 

(hA) 


(Heat generation 
factor) 


*Rotor or stator flag (zero for rotor, nonzero for stator). 


Computer Input 


The data input of a problem is submitted to the computer on the cards described in 
this section. Unless specifically stated, no input card may be omitted. 

Item 1: Title card (72H) 

word 1 columns 1 to 72 Title of problem being run 

Item 2: Output control (L. 7 ?) card (1016) for obtaining printouts in addition to that listed 
in the section Computer Output 

word 1 columns 1 to 6 A nonzero integer causes printout of the conduc- 

tion matrix (TD1 table) input 

word 2 columns 7 to 12 A nonzero integer causes printout of the bound- 

ary conductance matrix (H(NBPT, J), J=l, 9) 

word 3 columns 13 to 18 A nonzero integer causes printout of the element 

temperatures for the next to last iteration (un- 
less execution terminated by reaching maximum 
number of iterations) for comparison with tem- 
peratures of last iteration 

word 4 columns 19 to 24 Leave blank. This word is set within the pro- 

gram to cause printout of the convergence cri- 
terion (Item 5, word 4) and of the temperature 
change for each element at each N iteration; 
N is specified by word 2 of the Run data card 
(Item 4). 


word 5 

columns 25 to 30 

A nonzero integer causes printout of the bound- 
ary condition input matrix (TDA and TD2 lists) 

word 6 

columns 31 to 36 

Not used 

word 7 

columns 37 to 42 

Not used 

word 8 

columns 43 to 48 

A nonzero integer causes the final element tem- 
peratures to be punched on cards for use as ini 
tial estimates for another run 

word 9 

columns 49 to 54 

A nonzero integer causes thermal expansions to 


be punched on cards for use in other programs 
(e.g. , a stress analysis) 
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word 10 columns 55 to 60 A nonzero integer causes printout of element 

temperatures and heat fluxes after each itera- 
tion. Also thermal expansion data (a, T, a AT) 
are printed. This option is meant only for de- 
bugging data and should not be included for a 
standard run. 


Item 3: Miscellaneous options (NT) card (516) 


word 1 

columns 1 to 6 

Number of cases (designs) to be analyzed. If 
only one case, this word may be left blank. 

No other words now used; leave blank. 
Item 4: Run data (RD) card (5F6. 0) 

word 1 

columns 1 to 6 

Over- relaxation factor w (1 <w < 2) 

word 2 

columns 7 to 12 

Reference temperature for thermal expansion 
(a AT) 

word 3 columns 13 to 18 

No other words now used. 

Item 5: Case data (CD) card (10F6.0) 

The interval N required for word 4 of the LD 
card (Item 2) 

word 1 

columns 1 to 6 

Number of elements (400 maximum at present) 

word 2 

columns 7 to 12 

Number of boundary elements (200 maximum at 
present) 

word 3 

columns 13 to 18 

Maximum number of iterations allowed for a run 

word 4 

columns 19 to 24 

Temperature convergence criterion. The maxi- 
mum change in temperature of any element from 
one iteration to the next that is considered an 
acceptable error. 

word 5 

columns 25 to 30 

Not used 

word 6 

columns 31 to 36 

Starting temperature for varying temperature 
flow 1 (see varying temperature under boundary 
condition input (Item 10)) 

word 7 

columns 37 to 42 

Starting temperature for varying temperature 
flow 2 
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word 8 

columns 43 to 48 

Not used 

word 9 

columns 49 to 54 

To avoid numerical instability due to low flow 
rates in the varying temperature calculation, 
the boundary elements are subdivided. The 
number of subdivisions is specified between 1 
and 100 in this word. 

word 10 

columns 55 to 60 

Not used 


Item 6: Thermal conductivity data cards 

The thermal conductivity of each material is treated as a power series function of 

0 1 m 

temperature (k = AT + BT + . . . + NT ). Provision is made for curve fitting up to 
15 materials. The coefficients required are computed by an auxiliary program (or, of 
course, by hand) from available data. The degree of fit of a curve is defined as the high- 
est power of temperature appearing in the function. Thus, a linear fit is degree 1; a 
constant conductivity is degree 0; and so forth. 

The auxiliary program is designed to produce fits of all degrees up to 11 from which 
the best agreement with a set of average values of k is selected as the function to be 
used. Therefore, words 4 and 5 were included to avoid exceeding the computer capacity 
when the higher powers of temperature occurred. Later it was somewhat arbitrarily 
decided to limit the maximum degree of fit to five in view of data scatter due to variation 
between material lots. The primary reason was to increase the number of materials 
possible and to save computer storage space. 

Card 1 (5F10.0) 


word 1 
word 2 


word 3 
word 4 

word 5 


columns 1 to 10 
columns 11 to 20 


columns 21 to 30 
columns 31 to 40 

columns 41 to 50 


Upper temperature limit of curve fit 

Material code number. Each material must be 
assigned a number so that the correct proper- 
ties will be used in the program. 

Degree of fit (0 to 5) 

Thermal conductivity scaling factor. At present, 
this is 1. 

Temperature scaling factor. At present, this 
is 1. 


Card 2 (6E13.6) - The coefficients of the equation: A in columns 1 to 13, B in col- 
umns 14 to 26, and so forth. 
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Card 3 (blank card) - If less than 15 materials are used, a blank card must be inserted 
following the last conductivity card in order to terminate the reading of data. 

Item 7: Mean coefficient of linear thermal expansion cards 

The remarks made about conductivity (see Item 6) also apply to expansion. 

Card 1 (5F10.0) 

word 1 columns 1 to 10 Upper temperature limit of curve fit 

word 2 columns 11 to 20 Material code number 

word 3 columns 21 to 30 Degree of fit (0 to 5) 

word 4 columns 31 to 40 Thermal expansion scaling factor. At present, 

this is 1. 

word 5 columns 41 to 50 Temperature scaling factor. At present, this is 1. 

Card 2 (6E13.6) - The coefficients of the equation as in Item 6. 

Card 3 (blank card) - To terminate data read in for less than 15 materials (see Item 6). 

Item 8: Conduction matrix data (TD1) card (F4. 0, IX, F7. 4, 2F5. 3, IX, 4F4. 0, F3. 0, 
4F7.0) - 1 card per element; element numbers in numerical order 


word 1 

columns 1 to 4 

Element conduction number 


word 2 

columns 6 to 12 

Mean element radius, r ay , in. 


word 3 

columns 13 to 17 

Element radial width, Ar, in. 


word 4 

columns 18 to 22 

Element axial length, A l , in. 


word 5 

columns 24 to 27 

East neighbor conduction number, jj 

word 6 

columns 28 to 31 

North neighbor conduction number, jg 

word 7 

columns 32 to 35 

West neighbor conduction number, jg 

word 8 

columns 36 to 39 

South neighbor conduction number, 

word 9 

columns 40 to 42 

Material code number (see Item 6, word 2) 

word 10 

columns 43 to 49 

East contact resistance, R , 

c l,i 

North contact resistance, R , 

c 9 i 

(sec)(°F)/Btu 

word 11 

columns 50 to 56 

(sec)(°F)/Btu 

word 12 

columns 57 to 63 

West contact resistance, R , 

Co 4 

(sec)(°F)/Btu 

word 13 

columns 64 to 70 

O, 1 

South contact resistance, R , 

^ A i 

(sec)(°F)/Btu 
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Item 9: Initial temperature cards (16F5.0) - Put 16 initial element temperatures per 
card. 

words 1 columns 1 to 5, 
to 16 6 to 10, 

11 to 15, 
and so forth 

Item 10: Boundary heat transfer cards 

Each boundary element requires an indicator card. Following the indicator card is a 
set of one or more cards for each boundary condition associated with the element. A par- 
tial outline is presented in table n. This sequence is repeated for each boundary element. 

Card 1: Boundary indicator card (316) - 1 card per element 

word 1 columns 1 to 16 Boundary number of the element. In addition to 

the comduction numbers which are assigned to 
the elements, each boundary element is given a 
boundary number. Any convenient order of 
numbering can be used. This is for convenience 
in preparing data. 

word 2 columns 7 to 12 Conduction number of the boundary element. 

Most computation is based on the conduction 
number. 

word 3 columns 13 to 18 Boundary heat transfer type code number for 

first basic type associated with the element. 
Code numbers are as follows: 

(1) Specified heat flux - the net flux for all 
faces for which it is specified. If this type 
is associated with a node, it must appear 
first; but the other types need not appear in 
any fixed order. 

(2) Forced convection 

(3) Free convection 

(4) Radiation 

(5) Specified heat transfer coefficient 

(6) Specified temperature - By definition this 
will be the only condition for an element; 
the program considers only its effect on 
adjacent elements. 


Initial element temperatures (°F) are listed in 
sequence from 1 to the total number of ele- 
ments: elements 1 to 16 on the first card, ele- 
ments 17 to 32 on the second, and so forth. 
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Card 2: Index card (for each condition: basic type or subtype) 

This card has two forms. Form A is used for specified parameters (types 1 and 7) 
which require no additional cards. Form B (TDA) is used for the remaining types which 
do require additional cards. 

Form A for specified heat flux (15, E13. 6, 16) 

word 1 columns 1 to 5 Conduction number of element 

word 2 columns 6 to 18 Heat flux, Btu/sec 

word 3 columns 19 to 24 A nonzero integer if other types or subtypes 

associated with the node follow; blank (or 0) if 
no types or subtypes follow. 

Form A for specified temperature (15, IX, F6. 0) 

word 1 columns 1 to 5 Conduction number of element 

word 2 columns 7 to 12 Temperature of element, °F 

Form B (6F6. 0, F12. 0) - These (TDA) cards are the same for all types using them. 

word 1 columns 1 to 6 Conduction number of element 

word 2 columns 7 to 12 Heat transfer type (2 to 6) 

word 3 columns 13 to 18 Heat transfer subtype code number. (For additions 

or substitutions see limitations and scope. ) 

Subtype code numbers for type 2 (forced convec- 
tion) are as follows: 

(1) Duct flow, liquid - for explanation see 
headings under additional cards 

(2) Sides of rotors 

(3) Radial seal gap 

(4) Concentric cylinders 

(5) Flat plate 

(6) to (9) for additional subtypes 

Subtype code numbers for type 3 (free convection) 
are as follows: 

(1) Horizontal cylinder 

(2) Vertical cylinder or plane 

(3) Liquid film cooling, gravity flow 

(4) to (8) for additional subtypes 
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Subtype code numbers for type 6 (varying temper- 
ature forced convection) are as follows: 

(1) Duct flow, liquid 

(2) Radial seal gap 

(3) to (4) for additional subtypes 

Types 4 and 5, at present, have no subtypes; so 
word 3=0. 

word 4 columns 19 to 24 Nonzero number if other types or subtypes follow 

word 5 columns 25 to 30 Number of items to be read from the remaining 

cards for this condition 

word 6 columns 31 to 36 Area indicator code number specifies to which 

face of the element the condition applies; these 
numbers are as follows: 

(1) East face 

(2) North face 

(3) West face 

(4) South face 

(5) Specified area - used for special cases such 
as a small duct passing through an element 
or a boundary at an appreciably nonparallel 
or nonnormal angle to the axis. Since con- 
ductive resistance is assumed to be negligi- 
ble compared with boundary film resistance, 
the problem of location does not occur. 

word 7 columns 37 to 48 If a specified area (code number 5 in word 6) is 

used, the value in square feet is placed here; 
otherwise, this word is left blank. 

For types 2 to 6, the following additional cards are needed as shown under types 4 
and 5 and the subtypes of types 2, 3, and 6. Card 1 is the indicator card described pre- 
viously. Card 2 is the TDA card for the type or subtype. The properties are evaluated 
at the tem perature indicated by the equations in table I. 

Steady-state forced convection (type 2): 

Subtype 1 - Duct flow, liquid 

Card 3 (7F11.0) 

word 1 columns 1 to 11 Fluid bulk temperature, °F 
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word 2 

columns 12 to 22 

Thermal conductivity, Btu/(ft)(sec)(°F) 

word 3 

columns 23 to 33 

Specific heat, Btu/(lbm)(°F) 

word 4 

columns 34 to 44 

Absolute viscosity, lbm/(ft)(sec) 

word 5 

columns 45 to 55 

Mass flow per square foot of cross section, 
lbm/(ft^)(sec) 

word 6 

columns 56 to 66 

Hydraulic diameter, 4A /P , ft 

A W 

word 7 

columns 67 to 77 

Flow length (distance from start of duct to center 
of element), ft 

Subtype 2 - Sides 

of rotors (see Use 

of Local Coefficients in appendix D) 

Card 3 (7F11.0) 

word 1 

columns 1 to 11 

Fluid bulk temperature, °F 

word 2 

columns 12 to 22 

Thermal conductivity, Btu/(ft)(sec)(°F) 

word 3 

columns 23 to 33 

Specific heat, Btu/(lbm)(°F) 

word 4 

columns 34 to 44 

Absolute viscosity, lbm/(ft)(sec) 

word 5 

columns 45 to 55 

Kinematic viscosity, ft /sec 

word 6 

columns 56 to 66 

Average radius of element, ft 

word 7 

columns 67 to 77 

Gap between rotor and stator, ft 

Card 4 (3F11.0) 

word 1 

columns 1 to 11 

Rotor angular velocity, rad/sec 

word 2 

columns 12 to 22 

Outer radius of rotor or stator, ft 

word 3 columns 23 to 33 

Subtype 3 - Radial seal gap 

Rotor or stator flag (zero for rotor, nonzero for 
stator) 

Card 3 (7F11.0) 

word 1 

columns 1 to 11 

Fluid bulk temperature, °F 

word 2 

columns 12 to 22 

Thermal conductivity, Btu/(ft)(sec)(°F) 

word 3 

columns 23 to 33 

Specific heat, Btu/(lbm)(°F) 

word 4 

columns 34 to 44 

Absolute viscosity, lbm/(ft)(sec) 


26 



word 5 

columns 45 to 55 

o 

Kinematic viscosity, ft /sec 

word 6 

columns 56 to 66 

Average radius of element, ft 

word 7 

columns 67 to 77 

Gap between rotor and stator, ft 

Card 4 (2F11.0) 

word 1 

columns 1 to 11 

Rotor angular velocity, rad/sec 

word 2 

columns 12 to 22 

Outer radius of rotor or stator, ft 

Subtype 4 - Concentric cylinders (inner 

cylinder rotating) 

Card 3 (7F11.0) 

word 1 

columns 1 to 11 

Fluid bulk temperature, °F 

word 2 

columns 12 to 22 

Thermal conductivity, Btu/(ft)(sec)(°F) 

word 3 

columns 23 to 33 

Specific heat, Btu/(lbm)(°F) 

word 4 

columns 34 to 44 

Absolute viscosity, lbm/(ft)(sec) 

word 5 

columns 45 to 55 

Kinematic viscosity, ft /sec 

word 6 

columns 56 to 66 

Outer radius of annulus between cylinders, ft 

word 7 

columns 67 to 77 

Inner radius of annulus, ft 

Card 4 (2F11.0) 

word 1 

columns 1 to 11 

Rotational velocity of inner cylinder, rad/sec 

word 2 

columns 12 to 22 

Axial flow velocity, ft/sec 

Subtype 5 - Flat plate 


Card 3 (7F11.0) 

word 1 

columns 1 to 11 

Fluid bulk temperature, °F 

word 2 

columns 12 to 22 

Thermal conductivity, Btu/(ft)(sec)(°F) 

word 3 

columns 23 to 33 

Specific heat, Btu/(lbm)(°F) 

word 4 

columns 34 to 44 

Absolute viscosity, lbm/(ft)(sec) 

word 5 

columns 45 to 55 

2 

Kinematic viscosity, ft /sec 

word 6 

columns 56 to 66 

Flow length (distance from start of plate to cen- 
ter of element), ft 

word 7 

columns 67 to 77 

Velocity of bulk fluid, ft/sec 
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Free convection (type 3): 


Subtype 1 (horizontal cylinder) or 2 (vertical cylinder or plane), gas - identical except 
for the characteristic length used 


Card 3 (6F11.0) 



word 1 

columns 1 to 11 

Fluid bulk temperature, °F 

word 2 

columns 12 to 22 

Thermal conductivity, Btu/(ft)(sec)(°F) 

word 3 

columns 23 to 33 

Specific heat, Btu/(lbm)(°F) 

word 4 

columns 34 to 44 

Absolute viscosity, lbm/(ft)(sec) 

word 5 

columns 45 to 55 

9 

Kinematic viscosity, ft /sec 

word 6 

columns 56 to 66 

Subtype 1: hydraulic diameter, ft 


Subtype 2: flow length (distance from start of 
flow to top of surface of which element is a 
part), ft 

Subtype 3 - Liquid film cooling, gravity flow 


Card 3 (4F11.0) 



word 1 

columns 1 to 11 

Fluid film temperature, °F 

word 2 

columns 12 to 22 

Mass flow per unit length (at a right angle to flow 



direction), lbm/(ft)(sec) 

word 3 

columns 23 to 33 

Specific heat, Btu/(lbm)(°F) 

word 4 

columns 34 to 44 

Flow length (distance from point where flow 


starts to where it leaves surface (e.g. , one- 
half circumference of horizontal cylinder)), ft 


Graybody radiation (type 4): 


Card 3 (4F11.0) 

word 1 
word 2 
word 3 
word 4 
word 5 


columns 1 to 11 
columns 12 to 22 
columns 23 to 33 
columns 34 to 44 
columns 45 to 55 


Environmental temperature, °F 
View factor, environment to surface 
View factor, surface to environment 
Emissivity of environment 
Emissivity of surface 
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Specified heat transfer coefficient (type 5): 

Card 3 (2F11.0) 

word 1 columns 1 to 11 Fluid bulk (or environmental) temperature, °F 

word 2 columns 12 to 22 Heat transfer coefficient, Btu/(ft^)(sec)(°F) 

Varying temperature forced convection (type 6): 

Subtype 1 - Varying temperature duct flow, liquid 


Card 3 (7F11.0) 

word 1 

columns 

1 to 11 

Fluid bulk temperature, °F 

word 2 

columns 

12 to 22 

Thermal conductivity, Btu/(ft)(sec)(°F) 

word 3 

columns 

23 to 33 

Specific heat, Btu/(lbm)(°F) 

word 4 

columns 

34 to 44 

Absolute viscosity, lbm/(ft)(sec) 

word 5 

columns 

45 to 55 

Mass flow per square foot of cross section, 
lbm/(ft^)(sec) 

word 6 

columns 

56 to 66 

Hydraulic diameter, 4A X /P W , ft 

word 7 

columns 

67 to 77 

Flow length (distance from start of duct to center 
of element), ft 

Card 4 (4F11.0) 

word 1 

columns 

1 to 11 

Conduction number of annular neighbor (element 
on other boundary of flow if flow is between two 
elements). If none, leave blank. 

word 2 

columns 

12 to 22 

Conduction number of flow neighbor (element up- 
stream of element under consideration). If this 
is the first element along the flow, leave blank. 

word 3 

columns 

23 to 33 

Combined flow indicator (see remarks at end of 
this subtype) 

word 4 

columns 

34 to 44 

Flow initiation indicator (see remarks at end of 
this subtype) 

word 5 

columns 

45 to 55 

2 

Cross section of duct, ft 


There is provision in the program for varying temperature flows of two starting tem- 
peratures. Any number of varying temperature flows can be handled subject to the fol- 
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lowing additional restrictions. The program will handle only one such boundary condition 
per element. The number of elements with varying temperature flows cannot exceed half 
the number of boundary elements (95) at present. There can be no more than five com- 
bined flows (explained later). 

The flow initiation indicator distinguishes the two starting temperatures. If 1 is 
used, the starting temperature is word 1, Card 2 of Item 5 (Case data). If 2 is used, the 
starting temperature is word 2, Card 2 of Item 5. 

If two flows (limit of program) combine at a boundary element to form a new flow, 
this element is a combined flow element. For a combined flow element, the flow initia- 
tion indicator must be zero. Five such combinations are permitted by the program. An 
additional card is then required. 

Card 5 (216) - only for combined flow elements 


word 1 

columns 1 to 6 

Conduction number of one upstream neighbor ele 
ment 

word 2 

columns 7 to 12 

Conduction number of other upstream neighbor 
element 

Subtypes - Varying temperature radial seal gap 

Card 1 (7F11.0) 

word 1 

columns 1 to 11 

Fluid bulk temperature, °F 

word 2 

columns 12 to 22 

Thermal conductivity, Btu/(ft)(sec)(°F) 

word 3 

columns 23 to 33 

Specific heat, Btu/(lbm)(°F) 

word 4 

columns 34 to 44 

Absolute viscosity, lbm/(ft)(sec) 

word 5 

columns 45 to 55 

9 

Kinematic viscosity, ft /sec 

word 6 

columns 56 to 66 

Average radius of element, ft 

word 7 

columns 67 to 77 

Gap between rotor and stator, ft 

Card 2 (7F11.0) 

word 1 

columns 1 to 11 

Rotor angular velocity, rad/sec 

word 2 

columns 12 to 22 

Outer radius of rotor, ft 

word 3 

columns 23 to 33 

Combined flow indicator 

word 4 

columns 34 to 44 

Flow initiation indicator 

word 5 

columns 45 to 55 

Annular neighbor conduction number 
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I 

columns 56 to 66 Flow neighbor conduction number 

columns 67 to 77 Mass flow, lbm/sec 

only for combined flow elements 

columns 1 to 6 Conduction number of one upstream neighbor 

element 

columns 7 to 12 Conduction number of other upstream neighbor 

element 

Computer Output 

Standard output of results . - The standard output of the program is obtained with all 
list options (Item 2 of Computer Input) set at zero with the possible exception of the 
punch options (words 8 and 9). This output is as follows: 

Item 1: Problem title 

Item 2: Printout options (see input Item 2) 

Item 3: Miscellaneous options (see input Item 3) 

Item 4: Run data (see input Item 4) 

Item 5: Case data (see input Item 5) 

Item 6: Thermal conductivity curve fit coefficients for each material 

Item 7: Thermal expansion curve fit coefficients for each material 

Item 8: Element starting temperature estimates 

Item 9: Varying temperature boundary element dimensionless coefficient, hA/WCp 

Item 10: Varying temperature boundary element fluid internal heat generation 
Item 11: Conduction number and preceding element conduction numbers for each 
combined flow element 

Item 12: Conduction number, varying temperature number, and flow code for vary- 
ing temperature boundary elements 
Item 13: Element steady-state temperatures 

Item 14: Boundary element forced convection, free convection, radiation, and vary- 
ing temperature convection conductances 
Item 15: Varying temperature boundary condition fluid temperatures 
Item 16: Element thermal conductivities 
Item 17: Element heat fluxes to neighbor elements 
Item 18: Element conductances to neighbor elements 

Item 19: Element boundary conductance sums and boundary heat flux sums 
Item 20: Element free thermal expansion, a AT 


word 6 
word 7 

Card 5 (216) - 
word 1 

word 2 
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An illustration of the output can be found in the listing for the sample problem in 
appendix C. 

Program error m essages. - The program includes checks of some of the input to 
ensure correct relations. First, the conductivity (TD1) cards are checked for numerical 
order of conduction numbers. Once this has been achieved, a check is made for neighbor 
agreement. For example, if element 10 has element 11 listed as east neighbor, ele- 
ment 11 must have element 10 as west neighbor, assuming this is the correct relation. 
When the neighbor agreement is complete, a check is made of boundary index cards to 
ensure numerical order of boundary numbers. Then a check is made to see that the con- 
duction number (s) on the element boundary index card(s) agree with the conduction num- 
ber on the indicator card. Also there is a check for agreement of the boundary number 
corresponding to an element conduction number with the conduction number correspond- 
ing to the element boundary number. 


AFN(NTF) 
ALPHA(15, 11) 
CE(NJPTS, 4) 
GFN(NTF) 
H(NBPTS, 9) 
JBNT(JBPT) 
JBP 
JBPT 

JNTF(JVPT) 

JPT 

JVPT 

MATL(15, 11) 

NAT(NDR) 

NBDIM 

NBND(NBP) 

NBP 

NBPTS 


Program Symbols and Definitions 

conduction number of varying temperature element annular neighbor 
thermal expansion curve fit matrix (input Item 6) 
neighbor conductance matrix 

conduction number of varying temperature element flow neighbor 

boundary conductance matrix 

boundary number of element JBPT 

conduction number of boundary element (in card check) 

conduction number of boundary element, NBND(NBP) 

varying temperature number of element JVPT 

conduction number of element 

conduction number of varying temperature element, NODE(NTF) 
thermal conductivity curve fit matrix (input Item 5) 
conduction number of one adjacent element upstream of NDR 
maximum number of boundary elements accepted by program 
conduction number of boundary element NBP 
boundary number of boundary element 
actual number of boundary elements 
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NBR(NTR) 


flow number (input Item 10, Card 4, word 3) of varying temperature 
element 

NBT(NDR) conduction number of other adjacent element upstream of NDR 

NDIM maximum number of conduction elements accepted by program 

NDR combination number of combined flow varying temperature element 

NITS maximum number of iterations allowed for analysis 

NJPTS actual number of conduction elements 

NM(NDR) conduction number of combined flow element NDR 

NODE(NTF) conduction number of varying temperature element NTF 

NTF varying temperature number of varying temperature boundary ele- 

ment 

NTR noncombining flow number of varying temperature element 

NVDIM maximum number of varying temperature elements accepted 

QF(NTF) fluid shear heat generation at varying temperature element NTF 

Q(JBPT) sum of boundary heat fluxes of conduction element JBPT 

SUMH(NBPT) sum of boundary conductances of boundary element NBPT 

SUMT(NBPT) sum of boundary heat fluxes of boundary element NBPT 

TBND(NTF) fluid temperature at varying temperature boundary element 

TCF(NTF) dimensionless coefficient of varying temperature element, kA/WC 

r 

TDA(7) index card (input Item 10, Card 2) list 

TD1(NJPTS, 19) conduction data matrix 

TD2(20) properties cards (input Item 10, Card(s) 3 and those following) list 

THEXP(JPT) element thermal expansion, a AT 

T(JPT) initial temperature estimate for element JPT 

TS(INC) fluid temperature calculation by subdivision of boundary element 

Conductance matrix CE(NJPTS,4) . - The symbols for conductance matrix 
CE(NJPTS,4) are as follows: 

CE(JPT, 1) conductance between conduction element JPT and east neighbor 

CE(JPT,2) conductance between conduction element JPT and north neighbor 

CE(JPT, 3) conductance between conduction element JPT and west neighbor 

CE(JPT,4) conductance between conduction element JPT and south neighbor 
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Boundary conductance matrix H(NBPTS, 9). - The symbols for boundary conduc- 
tance matrix H(NBPTS, 9) are as follows: 

H(NBPT, 1) specified heat flux 

H(NBPT, 2) sum of forced convection (including specific coefficient) conductances 

H(NBPT, 3) sum of forced convection conductance x fluid temperature 

H(NBPT, 4) sum of free convection conductances 

H(NBPT, 5) sum of free convection conductance x fluid temperature 

H(NBPT,6) sum of radiation "conductances” 

H(NBPT, 7) sum of radiation "conductance" x environmental temperature 

H(NBPT, 8) sum of varying temperature conductances 

H(NBPT, 9) sum of varying temperature conductance x fluid temperature 

Conduction data matri x TD1(NJPTS, 19). - The symbols for conduction data matrix 
TD1(NJPTS, 19) are as follows: 

TD1(JPT, 1) element conduction number 

TD1(JPT,2) mean element radius 

TD1(JPT, 3) element radial width, Ar 

TD1(JPT, 4) element axial length, A l 

TD1(JPT, 5) east neighbor conduction number, 

TD1(JPT, 6) north neighbor conduction number, jg 

TD1(JPT,7) west neighbor conduction number, jg 

TD1(JPT, 8) south neighbor conduction number, 

TD1(JPT, 9) material code number 

TD1(JPT, 10) east contact resistance, R 

c li 

TD1(JPT, 11) north contact resistance, R 

c 2i 

TD1(JPT, 12) west contact resistance, R 

c 3i 

TD1(JPT, 13) south contact resistance, R c 

TD1(JPT, 14) east and west boundary areas (added during run) 

TD1(JPT, 15) north boundary area (added during run) 

TD1(JPT, 16) south boundary area (added during run) 

TD1(JPT, 17) boundary condition type code number (added during run) 
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TD1(JPT,18) thermal conductivity (added during run) 

TD1(JPT, 19) element temperature (added during run) 

TD2 list for subroutine SN1, type 2 (forced conve ction), subtype 1 (duct flow, 
liquid) . - The TD2 list is as follows: 

TD2(1) fluid bulk temperature, 

TD2(2) boundary film thermal conductivity, kj 

TD2(3) boundary film specific heat, C f 

TD2(4) boundary film absolute viscosity, p f 

TD2(5) fluid mass flow per square foot of cross section, W 

cL 

TD2(6) duct hydraulic diameter, D 

TD2(7) length from start of flow to center of element,' X 

TD2(8) not used in order to obtain compatability with VT1 

TD2(9) not used in order to obtain compatability with VT1 

TD2(10) not used in order to obtain compatability with VT1 

TD2(11) not used in order to obtain compatability with VT1 

TD2(12) not used in order to obtain compatability with VT1 

TD2(13) Reynolds number (calculated in subroutine), Re w 
TD2(14) Prandtl number (calculated in subroutine), Pr^ 

TD2(15) heat transfer coefficient (calculated in surboutine), hj 
TD2(16) conductance across boundary film (calculated in subroutine), hA 
TD2(17) hydraulic diameter/flow length (calculated in subroutine), D/X 
TD2(18) conductivity/hydraulic diameter (calculated in subroutine), k^/D 
TD2(19) 0. 7 power of TD2(16) (calculated in subroutine), (D/X)°* 7 

TD2 list for subroutine SN2, type 2 (forced convec tion), subtype 2 (sides of rotors). - 
The TD2 list is as follows: 

TD2(1) fluid bulk temperature, T^ 

TD2(2) boundary film thermal conductivity, kj 

TD2(3) boundary film specific heat, C f 

TD2(4) boundary film absolute viscosity, 

TD2(5) boundary film kinematic viscosity, 
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TD2(6) average radius of element, r &y 

TD2(7) clearance between rotor and stator, s 

TD2(8) rotor angular velocity, a> 

TD2(9) outer radius of rotor, a 

TD2(10) rotor or stator flag (zero for rotor, nonzero for stator) 

TD2(11) Reynolds number (calculated in subroutine), Re^ 

TD2(12) cube root of Prandtl number (calculated in subroutine), Prjf^ 

TD2(13) heat transfer coefficient (calculated in subroutine), 

TD2(14) conductance across boundary film (calculated in subroutine), hA 
TD2(15) conductivity/ average radius (calculated in subroutine), k^/ r ay 

TD2(16) 0.1 power of clearance/ average radius (calculated in subroutine), 

( s / r av>°‘ 1 

TD2(17) fluid angular velocity/ rotor angular velocity (calculated in subroutine), K 

TD2 list for subroutine SN 3, type 2 (forced convection), subtype 3 (radial seal gap). - 
The TD2 list is as follows : 

TD2(1) fluid bulk temperature, 

TD2(2) boundary film thermal conductivity, kj 

TD2(3) boundary film specific heat, C . 

P? 1 

TD2(4) boundary film absolute viscosity, p, 

TD2(5) boundary film kinematic viscosity, 

TD2(6) average radius of element, r„„ 

TD2(7) clearance between rotor and stator, s 

TD2(8) rotor angular velocity, o> 

TD2(9) outer radius of rotor, a 

TD2(10) not used to obtain compatability with VT2 

TD2(11) not used to obtain compatability with VT2 

TD2(12) not used to obtain compatability with VT2 

TD2(13) not used to obtain compatability with VT2 

TD2(14) not used to obtain compatability with VT2 

TD2(15) Reynolds number (calculated in subroutine), Re^ 
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1/3 

TD2(16) cube root of PrandtL number (calculated in subroutine), Pr^/ 

TD2(17) critical Reynolds number (calculated in subroutine), Re c 
TD2(18) heat transfer coefficient (calculated in subroutine), h f 

TD2(19) conductance from center of seal gap to surface (calculated in subroutine), hA 

TD2(20) heat generation factor (calculated in subroutine) 

TD2 list for subroutine SN4, type 2 (forced convection), subt ype 4 (concentric cylin- 
ders) . - The TD2 list is as follows: 

TD2(1) fluid bulk temperature, 

TD2(2) boundary film thermal conductivity, k f 

TD2(3) boundary film specific heat, C f 

TD2(4) boundary film absolute viscosity, |U £ 

TD2(5) boundary film kinematic viscosity, 

TD2(6) outer radius of annulus between cylinders, r Q 

TD2(7) inner radius of annulus between cylinders, r £ 

TD2(8) angular velocity of inner rotating cylinder, c o 

TD2(9) axial flow velocity, V a 

TD2(10) Reynolds number (calculated in subroutine), Re y 

1/3 

TD2(11) cube root of Prandtl number (calculated in subroutine), .Pr^. 

TD2(12) heat transfer coefficient (calculated in subroutine), h £ 

TD2(13) conductance across boundary film (calculated in subroutine), hA 

TD2(14) linear velocity of inner cylinder (calculated in subroutine), V R 

2 

TD2(15) square of total fluid velocity (calculated in subroutine), V tot 
TD2(16) total fluid velocity (calculated in subroutine), V tot 

TD2(17) annular clearance (calculated in subroutine), s 

TD2 list for subroutine SN5, type 2 (forced convection), subtype 5 (flat plate). - The 
TD2 list is as follows: 

TD2(1) fluid bulk temperature, T fe 

TD2(2) boundary film thermal conductivity, k f 

TD2(3) boundary film specific heat, C f 

TD2(4) boundary film absolute viscosity, 
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TD2(5) boundary film kinematic viscosity, 

TD2(6) length from start of flow to center of element, X 

TD2(7) bulk fluid velocity, 

TD2(8) Reynolds number (calculated in subroutine), Re x 

TD2(9) Prandtl number (calculated in subroutine), Pr^ 

TD2(10) heat transfer coefficient (calculated in subroutine), h 

TD2(11) conductance across boundary film (calculated in subroutine), hA 

TD2 list for subroutine FC1 , type 3 (free convection), subtype 1 (horizontal cylin- 
der). - The TD2 list is as follows: 

TD2(1) fluid bulk temperature, 

TD2(2) boundary film thermal conductivity, k f 

TD2(3) boundary film specific heat, C f 

TD2(4) boundary film absolute viscosity, /u f 

TD2(5) boundary film kinematic viscosity, 

TD2(6) diameter of cylinder, D 

TD2(7) Prandtl number (calculated in subroutine), Pr^ 

TD2(8) Grashof number (calculated in subroutine), Gr D 

TD2(9) heat transfer coefficient (calculated in subroutine), h^ 

TD2(10) conductance across boundary film (calculated in subroutine), hA 

TD2(11) Prandtl number x Grashof number (calculated in subroutine), Pr f Gr D 

TD2(12) fluid to surface temperature difference (calculated in subroutine), AT 

TD2(13) average (or film) temperature (calculated in subroutine), T &y 

TD2(14) conductivity/ cylinder outside diameter (calculated in subroutine), kj/D 

TD2 list for subroutine FC2, t ype 3 (free convection), subtype 2 (vertical cylinder 
or plane) . - The TD2 list is as follows: 

TD2(1) fluid bulk temperature, T^ 

TD2(2) boundary film thermal conductivity, 

TD2(3) boundary film specific heat, C f 

TD2(4) boundary film absolute viscosity, 

TD2(5) boundary film kinematic viscosity, 
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TD2(6) total height of surface, L 

TD2(7) Prandtl number (calculated in subroutine), Pr f 

TD2(8) Grashof number (calculated in subroutine), Gr^ 

TD2(9) heat transfer coefficient (calculated in subroutine), h f 

TD2(10) conductance across boundary film (calculated in subroutine), hA 

TD2(11) Prandtl number x Grashof number (calculated in subroutine), Pr^Gr^ 

TD2(12) fluid to surface temperature difference (calculated in subroutine), AT 

TD2(13) average temperature (calculated in subroutine), T ay 

TD2(14) conductivity/surface height (calculated in subroutine), kj/L 

T D2 list for subroutine FC3, type 3 (free convection), subtype 3 (liquid film cooling, 
gravity flow). - The TD2 list is as follows: 

TD2(1) fluid bulk temperature, T^ 

TD2(2) mass flow of liquid per unit length normal to flow, W/L 

TD2(3) fluid specific heat, C 

tr 

TD2(4) flow length (e.g., one-half circumference of horizontal cylinder), 

TD2(5) heat transfer coefficient (calculated in subroutine), h 

TD2(6) conductance from surface to liquid (calculated in subroutine), hA 

TD2 list for subroutine RA1, type 4 (radiation) . - The TD2 list is as follows: 

TD2(1) environmental temperature, T g 

TD2(2) view factor from environment to element surface, f 

6S 

TD2(3) view factor from element surface to environment, f 

s e 

TD2(4) emissivity of environment, e g 

TD2(5) emissivity of element surface, e_ 

TD2(6) radiation coefficient (calculated in subroutine), h 

TD2(7) "conductance” from element surface to environment (calculated in sub- 

routine), hA 

TD2(8) reflectivity of environment (calculated in subroutine), 1 - e g 

TD2(9) reflectivity of element surface (calculated in subroutine), 1 - e 0 

s 

TD2(10) difference between fourth powers of temperatures/ temperature difference 
(calculated in subroutine), (t^ - T g)/( T e “ T s ) 
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TD2 list for subroutine SCI, type 5 (specified heat transfer coefficient). - The TD2 
list is as follows: 

TD2(1) fluid bulk (or environmental) temperature, (or T g ) 

TD2(2) specified heat transfer coefficient, h 

TD2(3) conductance between surface and bulk fluid (calculated in main program), hA 

TD2 list for subroutine VT1, type 6 (varying temperature forced convection), sub- 
type 1 (duct flow, liq uid). - The TD2 list is as follows: 

TD2(1) fluid bulk temperature, T fe 

TD2(2) boundary film thermal conductivity, k^ 

TD2(3) boundary film specific heat, C f 

p, i 

TD2(4) boundary film absolute viscosity, 

TD2(5) fluid mass flow per square foot of cross section, W & 

TD2(6) duct hydraulic diameter, D 

TD2(7) length from start of flow to center of element, X 

TD2(8) annular neighbor conduction number 

TD2(9) flow neighbor conduction number 

TD2(10) combining flow indicator code number 

TD2(11) flow initiation indicator code number (must be zero for combined flow ele- 
ment) 

TD2(12) duct cross section, A 

TD2(13) Reynolds number (calculated in subroutine), Re w 

TD2(14) Prandtl number (calculated in subroutine), Pr^ 

TD2(15) heat transfer coefficient (calculated in subroutine), hj 

TD2(16) conductance across boundary film (calculated in subroutine), hA 

TD2(17) hydraulic diameter/flow length (calculated in subroutine), D/X 

TD2(18) conductivity /hydraulic diameter (calculated in subroutine), kj/D 

TD2(19) 0.7 power of TD2(16) (calculated in subroutine), (D/X)®’ ^ 

TD2 list for subroutine VT2, type 6 (varying temperature forced convection), sub- 
type 2 (radial seal gap) . - The TD2 list is as follows: 

TD2(1) fluid bulk temperature, T fe 

TD2(2) boundary film thermal conductivity, k f 
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TD2(3) 

TD2(4) 

TD2(5) 

TD2(6) 

TD2(7) 

TD2(8) 

TD2(9) 

TD2(10) 

TD2(11) 

TD2(12) 

TD2(13) 

TD2(14) 

TD2(15) 

TD2(16) 

TD2(17) 

TD2(18) 

TD2(19) 

TD2(20) 


boundary film specific heat, C f 
boundary film absolute viscosity, fXj 
boundary film kinematic viscosity, 
average radius of element, r &y 
clearance between rotor and stator, s 
rotor angular velocity, o> 
outer radius of rotor, a 
combining flow indicator code number 

flow initiation indicator code number (must be zero for combined flow ele- 
ment) 

annular neighbor conduction number 
flow neighbor conduction number 
fluid mass flow rate, W 

Reynolds number (calculated in subroutine), Re^ 

1/3 

cube root of Prandtl number (calculated in subroutine), Pr f ' 
critical Reynolds number (calculated in subroutine), Re c 
heat transfer coefficient (calculated in subroutine), hj 

conductance from center of seal gap to surface (calculated in subroutine), hA 
heat generation factor (calculated in subroutine) 


Flow Charts and Program Listing 

The flow charts for the main routine (SEAL2) and subroutines CQN and QT are pre- 
sented in figures 8, 9, and 10. A listing of the program is as follows: 
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n n o n o 


$IBFTC S t AL 2 


STEADY STATE HEAT TRANSFER ANALYSIS OF AXI SYMMETRl C SYSTEMS TO 
INCLUDE VARYING FLUID TEMPERATURES ALONG BOUNDARIES APPLIED TO 
ADVANCED AIR BREATHING ENGINE SEAL DESIGN 

DIMENSION Tul(400,L9) ,T<400) ,CE(400,4) ,Q(403) , THEX P (400 ) , T 02( 20 ) 
DIMENSION H( 20 0,9} , SJMHI200) ,SUMT( 200 ) ,NBND { 2 3 0 I , T BND( 9 5 I 
DIMENSION L0( 10) ,NI (5) ,RJ(5) ,C0(10) ,TDA(7) ,MATL(i5, 11), NODE! 95) 
DIMENSION AL 5 HA ( 15, 11) ,JNTF (400) 

COMMON/B^K 1/J3nT< 400) , KA N (200) , NTF , NOR » NTR , QF ( 9 5 ) , NQR 

COMMON /B.K2/TCF( 95) , G FN l 95 ) , AF N { 95 ) , NBR ( 9 5 > , NM ( 5 ) , NAT ( 5 ) ,NBT(5 ) 

REAL MATL 

101 FOR IMA T ( 1H1) 

10 2 FORMAT (F4.0, IX,F7.4,2F5. 3,LX,4F4.J,F3.0 , F 7 . 3 ) 

103 FORMAT ( 1H0, 10X, 43HC0N3UC TIGN (TD1) CARO ORDER ERROR, CARD NO., 13, 2 
2X , BHwRONG ) 

104 FORMAT (1H0, 10X,23HNAB0R ERROR FOR ELE MENT , I 4 , 3X , 6H AND/ OR, I 4 ) 

105 FORMAT ( 1H 0 , 45X , 3UHOONDUCTION ELEMENT TABLE ( T D 1) / / ( 5X , I 4, 2 X , 7{ 2X , I 
22, 1H, , 1PE L1.4) / 10X,6( 2X,I 2 , 1H, , 1PE11.4)// ) ) 

lOo FORMA T ( 16F5.U) 

107 FORMA T ( 1H 1 » 50X , 29HELE MEN T STARTING TEMPERATJRES// ) 

108 FORMAT! 10 ( 2X, I 3, IX, F6. 1) ) 

109 FORMAT! 25X, 74H80UNUARY COEFFICIENT MATRIX -( H* A ) BTU/F . S EC AND ( H*A 
2*T-B) JR rt*C-P 3TU/SEC/ /1X»1LHBNDRY COND . , 2X , 1 3 HF I X ED FLJX.6X, 17HF 
30RCEU CON VEST ION, 10X, 1 SHF REE CONVtC T I ON , L 4 X , 9 HR ADI AT I ON, 12X, 19HVAR 
4. TEMP .CON VECTION/iX, l IHELt M. ELEM. , 6 X, 1 H W , L 1 X , 3 HH* A , 3 X , 7HH*A*T -8, 3 
5X , 3rth*A,SX,7Hri*A*T-3, 6X,3HH*A,8X,7HH*A*T-E,3X,3HH*A,8X,5HN *C-1V/( 1 
6X,2( I4,1X),1X,9(1X,E12.5) ) ) 

110 FORMAT (1H0, 3 OX, 58HVAR.TEMP. CONVECTION DIMENSIONLESS COEFFICIENTS ( 
2H*4/W*C-P )//( 6( IX, 214, 1H, ,1PEL I .4) ) I 

111 FORMAT ( LM0,<t5X,3t>HFLUlL) SHEAR HEAT GENERATION, BT J / S EC// ( 6 ( IX, 2 1 4, 
21 H» , 1PE11 .4) ) ) 

112 FORMAT { 1H0, 2X , 1 OHDA T4 TIME =,F 7. 4 ,2 X,7HMINUTES ) 

113 FORMaTI 1H1,2UX,6HTEMPERATURES (F) AND HEAT FLJXES (BTJ/SEC) FOR E 
2ACH ITERATION) 

114 FORMAT! 3( 4X , I 3 , 1H , , 1 PE 14. 7 ,2 X , 1 PEI 4 . 7 ) ) 

115 FORMA T ( 1H l, 30X, 32H C ONVERG b NCE CRITERIA, I T ER AT I 0N= , 1 4/ / 23< , 8HC0N D. 
2NJ . ,6X,9HNEw TEMP. ,L0X,9HGLD TEMP. , 9 X , 1DHDI F F ERENCE , 8X , 9HM AX .OIFF. 
3//) 

116 FORMA T ( 1H0, 2X , 42HC ONVERGE NCE CRITERION NOT MET, I T ERAT IONS = , 16 ) 

117 FORMAT! HO, 2X , 43H C UNVi RGE NCE CRITERION ACHEIVED, I T ERAT ION S = , I 6 ) 

118 F0RMAT(1H1,20X,56HELEMENT( AT CENTER) TEMPERATURES (DEGREES F), ITE 
2RAT 1 3 N S= , lo//) 

119 FORMAT (6( 2X, I 4, 1H, , 1PE 14. 7) ) 

120 FORMAT! IH L,9X, 102HELEMENT FORCED CONVECTION, FREE CONVECT., RADI AT 

21 UN AND VAR. TEMP .CONVtCT. C OND UC TA NC E S ( H*A) BTU/HR.F//) 

121 FORMAT! 16X, 21 6,4519.7) 

122 t-ORMAT ( 1H i,45X, 38HVA.R . TEMP. C0NV6CT1 ON BOUNDARY T EM PS . ( F ) // ( 6( 1 X , 2 1 
24, 1H, , 1PE L 1.4) ) ) 

123 FORMAT ( lri 1, 42X, 43HELEhcNT THERMAL CONDUCTIVITIES BT J/ HR. FT . F/ / J 
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124 FORMAT T 3X , 3 5H ELEMENT ENrfS CONDUCTION HEAT F L J X ES , 2 X , I 3 , 4( 2 X , EL 4 . 7 ) 
2, 2X,6HBTU /HR ) 

125 FORMAT! l*X,45Ht Le ME NT ENtfS CONDUCTANCES TO NA BURS ( BT J / HR . F )//< L OX, 
214, 44 2X,E 14.7) ) ) 

126 FORMAT! 1HO, 10X, 104-1 ELE MEN T BOUNDARY NO, CONDUCTION NO, BND RY CONDJ 
2C T ANC E { H* A ) SUM ( 3 TU/H R • F ) AND BNURY HEAT FLUX SUM! BTU/HR )/ / < 3( 3X f 2 
314, 2X , h 14 . 7, 2X,E14.7) ) ) 

127 FORMAT! 1H 1 f 40X » 43HELE ME NT FREE THERMAL EXPANS I UN» ALPHA*( T-TO )/ / l 

128 FORMA T ( 6E 13.6) 


CALL TIME 1 ( START) 

READ RUN PARAMETERS 
60 DO 1 J=i, 10 
L D { J ) = 0 

1 CD( J ) ■= 0 • 

00 2 J =1, 5 
N H J ) = 0 

2 RD( J ) =0. 

DO 3 J=l, 16 
DO 3 K=l, 11 
MATH J ,K ) =0. 

3 ALPHA! J , X ) =0. 

NX = 0 

CALL CDAl LD, MI ,KJ ,C3 , MATl » ALPHA) 
IF (NX.NE.O) SO TO 99 

N D IM = 4 00 
NBD IM =200 
NVDIM=S5 
iMJPTS=C'J( 1) 

N3P T S =C J ( 2) 

\l I T S = C D( 3) 

Ni< = 0 

ZERO COMMON STURAbt 
DO 5 J PT = 1 , NJPTS 
Ti JPI ) =0 . 0 
0 ( JPT ) =0. 0 
J BN T( JPT) =0 
JNTF! JPT) =0 
THEXP { JPT) =0. 0 
DO 4 J = l f 4 

4 CE( JP T, J ) =0.0 
DJ 3 J = L , 19 

3 TDK JP T, J )=0. 0 

DO 6 NBPT =1,N3PTS 
N 3 N D ( NBPT ) =0 
SUM H( NBPT ) =0. 0 
SUM T( NBPT ) =0. 0 
DO 6 J =1, 9 
6 H4N6P T , J ) =U.O 

DO 7 K=i, NVJI M 
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GFf K )=0.0 
TCF(K)=0.0 
CFN(K)=0.0 
AFN(K) = 0.0 
NBR(K)=0 
N ODF ( K ) =0 
T END ( K >=0.0 
C 

on B K = 1 . 5 
NM(K ) = 0 
M AT ( K ) = 0 

8 N PT ( K ) =0 
C 

C READ ELEMENT GFOME TRY DATA ♦ CHECK CARD ORDER 

on 10 JPT=1,NJPTS 
R F AD ( 5, 10 2 ) (T n l(JPT,J),J=l *13) 

C 

N PC = TD H JPT* 1) 

IF ( JPT-NPC ) o, 10,9 

9 WRITF( 6* 1 03 ) JPT 
NK=1 

10 CONTINUF 

IF (NK.NF.O) GO TO 99 
C 

C C h FCK FOR CONDUCTANCE NABOR ERRORS 

NK = 0 

DO 1 5 JPT=1,NJPTS 
DP 15 J =5 , 8 
N =T D H JPT. J > 

IF (Ml 14.15.11 

11 IF ( JPT-N > 12. 14.12 

1? RJPT=JPT 

J N=J~4 

GO TO ( 13.212,3 13.413) .JN 
18 IF (RJPT-TD1 (N, 7>> 14,15,14 
218 IF (RJPT-TD1 <N» 81) 14,15,14 

818 IF (PJPT-TD1 <N,5) ) 14,15,14 

418 IF (RJPT-TD1CN.6)) 14,15,14 

14 WRITE! 6, 104 ) JPT.N 
NK= 1 

15 CONTINUF 

TF (NK.NF.O) GO TO 99 
C 

TF (LD( 1) .EO.O) GO TO 16 
W P I T F( 6, ion 

WRIT F( 6, 105) { JPT,(J,TD1(J°T,J ) , J= 1 ,18) , JPT= 1 , N J PT S ) 

C 

C C PM PUT F ELEMENT BOUNDARY AREAS, CONVERT ID 1 DATA TO FEET 

16 DO 17 JPT = 1 ♦ NJPT S 
RO=TDU JPT, 2 )+TDl( JPT, 3) /2. 

P T = TDl( JPT, 2 )~TD1( JPT,3)/2. 

TTK JPT, 1 4 > = .048 6332 3* TD1 ( JPT, 2)*TD1 (JPT, 3) 

TCU JPT, 15 >=. 0436332 3*TD1 (JPT, 4) *R0 
TC1( JPT, 1 6 > = . 04363323^ TDl(JPT,4) 5{c RI 
T ri( JPT, 2 ) = TD1 (JPT, 2 1/12, 

TC1( JPT,3I=TC1 ( JPT,3)/12. 
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r> r> n n r> on on on no 


17 Tf:i{JPT,AhTai{JPT,4)/12. 

READ INITIAL TEMPERATURE GUESSES, CURVE FIT CONDUCTIVITIES 
P FAD! 5*106) ! T ( JPT ) * JPT-1 *N JPT S) 

WRITE! 6, 107) 

WRITE! 6* 1 0 R ) { JPT* T(JPT) * J P T = 1 *NJPTS) 

00 IB JPT=1*NJPTS 
TCI! JPT* 19 ) = T( JPT) 

IS CALI FIT! TD1,JPT,NDIM*MATL ) 

READ BOUNDARY CONDITION DATA AND COMPUTE COEFFICIENTS 
IF (LD!5).NE*0) WR I TF ( 6 * 1 01 ) 

NK=0 
N TF=0 
N CR = 0 
N TR = 0 

DO 19* N BP T= 1 *NBPTS 

19 C At L CON! NBPT*NBOIM*HtNBND*TBND*TDI *TDA . TD2 * NDI M * T * LD, NK* NV DI M * NOD 
IF.JNTF) 

I F (NK.NE.O) GO TO 99 

CCN PRINTOUT 

IF ( LD< 2) •FO.O) GO TO 20 

WRITF! 6*101) 

WRITE! 6,109) C NB PT * NBND ! NB P T) , ! H ( NBP T * J ) *J=1,9) * NB PT - 1 * N BPT S ) 

20 I F ! NT F • EG .0 ) GO TO 21 

WPTTF! 6,101 ) 

WRITE! 6, 1 10) l J, NODE I J) *TCF ( J) *J = 1 * NTF ) 

WRITE! 6* 11 1 ) ! J , NODE ( J ) * QF ( J ) ,J=1*NTF) 

PERFORM ITERATIVE ANALYSIS 

21 CALL TIME l (TFX) 

T CAT = ( TEX- START )/360C, 

WRITE! 6,1 12 ) TDAT 
NCC = 0 
NPRT=RD( 3) 

DP 29 NIT-UMTS 
NK = 0 
N CP = 0 
N T F= 0 
NCR=1 
N TR = 0 
Nf C=NCC+l 
L D( A ) = 0 

DP 22 NBPT=1 ,NB P TS 

22 K AN ( N B P T ) = 0 

BOUNDARY CONDITION ITERATION 
DO 23 NBP T= 1 ,NRPTS 
j PPT^NBND! NBPT ) 

C AH OT!NBPT,NBDIM,NDIM, JB PT * H , SUM T , SUMH , TD1 ,T BNO , T * CD, NVO I M* N IT * N 
lOCF, JNTF) 

23 0(JBPT)=SUMT (NBPT) 

0 VFR-RELAXAT ION 
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MCR=0 

no 24 jpt=i,njpts 

24 TAIL TAT.S!TO]*JPT*NDIM*NBND,NBDIM*Q*T,H*RD*LD*€£*NK) 

IF CNK.NP.O) 0 0 TO 99 

C 

IF (LD< 10) , F 0 • 0 ) GO TO 25 
WRITF! 6*113) 

WRITE! 6*114) ! JPT* T( JPT) *0( JPT) *JPT=l ,NJPT$) 

C 

C CHECK FOR CONVERGENCE * CORRECTION OF TOi TEMPERATURES 

25 E RP =CD ( 4 ) 

IF (NIT.FO.N TTS ) GO TO 30 
IF INCC.NF.NPRT ) GO TO 26 
N CC = 0 
LCI 4)=1 

WRITE! 6, 115) NIT 

26 DO 27 JPT=1*NJPTS 

27 CALL CONC HK < JP T * ND I M * TD1 *ERRtLD*NK*T) 

IF (NK) 28*31*28 

RECAIC. K FROM NEWEST TFMP* 

28 DO 29 JPT=UNJPTS 
TT1I JPT* 19) = T( JPT) 

29 CALI FIT! TD 1 * JP T*NDI M* MATL > 

30 WRITE! 6* 1 16) NIT 
GO TO 33 

31 WRITE! 6* 1 17) NIT 

TFMPERATURE PRINTOUT 
IF !LD! 3) ) 32*33*32 

32 N ITM I =N IT— 1 
WRITE! 6*118) NTTMl 

WRTTF! 6*1 19) ! J PT* TO 1 1 JPT * 1 9 ) * JPT= 1 * NJ PT S ) 

33 WRITE! 6*118) NIT 
WRITE! 6*119) { JPT* T! JPT) , JPT=l *NJPTS) 

IF (LO! 8) .NF.O) PUNCH 106, (T(JPT) * JPT=1 *N JPTS ) 

CCNVFRT FROM SFC.TO HR S* 

DO 34 NBPT=1*NBPTS 
$UMT!NBPT)=SLMT< NBPT >*3600. 

SUMH!NBPT > =SUMH! NBPT )*3600. 

DO 34 J =2 * 8 * 2 

34 H!NBPT*J )=H(N8PT*J)*3600. 

PRINT OUT 80LNDARY CONDUCTANCES(H*A) 

WRITF! 6* 120) 

WRTTF! 6* 121) ( NBPT * NBNO ( NBPT) , <H(N8PT*J) * J= 2*8 *2 ) * NBPT= 1 , NBPTS ) 
PRTNT OUT VAR* TEMP * BOUNDARY TEMPS. 

IF (NTF.NF.O) WRITE !6* 122) ( J , NODE < J ) * TBND ! J ) , J= 1 * NTF) 

PRINT OUT ELEMENT THERMAL CONDUCTIVITIES 
DO 35 J PT= 1 * NJ PT S 

35 T ri( JPT, 18) = TDK JPT* 18) *3600. 

WPITF! 6,123) 
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WRITFC 6, 1 19 ) (JPT,T01(JPT f 18), JPT=1 ,NJPTS) 

COMPUTE AND ViP I TE HEAT FLUXES 
DO 36 JPT=1,NJPTS 
OC 36 K=l,4 

36 CPI JPT,K)=CE< JPT, K >*3600. 

WRITF(6,]01) 

OP Aft JPT=1,NJPTS 
NE=TDl( JPT, 5 ) 

IF CNF) 38, 37, 3ft 

37 0 JF = 0 • 0 
GO TO 30 

38 0 JE=CE ( JP T , 1 )* t TCNE)-T(JPT) ) 

39 N N=TDl ( JP T , 6 ) 

IF (NN) 41, AC, 41 

40 0 JN = 0. 0 
GO TO 42 

41 OJN=CE< JPT,2 )*< T(NN)-T< JPT) ) 

4? N W= TOl ( J D T , 7 ) 

IF ( NW ) 44.42,44 

43 QJW=0.0 
GO TO 45 

44 OJW=CF( JPT, 3 )*( TCNWJ-TC JPT) ) 

45 NS=Tni( JPT, 8 ) 

IF (NS) 47,46,47 

46 OJS=0.0 
GO TO 4ft 

47 0 JS=CE( JPT, 4 )*< T(NS)-T! JPT) ) 

4 8 WRTTEC 6, 124) J P T ,0 JE , 0 JN , OJ W,Q J S 

WRITEC 6,101 ) 

WRTTEt 6, 125) ( J PT, ( CE ( JP T , J ) , J = 1 ,4) , JPT = 1 , N JPT S ) 

PRINT OUT BOUNDARY CONDUCTANCE SUMS AND HEAT FLUX SUMS 
WRITFC 6,101) 

WRITF( 6, 126 ) (NBPT,NBND(NBPT) • SUMH (NBPT) , SUMT ( NBPT ) , NBPT =1, NBPTS ) 

TFFRMAL EXPANSION CALCN. AND PRINTOUT 
IF (LD< 10) .NE.O) WR I TF ( 6 , 1 01 ) 

DO 49 JPT=1,NJPTS 

C ALL F IT?( TD ltNDIM, JPT, ALPHA,THEXP) 

49 CALL DEFORM! JPT,NDIM, TD 1,RD*LD,THEXP) 

WRITE! 6, 127 ) 

WRITEC 6, 1 19 ) (JPT, THEXP ( JPT) , JPT=1 ,NJPTS) 

IF ( LD ( 9 ) • N E #0 ) PUNCH 128, (THEXP(JPT) ,JPT-1,NJPTS) 

CHFCK FOR ADDTL .CASES 
IF (NIC 11-1 ) 99,99,50 

50 GO TO 60 

99 CALL EXIT 
STOP 
END 
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$ IBFTC CL) AM 

SUBROUTINE CDA!LJ,NI , RO ,C0 , MAIL , ALPHA , NK) 

PROBLEM PARAMETER INPUT, CONDUCTIVITY AND EXPANSION INPUT 

DIMENSION LD! 10),NI (5) ,RU!5) ,CD(10> , MATL! 15 , II ) , ALPHA! 15,11) 
DIMENSION! FILEIll) ,WCOM{12) 

REAL MATL 

10 1 FORMAT ( 12A6) 

102 FORMAT ( 1H i. * 2 9X , 12A6//) 

103 FORMAT! 1016) 

104 FORMAT 1 1 H 0 , 1 3 X , 3HLD = , 1 0 ( 3 X , I 2 , 1 rl , I 3) ) 

10 5 FORMAT! 1H 0 , 1 3X , 3riN I =, 5 ( 3 X , I 2 , 1 H , I 3 ) ) 

106 FORMAT! 10F6.0) 

107 FDRMAT!lHO,L3X,3r!RO=,5I3X,I2,lH,lPE13.6)) 

10 8 FORMAT! 1HO,13X,3HCO=,5(3X,I2,H,1PE13.&)/16X,5!3X,I2,1H,1PE13.6)) 

109 FORMAT! 1H0, 2X , 58HND. CONDUCTION OR BOUNDARY ELEMENTS E< CEFDS ALLOTS 
2D STORAGE ) 

110 FORMAT ! 5F 10.0) 

111 F0RMAT(6E 13.6) 

112 FORMAT ( 1H0,59X, ILHDATA FOR K. / ( 17X ,6 < 2X , I 2, l PE 14. 7 ) ) ) 

113 FORMAT (1H0,59X, 15HDATA FOR ALPH A./ < 1 7X ,6 ! 2 X , I 2 , l PE 1 4. 7 ) ) ) 

READ! 5,10 1) WCOM 
WRITE 16, 102) WCOM 

GENERAL DATA 

R EAD( 6, 10 3) (LD! J) ,J=i ,10) 

WRITE! o,l 04) { J,LD!J> ,J=1 ,10) 

READ! 5, L03) <NUI ) ,1=1 ,5) 

WRITE! 6, 10 5) (I, N 1(1) ,1=1,5) 

READ! 5,106) (RD(I ) ,1=1,5) 

WRITE! 6, 107) ( I ,RD< I ) ,1 = 1 ,5) 

READ! 5,106) !C0< J ) ,J=1 , 10) 

WRIT£!fa,108) ( J,CD(J) ,J = l,10) 

IF (400. -CD! 1) ) 2,1,1 

1 IF < 200. -CD! 2) ) 2,3,3 

2 WRITE! 6, 109) 

NK = 1 

GO TO 11 

CONDUC TIV ITY DATA 

3 DO 6 L =1, 15 

READ! 5,110) (FIlE(NCO) »NCU=1*5) 

IF IF ILEl 1) ) 4, 7,4 

4 MN = FILE(2 ) 

NP 6=F ILEt 3)^o. 

READ! 5,111) (FILEINCJ) ,NCQ=6,NP6) 

00 5 N C J = 1 , NP 6 

5 MATL(MN»NCU)=FIL£(iNCQ) 
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6 WRITF(6,112) (NCOt FILE(NCO) ,NCQ=l,NPM 
EXPANSION DATA 

7 nr 10 i =un 

RF AO (5, 110) (FTLEINCQ) ,NC0=1 *5) 

T F ( FILFf 1 ) ) 8,11,8 

8 M N= F TL F ( 2 ) 

N P 6 = F I L F ( 3)+6. 

P FAPI 5* 11 1) (F TLE ( NC G 1 , NC Q=6*N P6) 

no q mco= l,MPft 

9 ALPHAfMN, NCO ) = F ILE ( NCQ I 

10 W R I TF( 6*1 17) f NCQ , F I LE ( NC Q) , NC Q= 1 , NP6 1 

11 RFTtJPN 
FND 


: FTTM 

S UBRniJT IN E F IT ( Till* JPT , ND I M , MATL ) 

FIT CONDUC T I V I TY (K) FRO* PPL YNOM I N A L CURVE 

0 TMFNS TON TCHNOIM* 1 9) ,C OS ( *5) , CA l C (6) ,MATL<15, 11) 

RFAl M A T L 

C 

101 FORMAT! 1H 0 , ? CX , 4 1HTE MPFR A TURF EXCEEDS CURVE FTT LIMIT FOR K) 
C 

TFMP= TDKJPT, 191 
MAM = TD1(JP1,9) 

DO 1 Mf,0=l,5 

1 C GS ( NCO ) = MATE (MAN* NCO) 

1 F< TFMP-CGSI 1)1 3,3,2 
? WPTTF! 6, 1 01 1 

CO TO 9 
3 NOP = COS( ? 1 
NOP 1 =NOP + 1 
DC 4 NCQ=l,NCP 1 
Lf = 5+NCQ 

<+ CA! C (NCO ) - M ATL ( MAN , LC ) 

1 F ( f OS ( 2 ) ) 6 , 5 , 6 

5 Y = CAT C( 1 )*r.QS( 4 1 
on TO B 

6 Y = CALCtNOP 1 1 

T FM l = TFMP*C0S<5) 

no 7 ’mp y=i,ncp 

J 2 = NOPl-NP Y 

7 y = CAL C ( J 2 ) +Y*TEM1 
y = Y* CO S ( A ) 

8 TC1IJPT, 1 8 ) = Y 
C 

9 R FTURN 

fno 
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SUbROU TIN E C«!N(N8PT,NBDlM f H , NBNO ,T8ND,TDL ,TDA,TD2, NDI M , T , L 0, NK , N V D 
IIM.NODE, JNTF ) 

COMPUTE HE A T TRANSFER CUfcFF ICI tNTS 

DIM EN S IUN HI N3D IM,9) ,NdND ( NbDI M) , T3 N0{ NVD I M) , T D A I 7 ) ,T 02 ( 23 ) 

DHEMS ION TDK NUiM, 19) , TINDI M) , LD( 10) , NODE ( Ns/D I Ml , J NT F ( NDI M ) 

CJMM0iN/BLKl/J3NT(400) ,RA,N(200) ,NTF , NOR , NT R, QF ( 9 5 ) , NQR 

COMMON /HLK2/TCF ( 95) ,GFN! 95) ,AFN(95) ,N3R(95 ) , NM(5), NAT ( 5 ) f NBH 5 ) 

101 FORMAT! 31 o) 

102 FORMAT ( IHO, 15X, 6JH3UJNDARY I INDICATOR CARD URDE R ERROR, BOJNDARY 
2EL EMEN T NO. = , 14, 11H CARD kEADS,14) 

103 FORMAT ( 13, IX, F6. 0) 

104 FORMAT ( HO, 10X, 7lHaOUNDARY INDEX (TYPE l OR 7) CARD ORDER ERROR, 
2CJNDUCTIJN ELEMENT NJ.=,I4,11H CARD READS, 14) 

105 FORMAT! Id, Eli. o,I 6) 

106 FJR:MAT(6F6.0,F12.0) 

107 FORMAT! 7F 11. 0) 

1C 8 FORMAT ( 1H0, 10X, 71HBUUNDARY INDEX {TYPE 2 TO 5) CARD ORDER ERROR, 
2C0NDOC riON ELEMENT N3.=,I4,IH CARO READS, 14) 

109 FORMAT ( 1H0, 2X , 8HELEME N T= , F 4. 0 , 2 X ,5 H T YPE= , F 2 . 0 , 2X , 3 HS J BT YPE = , F2 . 0, 2 
2X»5HNEXT=,F2.0,2X,llHITbMS RE AD = ,F i . 0 , 3 X , 1 OH ARE A CODE=, F2.0, 2X, 1 3H 
3SPEC IAL AREA = ,E 12. 5,1 X,7H SU.FEE T) 

110 FORMAT! 1H0,2X, 17HD0CT FLOW, LI U Ul D /2 X ,4HT-B= , F 5 . 0 , 2X , 4 HK-F = , El l .4, 
22X , 4HC-P=,El 1.4,2 X,3HMU=, El 1.4 ,2 X,10HMASS FLOW = , El 1 . 4 , 2X , 7 HD-HY 0R = 
3, E 1 1 . 4 , 2X , 7HL-FLJW = ,E l 1.4/2 X,3HRE= ,E1 1 .4 , 2 X , 3HPR= , E 1 1 . 4, 9< , 8HH COE 
4FF = »C 14. 7 , 2X » 1 IHC OE FF * ARE A= , E 1 4 . 7) 

111 FORMAT! lr!0,2X,15HSIDES OF ROTORS/2 X ,4HT-B= , F5 . 0 , 2X , 4HK -F=, E 1 1 . 4, 2X 
2, 4HC-P =, E 1 1.4, 2X, 3HMU = ,E 1L . 4 , 2 X , 3H NU= , El 1 ..4 , 2 X , 5HR- AY = , Ell .4, 2X, 4H 
3GAP=, E 11. 4, 2X, 6H0MESA = ,F6. 0/4HR-0= ,E1 l . 4 , 2 X, 3 HFL AG= , F2 . 0, 9 X , 3HRE=, 
4E11.4, 2X, bHPR** 1/3 = ,;1 l.4,2X,3HH CQEFF=,E14. 7 , 2X , 1 1 HCOEFF* AR EA=, El 
54.7 ) 

112 FORMAT ! 1H0,2X, 1 5HRADI AL SEAL G A P/2 X ,4H T— B= , F5 . 0 , 2X , 4HK-F= , E 1 1 . 4, 2X 
2, aHC-P = ,c 1 1.4, 2X, •JHMU = ,E1 1. 4,2 X,3HNU= , El l . 4 , 2 X , 5HR- AY = , E 1 1 . 4, 2 K , 4H 
3GAP =, E 11. 4, 2X, 6H0MEGA=, F 6. 0/2 X,4HK-0=, El 1 . 4 , 1 5 X , 3HRE= , E 1 1 . 4, 2X , 8H? 
4R ** l/3 = ,E 11.4, 2X, 5 HRe-C = ,E1 1.4/2 X.8HH C UEF F= , E 14 . 7 , 2X , 1 1HCUEFF*ARE 
5A = , E 1 4 . 7) 

113 FORMAT! 1H0, 2X , 20HC0NC ENTKi C C YL I NDERS/2X , 4 HT -B= , F5 . 0 , 2X , 4HK-F= , E 1 1 
2.4, 2X, 4HC-P = ,E1 1.4,2X, 3HMU=,E1 1 . 4,2 X,3HNJ= , E 1 1 . 4 , 2 X , 4HR-0= , Ell .4, 2 
3X.4HR- I=, Eli. 4,2X ,6H0MeGA =, F6. U/2 X , 5H Y- AX= ,E1 1 . 4 , l 5X , 3HRE= , Ell. 4, 2 
4X, 8HPR** 1 / 3= , E 11.4»2X,<JHH C0EFF = ,E14.7 ,9 X , 1 1HC0EFF*ARE A= , E 14.7 ) 

114 FORMAT ! 1H0,2X, 10HFLAT PLATE /2 X,4HT-B= ,F5. 0 , 2 X , 4HK-F= , E l 1 .4 , 2X , 4HC- 
2P = , El 1 .4, 2X, 3HMU = » E l 1. 4 ,2 X,3HNU= ,tl l. 4 ,2X , 7 HL-FLOW = , E 1 1.4, 2X, 7HY-F 
3L0rt=, E 11. 4/2X, 3HRE = ,E1 1 . 4 ,2 X , 3HP R= , E 1 1 . 4 , 9 X, 3 HH COEFF= , E14.7,2X, 11 
4HC0EFF*AR EA = ,E 14. 7) 

FUTURE FORCED CON VEC . CCE= F . FORMATS (115-8) 

119 FORMAT! IHO, 2X, 39HHURI ZONTAL CYLINDER OR VERTICAL SURFACE/ 2X, 4HT-B= 
2, F5.0, 2X, 4HR-F = ,E 11. 4 ,2X,4HC-P=,E1 1 . 4 , 2X , 3 HM J = , El 1 . 4 , 2 X , 3HNU = , Ell. 
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34, 2 X, l 7HL-FLUW JR J-HYUK = ,fcl L. 4 / 2 X , 3 iiPR= » C l 1.4 ,2X, 3 HG R= » El 1.4, 9X, B 
4HH CUEFF=,E14. 7,2X,1HCUEFF*AREA=,E14. 7) 

120 FORM A T ( 1H 0, 2X , 19HLIQUI0 = l LM C 00L1 NG /2 X ,4HT-b= , F5. 0 , 2X , 4HW / L = , E 1 l . 
24, 2X, 4HC-P = ,E 1 1.4,2X, 7 ML -F L 0W= , E 1 1 . 4 , !0X,8rH C0EFF= , El 4. 7, 2X, l IHC3 
3EFF#ARcA=,E 14. 7) 

FuTUkE FR EE CGNVEC .CUEFF. FURM ATS (121-3) 

12 6 FORMA T ( lHO,2X,9HR4JI4TIJN/2X,4HT-t=,F5.J,2X,l3HFACrOR-SE=, El 1.4, 2< 
2, 10HFACTJR-ES,E 1 1 .4, 2 X, 8HEM1 3. -E = , F 1 1 . 4 ,2 X ,3 ME MI S . -S= , El 1 . 4/ 2< , 8HH 
3 CUEFF =»E 14. 7, 2X, 1 IHC0EFF*AREA= ,E14. 7) 

127 FjRMA T ( 1H0,2X »21HSPECI FI ED COEFF IC I ENT/2X ,4HT-B= , F 5 .3 , 2X , 3HH COEFF 
2=,E 14. 7, 5X, 1 1HC0 Ei-F*ArEA=,E14. 7) 

128 FORMA T( HO, 2X,2o-IVAk. TEMP. DUCT FLOW, L I 00 1 0 / 2 X , 4HT -6= , F 5 . 0 , 2X , 4HK- 
2F = , El 1 .4, 2X,4HC-P=,E 1 1 . 4 , 2 X , 3HiMU = , E 1 1 . 4 ,2 X , 1 3HMAS5 FL0W=,E 1 1 .4, 2X, 
37HO-HYOR=»El 1.4,2X,7BL-FLOW=,fcl 1 . 4/ 2 X , 9HANN. NBR .= , F4. 3 , 2X , 9HFL OW N 
4BR = , F 4 . 0, 2X , 7HF L AG-C = , F 2. 0 ,2 X , 7HFL A3— I = ,F2.3,2 X,6HX— S EC= , E 1 l .4/2X, 
53HRE=, Ell .4,2 X, 3HPR=, El 1.4,6 X,8HH C OE FF= , E 14 . 7 , 2X , 1 1 HCUEFF *ARE A= , E 
614.7) 

129 FDkMATI IH0,2X,24HVAK. TeMP. RADIAL SEAL GAP/2X , 4HT-B= , F5 .0, 2X , 4H<-F= 
2, El 1.4,2X,4HC-P=,Eli. 4 , 2 X , 3HMU= , cl 1 . 4 ,2 X, 3 HNJ = , El 1 .4, 2X , 5HR-AV =, E l 
3 1 • 4* 2X , 4H GAP = , E 1 1 . 4 , 2 X , 6HUMEG A = , F6 . 0 / 2 X , 4 H k— G= , El 1 . 4 » 2 X » 7HFL AG— C=» 
4F2.0, 2X, 7HFLAG- I = , F 2. 0 , 2 X , 9HANN . NBn. = , F4. 0 , 2 X , 9 HFL G.NBR .= , F4.D, 2X, 
5 10HMA S S FLL);V = ,Lll.4/2X,3HRE = ,E11.4,2X,dHPR**l/3=,E11.4,2X, 5HRE-C = , 
6E11 .4, *X, BHM Cuc-F=,E 14. 7 ,2 X»1 l HCOLFF* ARE A= , E 1 4 . 7/ 2X, 14HHr .GEM. FAC 
7T3R=, E 12. 6) 

FUT OR E VAR. TEMP. CUEFF. FORMATS (133-1) 

132 FORMAT ( 1H0.2X , 33HCUM3 INEJ FLOW ELEMENT, CGNUJCT ION N0.= , 14 , 2X, 24HF 
2L0W FROM CJNU.tLEMtNTS , I3,4 HAnD ,14) 

FU1URE COEFF . SUBPrOGk. FORMATS (118-3 ( l-ORCEO) ,121 -5 (FREE), 130-1 
( VAR.TcMP . ) 

KEAO INDICATOR CARD, CHECK ORDER 
READ( 5,101) i'jBP, MtJNJl NBP ) , INDtJ 
IF ( NBPT-NoP ) 1,2,1 

1 rtR I TE ( -j, 102) NtJP T,.NBP 
NIK = l 

SET PARAMETERS 

2 J0PT=N3ND ( N3P ) 

J BN T ( JbPT )=N0PT 
TDK J3PT, 17) = l NOB 

SPECIFIED TEMPERATURE 
IF ( INOB.NE. 7) GJ TO 4 
REAO( 5,103) JBP , T ( J 0 P ) 

IF IJBP.NE.JBPT) GO TO 3 
TD 1 ( J 3 P T , 19) =T< JBP) 

GO TO 99 

3 WR ITE( 6, 104) JBPT, JBP 
N<= 1 

GO TO 59 

4 GO TO ( 5, 7, 7, 7, 7, 7) ,1 NOB 
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C SPLCIFIElI heat flux 

5. READ( 5,10 5) Jr) P » H ( NbP T , 1 ) , I AUUQ 
IF (JdP.tQ.JbPT) GU TG 6 
WRITE! 6,104) JdPT,j3P 
NK= 1 

6 IF ( IADDg.Ej.O) GO T3 99 

BOUNDARY AREAS 

7 AE= TU 1( J3PT, 14) 

AN= TD 1 1 J3 PT, 15) 

AS = TL) 1( JdPT, 16) 

READ DATA F JK TYPES 2 TO 6 
90 READ! 5,106) ( TD A ( J ) , J =1 , 7 ) 

N I T M= T OA ( 5) 

R fcA D( 5,107) ( TD 2 ( J ) , J = 1 , HI TH) 

oHECK. CGND JC T ION NG . M A TC H 
J BP =T Da( 1 ) 

IF (JbP-JbPT) G,9,o 
b WRITE! o, 103) JBPT ,JBP 
NK= 1 

9 IF (NR .Nt .0) oG TG 99 
ITYPL- = TDA ( 2) 

GO TG (99, 10,29,46,43,55) , I TYPE 

FORCED C3NVEC 1 I UN 

10 ISTYPE=TJA( 3) 

GG TG (11 , 13, 15,17,19) , I STYPE 

DULT FLGW, LIQUID 

11 CALL j,Il( TDA, TD2, At, AN, AS) 

H( M BP T,2) =H( NbP T , 2) «- TD2( 16) 

H{ NtiP T ,3) =H(N3PT,3)*-TJ2(16)*TD2(1) 

IF (LD(5).EJ.0) GG T3 12 
WR I TL ( fc, 1 09) (TJA(J) ,J=1,7) 

WR I Tel 6, 1 10) ( TJ2 ( J ) , J=1 ,7) , ( TD2 (K) ,K=13 , 16) 

12 IF ( T JA( 4 ) ) 90, 99, 90 

SIDES GF RGTGnS 

13 CALL SN2( TDA, TD2, Ac, AN, AS) 

H( NbP T ,2) =H( NdPT,2)<-TD2( 14) 

H( NbP T ,3) =H(Nl)PT, 3)+TD2< 14) * TJ2 (1) 

IF (LD(5).tJ.O) SU TG 14 

WR I TE ( 6, l D9) ( TD A ( J ) , J =1 , 7 ) 

wR I Tt ( 6,1 11) ( TJ2 ( J ) , J=1 , 14) 

14 IF { T DAI 4 ) ) 90,99,90 

RADIAL SEAL oAP 

15 CALL 3N3( TDA, TD 2, Ac, AN, AS) 

H( NbP T,2)=H(.NHPr,2H-TJ2(I9) 
H(NdPT,3)=H(N3PT,3)*TJ2( 19) *T02 (1) 

IF ( L D ( 5 > . El . 0 ) SU TD Lb 
WR I Tt ( 6, 109) ( TD A ( J ) , J =1 , 7) 

WR i Tt ( 6, 1 12) ( TJ2 l J ) , J = 1 ,9) , C TD2 (K) ,K=15, 19) 
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16 IF ( T DA l 4 ) ) ‘10,94, 90 


CONCENTRIC CYLINDERS 

17 CALL S N4( TJA,TD2,AL,AN»AS) 

H(NBPT,2) =H(NJPT,2)+TL)2<13) 

FUNBPT,3)=H(N3PT,3H-rj2{i3)*Tu2(l> 

IF (LD( 5) .tO.O) GO n 18 

WRITE! 6,109) ( TDA ( J ) , J = l , 7 ) 

WRITE ( 6, 1 13) (T02(J),J=l,I3) 

18 IF ( T0A( 4) ) 90, 99, 10 

FLAT PLATE 

19 CALL SN5! TDA ,Tu2,AE ,AN ,AS) 

HIM BP T,2) =H!NBPT, 2)+TJ2( 1 L) 

F!<NBP f ,3) =F!< NBPT, 3 H-T0 2! ilKTilZU) 

IF <LJ!5).EG.O) GO TO 20 

WR I TE ! 6, l 09) ( TJA I J ) , J=1 , 7) 

WR I TE ! j, 1 1 4) ( T02IJ) , J=1 , 1 1 ) 

20 IF ( T DAI 4 ) ) 90,99,90 

FOR FUTURE SUBROUTINES 14 MORE) 1ST AT E M. NOS 2 1-2, 2 3-4 , 26-6 , 27-8 ) 

FRbE CONVECTION 

29 ISTYP E=T0A! 3) 

GU TO ( 30 , 32, 34) , I STYPt 

HU* izontal cyl injek 

30 CALL FCU TJA , TD2,AE, AN, A3, T,J8P T.NDI M) 

HINBP r ,4) =H(NBP r, 4)4- TU2I LO) 

H(NBPT,5)=H(Ni3PT,5)*-Tj2(10)*T02(1) 

IF !LJ!5).EU.0> Gu TO 31 
WRITE! 6,109) ( TUA ( J ) , J =1 ,7) 

WR I TE ( 6, 1 1 9) (TJ2IJ) , J = 1 » 10) 

3 1 IF ( T OA ( 4) ) 90, 99, 90 

VERTICAL PLAVt OR CYLINUtR 

32 CALL FC2( TUA , TJ2, At , AN ,A 6, T , J3 P T ,Nj I M) 

HINBPT ,4) =H(N8PT,4) + TJ2( 10) 

HI 9 BP T ,6) =ri(NBPT, 5) + Tu2( 10) ATJ2 111 
IF ! L J! 5) .bu. 0) \j u U 33 
WR I TE ( 6, lU9) I TJA ( J) , J=i , 7) 

WR 1 TE ( 6, 1 19) ( T02 ( J ) , J=1 , 10) 

33 IF ( T J A ( 4 ) ) 90, 99,90 

LIojUlU FILM COOLInG, GkAVITY FLOw 

34 CALL FC3! TJA , T J 2 , At , AN , A S) 

FilNBP T ,*+) =Fi ( Nb P T , 4 ) + TO 2 ( o I 

FI! VdPT, 5) =H(NBPT, 5) + T 32! 6) * T02 1 1 ) 

IF (LJ(5) .Eoi.O) GO TO 35 
WRITE! 6, 1 09) ( TOW J) ,J = 1 ,7) 

WR I TE ( o,120) ( TU2(J) ,J=1 ,o) 

35 IF ( TO A! 4 ) ) 90, 99, 90 

FOR FUTURE SUBROUTINES 15 MOkE ) (STATEM. 3 6-7,38-9,43- 1,42-3,44-5) 
RAJIA T ION 
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c 

c 


c 

c 


c 

c 


c 

c 


46 CALL i< Al( IDA, Tu2fAEfAMfAS,T f JaPTfMUlM| 
H( M BH* T,o) =HI \I8PT,6) + TD2I 71 

HI NdP T,7) =HI M0P T, 7) + TD2I 7)*TD2 (1 ) 

IF I L D I 5 ) .EQ.O) 3D n 47 
wkI TEI 6, I 09) I TDAl J i , J = i t 7) 

WRITE <6, 126) (TD2(J),J=1 r 7) 

47 IF ( TDA ( 4) ) 90, 99, 90 

SPECIFIED COEFFICIENT 

48 N Ck = T DAI o ) 

GU TJ 149, 50, 49,31 ,52 ), NCR 

49 AF= AE 

GO TO 53 

50 AF=A^ 

30 TO 53 
5 1 A F = A S 
GO TO S3 

52 AF = TO A { 7) 

53 T U 2 ( 3)=TD2(2)*AF 

HI \ 8P T , 2 ) =H ( N 3 P T , 2 ) + T 0 2 ( 3 ) 

HI NBPT,3) =H< \ittP T, 3J+TD2I3 >*T02 ( 1 ) 

IF (LDI 5) . LQ . 0 ) 3U TO 54 
WR I T l ( 6,109) I TO A I J ) , J=1 ,7) 

WRITE! 6,127) ( TD2( J) , J=1 ,3) 

54 IF { TDA< 4 ) ) 90, 99, 90 

VARIA3LE TEMPERATURE FURCED CONVECTION 

55 N F F=N TF + 1 

IS! YP E =Tb A I 3 ) 

GJ TO I 56, 57) , I STYPE 

VAR .TEMP. DUCT FLOW, LIQUID 

56 CALL SNLI TDA , TD2, Ac, AN ,A S) 

HI N6P T, S) = TD 2 1 16 ) 

HI N8P T ,9) =Tj 2I 5)*rD2(3)*T02 ( 12 ) 

TBNDI NTF) = T3 2< l) 

TCMNTF) = HIM6Pr f 3) /HIN8P T,9i 
AFnIINTF>= TD2I 8) 

GF\ ( M IF ) = TD 2 ( 9) 

N3DEIMTF) =J3P T 
JXPT=NODc ( N T F ) 

JNT F{ JXPT )=NTF 

IF (LDI 5) .EU.O) GO TO 6U 

WR I TE ( 6, 109) ( T DA ( J ) , J-l , 7) 

WRITE! 6, 128) (TD2I J) , J = i, L6) 

GO TO 60 

VAK .TEMP. RADI AL StAL GAP 

57 CALL SN3I TDA, 102, AE, AN, AS) 

T b N D I N IF ) =TD 2 ( 1) 

HI ^0P T f b)i=TJ2( 19) 

Hi N BP T, 9) = To2 I 14)*TD2<3) 

TCFINTF) =H(iMdPT f 3) /HI NrtP T, 9) 

AFN I N TF ) = 102! 12) 

GFN ( N T F) = TD2 I 13) 

NiODtl N TF ) = J3 P T 
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JXPT = \IODE (NTF ) 

3NT F ( JXPT )=nTF 

FLUU SHEAR HEAT GENERATION 

JMRP = Tl> 2< 8)*T02(8)*n i ( JtlPT,2)-TDl (JBPT,2) 

Of ( MTF ) = T02< 20)*. OOOU3 9941 6*TD2 (4) *UMRP*T D L ( J B PT , 1 4 ) / T 02< 7 ) 

If <LD(5). ty.O) GO TO 60 
WRITE (6,1 09) { TO 4 ( J ) , J = L , 7 ) 

WR ITE ( 6, 12 9) ( TD2 ( J ) , J=L ,20) 

GJ TJ 60 

FDR FJTlRc SUBROUTINES (2 MORE) (STATEMENTS 58, 69) 

FLOW INITIATION ELEMENT 

60 IF (TO 2( 1 1 ) .EO. 0. ) GO TO 6L 
N TR =NJ TK + i 
N6R(NTR)=rD2I II) 

GO TO 62 

COMBINED Flow element 

61 IF (TO 2( 10) .£<J. 0. ) GO Tu 62 
NOR=NDH«*l 

REA31 5,101) NAT ( nDR ) ,NBT (NOR) 

WRITE(b,L32) JBP T,NAT( NOR) , N B T ( N UR ) 

NJ M ( NOR ) = J BP T 

62 IF (TOA( 4) ) 90, 99,90 

99 RETURN 
END 
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$ i BF rc SiY 1M 


SUBKG J FINE SIM 1 ( TDA , TO 2 ,AE ,AN,AS) 

Ht A T THAN SFER COEFFICIENTS FOR LIUUIDS IN DJCTS 

DIMENSION TD2I 20) , TDA I 7) 

CUM P J F fc REYNOLDS AND PR AND T L NUMBERS 
T0 2I 1 3 ) = 1D2I 5)* TD2 ( 6) / TD2 14) 

TU2( IN ) = FD2I 3)* TO 2 (A) /TD2 (2) 

LAMINAR jR TUkOULL.NT COEFFICIENT 
T U 2 1 i 1 )=TU2( b) /T.J21 7) 

TD2( 13 ) = TD2( 2) /TJ2( o) 

TD2I L 0 ) = T D 2 ( 1 7)**. 7 

IF ( TD2( 1 3)-2t>00. ) 1 , l ,2 

1 T 02 ( 15> = l .24* TJ2< Id )* ( ID2< 13) *T02t 14) *T02 < 17) ) **U ./3. ) 

GU TO ‘3 

2 IF I T J2< l 3 ) - 7000. ) 4,3,3 

3 TU2 1 15 )=. 0 2 3* TD 21 18)* ( T0 2( 13) **. 8) * (TD2 (HI#*. 4) *( l . + . 3*TD2< 19) ) 
o J TO 3 

4 TU2( 1 5 ) = . 5* ( ( 7000.- TO 2 113) ) A*500 . ) * l . 24* ( T D2 < 1 3 ) *T D2( 14 > *T 021 1 7 ) ) * 
2*1 l ./ 3 .) f .5* ( ( TJ2( 13) -2500. ) /4500. ) *. 023*T02 U 8 ) *< T D2 < 1 3 ) * * .8 ) *( T D 
32 ( 1 4) * * . 4 ) * { 1 . + . 3* TD2 ( 19) ) 

oJEFFICIENT*ARtA 

5 NCk =TD A ( o ) 

GJ TD ( c , 7, 6, 8 , 9) , N^ K 

6 AF=AE 

GO TD 10 

7 AF=AN 

Gu TO 10 

8 AF = AS 

30 TJ 10 
4 AF=TDA { 7) 

10 TD21 16 )=AF* TO 2 { 1'3 ) 

RtT JRN 
END 
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tlBFTC SN2M 


S LRR 01 J T IN E SN? < TOA * TD2 , AF ,AN,AS) 

H F AT TRANSFER COEFFICIENTS FOR FLUID MOVING RELATIVE TO RADIAL 
SURFACE SFPFRATE BOUNDARY I A YEP S 

OIMFMSION TO 2( 20 1* TO A ( 7 ) 

REYNOLDS AND PR AND TL NUMB FR S 
T C? ( 1 1 ) = TD2 ( f 1 * TD2 ( 6 ) * TD2 ( 6) /TD2 (5 ) 

TC2< 1? )=< T07(3)*TD2(4) /TD2C2) I **<l ./3. ) 

LAMTNAR OR TLRBULENT COEFFICIENT 
T D?( 15 ) =T D 2 ( 2) / TO 2(6) 

TD2( 1M=f T0?<7)/TD2(6> )**.1 
IF (TD2( 1 1 1—15 ROOD. ) 1,1,4 

L AM INAR 

1 IF (T02C 10) .FQ.D.) GO TP 2 

T C 2( 17) = * SOI F 5— *65C8*TD?(7)/TD?(9) 

GC TO 3 

2 TF2< m = , 4SE 15+.65C8«TD2( 7) /TQ2 (9) 

1 T P7( 1 1 ) = • 57 4*Tf)7( 15I*TD2( 16)*< TD2( ll)**« 51 *TD7 < 17) / TD2< 17) 

GC TO 7 

TURBO! FMT 

4 IF ( T97( 10) .FO.O.) GO TO 5 

TC?( 17 ) 4E C<: R-.157 7 8* TD? (7) /TD2 (9) 

GO Tn 6 

5 T T 2 ( 1 7 )=. 5144?+. 35778* TD2 ( 7 ) /T 02 (O ) 

6 TD7( 1 I )=.01B26*TD2( 15)*TD?(16)*( TD?( 11 ) **.3) *TD2(12)/TD2( 17 1 

CPFFFIC IEMT*ARFA 

7 N CR = TDA ( 6 ) 

GC TO { 8, 9, R, 10,1 l ) * NCR 
fl A F= A F 
GO TO 12 
9 A F = AN 
GC TO 12 

I P A F = A S 

GO TO 12 

II A F= T DA ( 7 ) 

1? TF2( 14)=AF*TC?( 11) 

P ETURN 
FND 


no n n n o n 


$ T R FTC SN3M 


SUBPOUTTNF SN3< TDA * TD? , AF , AN , A S ) 

C 

H FAT TP AN SF FR COEFFICIENTS FOR ROTOR AND STATOR WITH SMALL 
CLEARANCE MFROPf) BOUNDARY LAYERS 

OTMFNSmN T0;( ?0),T{?A{7) 

REYNOLDS AND RRANDTL NUMBFRS 
T 0 ?( 15 ) =T D 2 ( n*T02(6)*T02<6)/TD2(5> 

TC?( ]6)=< TO 2(3)* TD?(4)/TD ?(?))**(!. /3.) 

t £*TNAR OR HRRULFNT COEFFICIENT 
T r?t 1 7 ) = 3 66 • ?4* ( ( TD2 ( R ) /TD2 { 7) ) ** ( ] 0. /9. )) 

TF ( TD 2 ( l E ) - TD ? ( 17) ) 1*1*? 

1 T P? { 13 )=? .*TC?< ? )*T0?( 16) /TO?< 7) 

rr?( ? 0 )=i . 
or TO 3 

? Tr?(18)=*0?S7*(T02<?) / T D 2 ( 6 ) )* ( (TD? (6 ) /TD2 C7) > ** <1 ./6. ) ) *TD2( 16 ) *( 

IT r?( 15 }** .8 ) 

T C? ( ?0) =*000 5 26 6* ( TD 2 < 17 ) ** 1 . 5 ) * < ( TD2 (7) / TD? ( 6 ) ) ** ( 5 ./ 3 * ) ) 

C 

C CrFFFICIFNT*ARFA 

3 N CR~ T 0 A ( 6 ) 

GC TO ( 4, 5* 4,6, 7) f NCR 

4 A F= A F 
CP TO R 

5 A F= AN 

no rn « 

6 A F= A S 
GP TO 3 

7 A F= T C A ( 7 ) 

8 T T? ( 1R ) - A F* T P 2 ( IB) 

R FTURN 

END 
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r> n r> n n n a n 


tlBFTC ^NAM 


SURPniJTINF SN4< TDA,TD?.AP ,AN,AS) 

H F AT TRANSFER COEFFICIENTS FPR FLOW BFTWFEN CONCENTRIC CYLINDERS 
INNFP ROTATING* NO AXIAL FI OW FF FF r T ON BOUNDARY LAYFRS 

D IMF NS I ON T 0 2 ( 20) • TDA (7) 

TIJRBUI.FNT CPFFFICIFNT 

T R?( 11 )=( TD2(3)*TO?(4)/TD2(?) )**(l-/3.) 

TC?< 14>=TP2{ F)*.5*( TO?( 6) + TP? ( 7)) 

TO?( 15 )=T0 2( 14 )*T0 2< 14 ) + TD? < <3) *TD2 (9) 

T F 2 ( )4) = S CRT ( TO 2 ( 15) ) 

Tn?( I 7 ) =T D 2 1 O-T02I 7) 

T R?( 10) =2 „*TC?( 17)*TD?< 16) /T02 (5) 

tr?( 12 )=( TO 2(2)/ TO 2( 17) ) * ( . 01 5 * ( TO 2 (1 0 ) * * . B ) *T02 i 1 1 ) *{ (T0?(6)/TD2( 
17 )>**.46) +*r c?* <TQ2< 14)**<?./3.) ) * TD? <1 7 ) *T02 ( III /TO? ( 5 ) ) 

C C F F P I C IFNT*ARFA 
N fR= TO A ( 6 ) 

OP TP (1*2* 1*3*4 ) f NCR 
t A P= A F 
OP TP 5 

2 A F = AN 
or TP 5 

3 A F = A S 
OP TO 5 

A A F = T P A ( 7 ) 

5 T r?( 1 3 ) = A F* T P 2 ( IP) 

R FTURN 
^ N 0 



r> n on non 


i T B FT f SN5M 


S UPROOT INF SN5( TDA * TD? . AE ♦ AN * A S ) 

HFAT transfer coefficients for flow gvfr flat plates 

n TMFNS ION TD2( 20 >♦ THAI 7) 

REYNOLDS AND PRANDTL NUMBERS 

rn?( P ) = TO 7 ( 7 )*T02( 6) /TD2( 5) 

rr2( <3 1 = (TD2( ^M f T02«4)/TD2(?) )**n . /B. ) 

LAMINAR OR TURBULENT COEFFICIENT 
T F( TO?{ P ) — *SC COO0. ) Mi? 

1 TT2( 10 )=.?2?*TQ7(2)*(TD2<8)**. 5)*T02t9) /T02I6) 

CP in B 

7 TP2C 10)=* 0 28 f*Tr>?< 2)* <TD2< 0) **. 81* TD2 (9) /TD? 16) 

C C F F F I C IENT^AP FA 
B M CR = TO A f 6 ) 

00 TO ( A* 5* 4,6, 7), NCR 
A A F = A F 
or T 0 R 
6 A F“ AN 

or to p 

6 A F= A S 
or TO R 

7 A E= T C A ( 7) 

P T C?( 11 ) = A F * T C 7 ( 10) 

RFTtJRNj 

FNO 
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-n r> n n n n ^ n n 


SIR<=TC fCIM 


SURPPUT IN F F C 1 C TOA « TD? * AF * AN * A S * T , J6 P T * NO I «) 

HEAT TP AN SF FP COFFF tC I EN T S FOR HOR 1 1 . C YL T NOE R FRFF CONVECTION 

H IMFNSTON TCA< 7 ) * TO 2 (20) . T(ND[ M) 

TT?( 1?)=ABS( 10 ?< 1 )-T( JBPT) ) 

TC2( 11 ) = „ 5*< TD2( 1 ) + T< JRPT) ) +4 5 9. 7 
T C?( 14 ) = T D ? { 2) /T 02< 6) 

PPANDTL ANO CRASHDF NUMBFRS 
T C? ( 7 ) = Tf) ? ( 3 )*T02(4 )/TD2(?> 

T F 2f ft) =T0 2( 6 )* T0 2< 6)* TO? ( 6 ) * 32 * 1 74* T02 (12) /(TO? ( 5 ) *T n2 ( 5 ) *T 02 ( 1 3 ) ) 
T C 2 ( 1 l )=T0?( 7)*TQ2( 8 ) 

1 AMTNAR OR H.RBULANT COEFFICIENT 
T F ( Tf)2( 1 1 )- .IF + 10 ) 1*1,? 

1 T r?( 9)=.5T*TC?t 14)*T02( 11 )**(! ./4* ) 

GO to 3 

? T r?( <5 J = . I Q * T C ? ( 14)*TD2(U)**(l./3.) 

CP FFFIC IENT* AR FA 

3 N f R=TD A { 6 ) 

or TO ( 4* E. 4*6* 7) .NCR 

4 A F = A F 

or TO 8 

6 A F= AN 
Of TH ft 

6 A E = A S 
Of TO ft 

7 AE=TOA( 7) 

ft T C?( 10 )=AF*TC?( 0 ) 

R FTUPN 
FP^D 
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M B FTC FC7M 


S URROI jT IMF FC? ( TQA , TD2,AE , AN ,A S , T , JBPT , NOT M) 

H f A T TRANSFER COEFFICIENTS FOR VFRT. SURFACE FREE CONVECTION 

n IMFNS ION TC fi( 7) , TO? i?Q) ,T(NDI H) 

TC2M2)=ABS(T0 7(1 }*T ( JBPT ) ) 

T T 7 1 1 3 ) = • S* ( TO 7 ( 1) + TIJBPT) )+45 9,7 
T r 7 ( 14 )=TD?( 2 ) / TO 2 ( 6 ) 

PPANOTL ANO FRA SHOP NUMBERS 
T D? ( 7)=T07( 2 ) *T02 ( 4 ) /TD7( 7) 

T P7I 8 ) =TO 7 ( 6 >*T0 2< 6)*TD2{6)*32.174*TD?f 12 > /( T07( 5) *T D2 ( 5 ) *T D2 ( 13 )) 
t r 7 { ii ) = t n 7 { 7 )*to 2( a) 

1 A MINAR OR 7 LRBIJLA NT COEFFICIENT 
I F ( T07 ( J I )-.i*= + 10> 1,1*7 
1 TC7(9)=.44*TC7( 14)*TD2<1 1 )♦*(!• /4 • ) 

GO TO 3 

7 Tr?<Q) = *l3*Tr7UM*TQ?Ul)*Ml./3.) 

C CfFFF [C I E N 7 # A° FA 

3 m CR = T 0 A ( A ) 

GO TO ( 4* 5,4,6* 7), NCR 

4 A F= A F 
GO TO 8 

6 A F= AN! 

GO TO 8 

6 A F = A S 
GO TO 8 

7 A F= T OA f 7) 

8 TC2( 10 ) = A F* T T7 ( 9) 

R FTlJRN 

FNO 
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n r> n n n 


$ F B FTC FC3M 


SUBROUTINE FC3( TDA , TD? * AE * AN * A S ) 

H F AT TRANSFFP COEFFICIENTS FROM HORIZONTAL CYLINDER TO L I QU 1 0 
FILM ( GRAVITY FLOW) 

OTMFNSION TO ?( 20 ) * TD A ( 7 I 

T C2f 5 ) =TD 2 ( 2 )*T02(?) /T02( A) 

N CR = TO A ( 6 ) 

GO TO (1,?,1,3,4 ) • NCR 

1 A F= A F 
GC TO *5 

2 A F= AN 
GC TO S 

3 A F- A S 
r,n TO 5 

A A F= T CA ( 7) 

*5 rr?( 6) =AF*TO?( F ) 

R FTIJRN 
FND 


TC RA1M 

S UPROOT IN E R AH TOA*TO?,AF f AN * A S * T f JB PT f NDI M ) 

R AOI AT ION ( THERMAL l-GRAY BOOY 

n IMFNS ION TC A( 7) ,TD2( 20) f T( NDI M) 

T AF = TD2 ( 1 ) 

T A S = T ( JBPT ) 

T C 2 ( fi 1 = 1 • — TD 2 ( A ) 

T 02 ( 9 ) = 1 •— TO 2( 5) 

T F 2 ( 10 )=T AE*TAF*TAE+TAE*TAF*TA S+TAE*T AS* TAS + T AS*T AS*TAS 
T C 2 ( 6)=A*B3F-13*TD?( A)*TD2(9)*T02 <?)*T02(10)/(l.-<Tn?(2)*TD2(3)*TD 
1 2 ( R ) *T D?C 9) ) ) 

Nf C R = T D A { f> ) 

GO TO ( 1, 2 v l • 3 * A 1 ♦ NCR 

1 A F= A F 
GC TP 5 

2 A F= AN 
GO T D *5 

3 A F= A S 
GO TO S 

A A F=T CA ( 7) 

*5 T T? ( 7 ) = AF* T 0 2C 6) 

RFTURN 

FND 



$ I BFTC oiTM 


SUbROUTI NE JTINBP T.NB JIHtNOI M , JB PT , H , SUMT , SUMH, TD1 , T BN 0, T, CD, NVOIM 

1, M IT, NODE , JMTF ) 

C 

C COMPUTE BOUNDARY HEAT FLUXES 

C 

DIMENSION HI NBQl.M , 9) , S UM T ( NBDI M ) ,SUMH I NBDI M) »TDL(NQIM, 19) 

DIMENSION TdNJl NVOIM) , I (NJIM) ,CD 110) ,TS(130) , JNTFI NDIM) 

COMMON /8LKI/J3 NT < 4 00) , RAH (2 00) , N TF , NOR , NTR , OF I 95 ) , NOR 

COMMON /BLK2/TCF 19 5) ,S F N I 95 ) , AF N ( 95 ) , NBR 1 9 5 ) , NM I 5 ) , NAT I 5 ), N8TI5 ) 

C 

101 FORMAT I IH0.2X, 39HCJM3 I NE J FLUV-.S ELEMENT, COMBINATION NO., I 2, 3< , L2H 
2CJNDUC T.NU . = , I 3,3 X, 10H3NJRY.NO. = ,1 3 , 3 X, 31 HPREV I OUS ELEMENT VAR .T EM 
3P .NuS . =, I 3, 2 X , 3HAN0 , 1 a ) 

102 FORMAT I 1H0,2X, 1 9HC UNO UC T. ELEMENT NO . , I 4 , 3 X , 2 OHV AR. T EMP . ELEM ENT NO. 

2, 14, 3X , 10HFL0W CODE =, 1 2) 

C 

SUMTI \| BP T ) =0 . 0 
SUMH(N3PT)=0.0 
I IYPE = TD1 ( JdP l , 17) 

IF < I TYPE .NE. 7) SO TO 1 
Go TO 99 
C 

c specifies flux 

1 IF ( H I NBP T , 1 ) . f J . 0 . )' jO TO 2 
SUM T( N 3P T ) =H I NB P T , l ) 

C 

C FORCED CONVECTION 1 1 NCL. SPEC I FI ED COEFF. ) 

2 IF (ril.NB? T,2).EO.O. ) GJ TU 3 
SUM HI NBPT ) =H(NBPT,2) 

SUMTI NBPT ) =SU(MT( NBP T) H ( NBP T , 3 ) -H ( NBPT ,2) *T( JB PT ) 

C 

C FREE CONVECTION 

3 IF (H(Na 3 T,4) .EO.O.) Si) TO 4 

SUM HI M BP T ) = SUMH I NB P T ) M(NBPT,4) 

SUMTI NBPT ) =^UM T I NBP T) i-H I NBPT ,5) -H (NBPT ,4)*TIJBPD 
C 

C RAO IA T IUN 

4 IF IHI NBP T, 6) .to. 0. ) SU TO 5 

SUM HI N BPT)=SUMH(N3PT)+HINBPT,6) 

SUMTI NBP I )=SUMTI NBPT) ■MINBPT ,7)-H(NBPT ,6) *1 ( JBPT) 

C 

C VARYING TEMPERATURE 

3 IF (HI NBP T ,3) . EG. 0. ) GO TU 99 

NTF=NTF+I 
NP T=0 
L H\ =0 
w F 1 =0 . 

SUM Hi NBPT ) =SUMH( NBPT) +HINBPT ,d) 

IF (KANIN3PT ).E0. 0) SU TO 6 

SUMTI N BP T ) =S JM T ( Nd P T) a H ( NBP T , 8 ) * I TBND ( NT F ) -T I J B PT ) ) 
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6 JNPT=AFN< NTF ) 

JlHN= GFN( NTF ) 

T S = 0 • 

IF (JNPT.EQ.Q) GJ TO 7 
NPT =J BNT( JNPT) 

7 IF (JLHN.EO.O) GJ TO 3 
LHN-JNTFI JLHN ) 

6 IF (.vlM(NDR)-JBPT) 14,9,14 

COMBINED FLOWS 
9 JNFR-NATi NDR ) 

JNFS-NBK NDR) 

IF ( JNFK. EQ.O) GJ FO ID 
NFR=J BNT ( JNFR) 

JAT=JN TF ( JNFR) 

GJ TJ 11 
LO M F R = 0 

TBMD( NFR) = U. 

H(mPK, 9 ) = 0 • 

11 IF (JNFS.Eg.O) GJ TO 12 
NFS-JBNTl JNFS) 

JB f=JN TP ( JNFS) 

GJ TJ 13 

12 NFS-0 

TBi\U( NFS) =0. 

Hi, NFS , ^ ) = 0 . 

13 IF ( { vi IT- l).EU.O) *vRi TE ( 6 » 1 01 ) NOR , J BP T , NT F , J AT , JBT 

STARTING TEMP. FOR COMBINED FLOWS ELEMENT TEMP. CALC. 
COMbF =Hl NFR, 9) +HCNF S, 9) 

PNFR-H (INF R, 9) * I bND( NEK ) /COMBF 
PNFS = H(NFS, 9)*TdNl*f Nr S) /COMBF 
STbMP=PNFK«-PNFS 
NUR=N JR + 1 
GO TJ lb 

START ING TEMP. ORDINARY FLOW ELEMENT 

14 NTk-NTR+1 
NEK=NBi<(N TR ) 

IF { ( N IT- i ) • b J • 0) WRI TE C 6 » 1 02 ) J6P T , NT R ,NBR{ NT R ) 

IF (LHN.NE.O) GO TO 17 
GJ TJ ( 15# 1 6 ) fNER 

15 STtMP -CD ( 6) 

GO TO 18 

16 STEMP-CJI 7) 

GJ TJ id 

17 STEMP=TBMD(LHN) 

18 IF UAN(NBPT) .ME. 0) GO TJ 99 

CALC . TEMP .CHANGE 
HI\lC=TCFlNTF)/CCH 9) 

IF (JF(NTF).EJ.O. ) GO TO 19 
QFl=wF(NfF)/iHCNBPT t 9)*Ci3i9n 

19 NC D=C 0(9) 
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on no 


IF (NPT.Nt.O) GO TO 22 

TEM P .CHANGE-NO ANN, IN Ad OK 

DO 21 INC = 1 r NC 3 

IF ( INC.NE.l) GO TO 20 

TS( INC )=STEMP+HINC*(T< JB P T) -STEMP ) +QF i 
G] TJ 21 

20 TS( INC J = T SI IMC— II +H l NC* ( T( J8 P 1 1 - TS ( I NC— I ) H-QF1 

21 CONTINUE 

7 BND( :N TF)=TS( INC) 

GO TO 25 

TEMP.CHANGL-ANNULAk N A 3 U i< 

22 HANC=H( NP T f 3)/<H(NPT#9)*CD(9)) 

DO 24 INC =1 f NC J 

IF UNC.NE.l) GO TO 23 

I Si INC )=S TEMP +hl NC* (T( J bP T) -STE MP) *H ANC* f T ( JNPT ) -STEMP) +QF I 
GJ TJ 24 

2 3 TS( INC )=rS(INC-l)+HlN:*(T{JBPT)-TSUNC-U ) +HANC* ( T (JNPT ) -f S ( INC-l) 
1)+0F1 

24 CONTINUE 

TdNDIN TF ) =TS{ INC) 

TbN D( NPT ) =T3ND I NTF ) 

K A ^ ( NP T ) = i 

2 5 SU4 T{ NdPT )»SUMT(NtiPT) + -i (NtiPTvd) * ( TtJND (NTF) -T ( JBPT) ) 

99 *ETU*N 
END 
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m n n n n -1 n n 


t T B FTC TAT SM 


SURROUT TNF T AT S ( TO 1 ♦ JPT , ND I M ,NBNO, NBD I M , Of T , H , R 0 ♦ L n , CE , NK ) 

OVFR-RFLAXAT ION AND NEW TE MPFR A TIJR F CALCULATf C-N 

DIMENSION T01(N0IM,19) ,CE (NO JM,4)»LO(lO) , BO (5 ) 

DIMENSION NB NO ( NBD I M > , 0 ( ND l M ) . T ( ND T M ) »H ( NflDT M 1 9 ) 

CrwiMON/RLK l/JBNTf 390 ) , KAN ( 190) * NTF , NOR , NTR , OF <9S ) , NOR 

101 FORMAT! 1HQ, 5 X, 5PHERR0R IN B NOR Y AND CONI). NO CnRPES P . ( T ATS CHECK OF 
2 I TF» AT TON ,16,1 3H ) COND# El FMENT t 14,2 5HG IVES BNDRY NO OF ELEMENT, 14) 


H SUM = 0 . 

FAST CONDUCTANCE 
N F= TO 1 { J PT , 5 ) 

TE (NF) 2,1, 2 

1 CF( JPT* 1)=0.C 
CP TO 3 

2 C F( JP f , l ) “ 1 • / ( TDK JP T , 4 ) / (2.*TDl ( JPT* 14) *T01 ( J PT * 1 8 ) ) + T 01 ( N E, 4 ) / ( 2 
U*Tni{NE, 14)*T0 1(NE'18) ) + TO 1 ( J PT * 1 0) ) 

NfRTH CONDUCTANCE 

3 NN=TPU JPT, 6 ) 

T E ( NN ) 5* 4* E 

4 rF(j°T f ?)=o.r 
or TO 6 

5 CF( JPT, ?) =1 . /( 2.*TD 1 ( JPT , 3 ) /(TD1 ( J P T , 1 B ) * { 3 • *TD1 (JPT ,L5)+T01(JPT, 1 
16 )) )+?,*Tni(NN, 3) /( TD1(NN* 18)*(3 # *TD1 <NN,16)+TD1 ( NN , 15 ) ) )+TDl(JPT, 
21 1 ) ) 

WEST CONDUCTANCE 

6 NW=TC1 ( JPT* 7 ) 

T F < NW ) B, 7, P 

7 CE( JPT, 3 ) = 0 • C 
00 TO 9 

B C F( JPT, 3) = 1 . /( TDK JP T, 4) /( ?**TD1 f JPT , 14 ) *TD 1 (JPT ,1 8) )+TDL( NW*4)/(2 
1. *TD1 ( NW* 14 )*TI)1 (NW, IB) ) + TDl ( JPT, 12) ) 

*C 

C SOUTH CONDUCTANCE 

9 NS- TDK JPT, f) 
r E (NS ) 1 1 , 1 C, 1 1 

in C El JPT, 4 ) - o . c 
DO TO 12 

11 CF(JPT,4)-l./( 2 .*T0 1( JPT *3 ) /( TO If JPT, 18) *(3. *TD1 ( JPT, 16 )+TDl( JPT, 1 
15 )))+7.*TDl (NS,3) /< TD1 (NS, 18 )* (3.*T01 ( NS, 15 )+T01 (NS, 16) ) )+T01(JPT f 
212)) 

C 

1 2 NTD=T01( J»T, 17) 

IE(NT0.E0.7) On TO 99 


67 



C COMPUTE RES TCUAL « OVER-RFLA* 

R FS= CE ( JPT * 1 >*<T(NE)-T{ JPT) ) +C E ( JPT , 2 ) * ( T i NN ) -T C JPT) ) *CE ( J PT , 3 ) *< T 
If NW)-T(JPT) |4CF(JPT#4)*(T<NS)-T( JPT) )+Q( JPT) 

CS UM=CF{ JPT. 1) +CF< JPT,?) + CE ( JPT,3)+CE< JPT.4) 

c 

WRES=RD{ 1 )*RFS 
C 

f COM PUT F NFW TFMPFR A TURF 

l F (Of JPT) ) 13, 15, 1 3 

1 3 N ePT=J8MT( JPT) 

J PPT=NBND<N8PT) 

IF ( JBPT.FO^JPT) GO TO 14 
W R T TEl 6 . 1 Cl } NTT.JPT.JBPT 
MK - 1 

GC TO 99 

1 4 MSUM=H(NRPT*2)+H(NBPT.4)+H<NBPT,6) +H(NBPT,8) 

! 5 TfJPT)=T(JPT ) + WRES / ( C SUM+HSUM) 

C 

99 RETURN 
FK ; 0 
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cue 


* T B FT C FITA 


SUBROUT INF FIT?! TO 1, NO IM, JPT, AL PH A, THFXP) 

FIT FXPANSICN COFFF 1C I ^NT FROM POL YNOMI NAL CURVE 

0 IM FUSION TO If NDIM.19) ,COS( 5) .CALC (6) .ALPHA! IS ,11 ) , TPEX P ( NO I M ) 

C 

101 FORMAT! 1H0, PCX, A5HTFMPERATURF EXCFEDS CURVE FIT LIMIT FOR ALPHA) 
C 

T FMP= TO) ( JPT, 19) 

MAN = mil JPT, 9) 

on 1 NCOS 1.5 

1 CCSINCQ)* ALPHA { MAN, NCO) 

1 F ( TFMP-C.OS! 1) ) 3,3,2 

2 WRITE! 6, 1 Cl) 

Of TO 9 

3 NOP = COS! 2 ) 

NOP 1 =NOP + 1 

00 A N CQ= 1 » NCP l 

1 r = 5+NCQ 

A CAIC!NCO)= ALPHA ( MAN ,LC 1 
TFIC0SI3)) 6,5,6 

5 Y = CAI C.I 1 )*CCS! A ) 
on TO R 

6 y = CALC! NOP 1 ) 

TFM1 = TFMP*C0S!5) 

00 7 NPY=1,MQP 

J? = NOPl-NPY 

7 Y= CALC! J7 1 +Y*TEM1 
Y = Y*COS ( A ) 

B TFFXP! JPT ) = Y 
C 

9 R FTURN 
F NO 



^ r> n r> r> r> 


tTBFTC CONCHC 


SUBROUTINE CONC FK( JPT.NDI M, fi.)l ,ERR , LO , NK , T) 

CHFCK FOR CONVERGE NIC. F OF OVER-RELAXATION CALC. 

n TMFNSTON TDl{N0IM,19),LD(10),T<NniM) 

101 POR^ATl 24X, I*, 415X.E14.7)) 

0 !FF=T( .JPT )-ini ( JPT, 19) 

IF { ABSIOIFF l.LF.FPR > GO TO 1 
NK = 1 

l IF { L0< 4) .EC .0 ) GO TO 9 

WRITP(6,101) JPT.TI JPT) , TDK JPT.19) ,OIFF ,ERR 

9 RFTUPN 
FNO 


: CFFRMM 

S URROIIT INF REFORM! JP T,NDI M, mi ,R0, ID, THEXP) 

COMPUTE FRFE THFRMAL EXPANSION Al PHA*DELTA TEMP 

n IMFNSION TCKNDIM. 19) *THFXP(NDTM) » RD ( 5 ) , LD C 10 ) 

C 

101 FORMAT MHO, 1CX, lfiHFl FMFNT CPNQ . Nf)= , 1 4 ,3 X ,6H AL PHA= , El 3. 6, 3X , 5HTEMP= 
?, F1P.5, 3X, 1 4 FAL PHA* ( Tl-T0) = .E13.6) 

C 

T IN =R Dt ?) 

0 FI TF-M = TO 1 (JPT, 191-TIN 

A=THFXP( JPT I 

A TT =THFX° f JPT) *OELTFM 

T FFXP( JPT ) =ACT 

IF ( I 0( 10 ) .FC.O) GO TO 99 

WPTTFl ft.101 ) JO T, A. TDK JPT, 19) .THEXP(JPT) 

C 

99 RFTIIRN 
T NO 
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Figure 8. - Flow chart for main route SEAL2. 
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Figure 10. - Flow chart for subroutine QT. 








APPENDIX B 


CHANGING OF BOUNDARY CONDITION COEFFICIENTS 

The program is designed to make the addition or substitution of equations for heat 
transfer coefficients relatively easy. The coefficients, the products of coefficient times 
boundary area and the products of coefficient times boundary area times bulk fluid (or 
environmental) temperature are calculated in subroutines SN1, SN2, SN3, SN4, SN5, 
FC1, FC2, FC3, and RA1. 

The varying temperature coefficients, at present, use subroutines SN1 and SN3. 
Other calculations for the varying temperature boundary conditions are performed in 
subroutines CQN and QT. 


Substitution of a Different Coefficient 

To substitute a new expression for one of those presently used, the appropriate sub- 
routine is rewritten using the present words in the TD2 list for the output (coefficient, 
etc. ). This assumes that the same properties or at least the same number of properties 
are required. If not, the number of items to be read in must be changed in subroutine 
CQN. Also the TD2 words used for output must be changed in the coefficient subroutine 
as well as in subroutine CQN. The appropriate formats of subroutine CQN must also be 
changed to reflect the TD2 changes. 

The TD2 list of a varying temperature condition must be compatible with the TD2 
list of the coefficient subprogram used. Thus, if subroutine SN1 is changed to require 
more input data, subroutine VT1 (see table II) must also be changed. The output, which 
is common to both, must use the same new TD2 words. 

TD2 word changes for a varying temperature subroutine must be made in subroutine 
QT as well as in subroutine CQN. Wherever a TD2 word (e.g. , TD2(15)) appears in the 
varying temperature sections of subroutines CQN and QT, the new word must be substi- 
tuted. 

Also the two indicator flag words (TD2(10) and TD2(11)) must be the same for all VT 
(varying temperature) conditions. The calculations involving the flags require this. The 
sections in subroutines CQN and QT handling the temperature change along the flow and 
the effect of flow combination use the flags to control the flow in the subroutines. If one 
or both of the flag words are changed for one VT condition, it (they) must be changed for 
all conditions and the new word(s) must be used in subroutines CQN and QT. 
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A word of warning should be added. As implied previously, the TD2 list is used for 
both input to and output from the coefficient subroutines. This table is dimensioned for 
20 items. If more are needed, the dimensions in SEAL2, CQN, and the coefficient sub- 
routines must be increased. 


Addition of New Coefficients 

If an additional boundary condition is desired, a new subroutine must be written. A 
varying temperature condition involves the precautions mentioned under substitution. 

Also the TD2 list capacity must be remembered. 

Also, for an added condition, a section must be added to subroutine CQN. This sec- 
tion will be analogous to the subtype sections already present. The products of coeffi- 
cient times area and of coefficient times area times bulk fluid temperature must be 
stored in the correct places in the conductance matrix. An output format for the TD2 list 
must be composed. Unused statement and format numbers are provided for this purpose. 
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APPENDIX C 


COMPUTER PROGRAM APPLICATION TO MAIN SHAFT SEAL 

Seal Description 

This appendix is presented as an application of the computer program to a specific 
seal design. Figure 11 shows a schematic of a mainshaft seal being developed for gas- 



turbine engines. This is a face-type seal with a self-acting gas bearing which acts to pre- 
vent rubbing contact when deformation of the sealing faces occurs. 

In figure 12(a) a modified version of the nosepiece assembly is presented and the 
boundary conditions are shown. This modification is for the purpose of concentrating on 
the mechanics of the program application. 

Figure 12(b) shows the element schematic including conduction numbers and dimen- 
sions. A portion of the seal plate is included since it influences the temperature of the 
airflow through the seal gap and therefore has an effect on the nosepiece temperature. 

Also a portion of the piston ring and its holder, which form a conduction path to the nose- 
piece carrier, are included. 

The numbering sequence starts with the nosepiece assembly since this is the region of 
interest for this problem. The desired printout is therefore available at the beginning of 
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(a) Modification of and boundary conditions for nosepiece assembly of seal depicted 
in figure 11. 



(b) Dimensions and conductance element numbering for nosepiece assembly depicted in figure 12(a). (All 
dimensions in inches.) 


Figure 12. - Example problem. 




(c) Boundary element numbering and boundary heat transfer conditions for nosepiece 
assembly depicted in figure 12(a). 


Figure 12. - Concluded. 


each section of output, and the remaining elements can be ignored. The piece by piece 
numbering is not necessary; any convenient order could be used. 

Figure 12(c) shows the boundary numbers and the subroutines used to calculate the 
heat transfer coefficients. Here again, the sequence is a matter of convenience and the 
nosepiece assembly is numbered first. Numbers 1 to 11 are assigned to the region where 
concentric cylinder flow (SN4) is assumed. Conduction element 56 was ignored since only 
a relatively tiny part is a boundary of the piece. Numbers 12 to 41 were assigned to the 
liquid film cooled (FC3) boundaries. The varying temperature duct flow (VT1) is num- 
bered next. Since conduction element 12 already has a boundary number (7), the sequence 
continues with 42 assigned to conduction element 11. The method used is an approxima- 
tion to avoid the complexities involved if conduction elements 18 and 14 are included as 
should be done. For the varying temperature seal gap (VT2) flow, numbers 48 to 50 were 
assigned to the nosepiece elements not already given boundary numbers. 

For the seal plate, the east face of conduction element 92 and the north face of 110 
were assigned coefficients. The south face of 92 was assigned a heat flux. These values 
were obtained from a thermal analysis of the seal shown in figure 11. Boundary numbers 
53 to 60 were assigned to the seal plate elements exposed to seal gap flow. The remain- 
ing boundary numbers, 61 to 72, were assigned to elements given specified temperatures 
obtained from the previous analysis mentioned. 

The title and parameter cards (Items 1 to 5 under Computer Input in appendix A) were 
filled out as indicated in the first section of the printout. Also in this section, Items 6 
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and 7, the thermal conductivity and expansion data, are indicated. 

Next are the conduction element cards (input Item 8). These are shown in the second 
section of printout (conduction element table). Words 1 to 8 of each card are obtained 
from the information in figure 12(b). In word 9, three materials were used. The nose- 
piece was material number 1. The piston ring and its holder were material number 3. 

The remaining parts were material number 2. The last four words are contact resis- 
tances. 

It should be noted that each contact resistance must appear twice since the program 
ignores neighbor contact resistance in computing heat fluxes. Thus for elements 42 
and 46, the contact resistance appears in word 11 (north contact resistance) of element 42 
and in word 13 (south contact resistance) of element 46. Contact resistances between 
parts were calculated using a "conductance coefficient" (h c = 0. 2) for conductance across 
the interface. The resistance between elements 90 and 60 was the difference between the 
resistance of the air gap and the resistance of an equal volume of material number three 
(element 90 uses the same mean radius and radial width as element 89). The resistance 
between mismatched elements (elements 52 and 54 were changed in axial length to approxi- 
mate the actual dimensions in the region) was calculated as illustrated in figure 13. It was 

ii. 

assumed that the larger i in element (e. g. , 52) has a conduction path length of d' rather 
than Ar/2 as in figure 5. The area was based on AZ m rather than AZ^. The contact 
resistance was the difference between d'/27rr m A Z m and Ar^/4jrr m AZ^. 

The initial temperature estimates (input Item 9) follow the conduction cards. These 
values were obtained from a previous analysis in order to reduce the computer time re- 
quired. However, any reasonable estimates can be used including the use of one temper- 
ature for all elements. 


r 2 


r 2i 



Figure 13. - Outline of calculation of contact res istance betwee n 

elements of different axial dimensions. d' = Va 2 + (Alj/2) 2 ; 
At m - (AI, + At 2 >/2. 
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-a 

CD 


The boundary heat transfer cards (input Item 9) are assembled element by element. 
The boundary numbers determine the order of the sets. A set consists of the indicator 
card for the element followed by the index card, plus any others required, for each ’’side" 
of the element that is a boundary. The indicator cards do not appear in the printout but 
the order of the other cards for each element is indicated in the section of printout follow- 
ing the initial temperatures. This section, as well as the conduction element table, are 
not standard output items but are included to indicate the cards required. Some calculated 
values (Reynolds or Prandtl numbers and the following items for each element side) are 
included in this printout. 

In preparing the varying temperature boundary cards, the region at which the flows 
joined required the use of approximations because the program permits only one varying 
temperature condition per element. Boundary element 47 (fig. 12(c)) was assumed to be 
affected only by the duct flow. For boundary elements 49 and 50, special areas (code 
number 5) were used with the sums of the appropriate faces as the area values. Ele- 
ment 56 was made the combined flow element with 47 and 55 (conduction numbers 19 
and 98) as the upstream elements. 

The film cooled (FC3) boundary was another area where approximations were in- 
volved. It was assumed that flows originated at the northwest corners of boundary ele- 
ments 16 (two 0.018 lbm/sec), 24 (two 0.018 lbm/sec), and 35 (two 0.012 lbm/sec) (see 
fig. 12(c)). The flows were assumed to be either toward or parallel to the axis except for 
elements 20 and 31 to 33. The arrows outlining the region indicate this. Elements 20 and 
31 to 33 had flow around the circumference and were the only elements where the equation 
properly applied. For the other elements, the flow length was taken as the sum of the 
element sides over which the flow occurred and the length perpendicular to flow was the 
average circumference for the side. 

After all cards were prepared and assembled, a preliminary run was made to check 
for errors. The LD card (input Item 2) options (or words) 1, 2, and 5 were used. Limits 
of 1 minute and 100 iterations were set. The results of the final run are presented in the 
sections of printout starting with the element center temperatures. These temperatures 
and the varying temperature boundary temperatures are also plotted in figure 14. Im- 
proved results could be obtained by using fluid properties based on the new film tempera- 
tures calculated from the results. This applies especially to the varying temperature 
flows where the properties were based on the initial flow temperatures (1300° F). 
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Computer Output for Sample Problem 


sample prooifm* mainshaft sfal ^osepiecf 


l 0 = 

1. I 2. 

1 

3, 

l A, -0 

5, 1 6,-0 

7, 

-0 8, -0 9, 

-0 10, -0 

N T = 

1. 1 2, 

-0 

3, 

0 

t 

sf 

c 

1 

5. -0 




R0= 

1. 1.5C00C0E 

CO 

Z\ 

, 7 . OOOOOOF 01 

3, 1. OOOOOOF 

02 

A*-0 . 

5,-0. 

CD = 

l« I.IPOOOOE 

02 

? i 

, 7.200000E 01 

3, l.OOOOOOE 

03 

A, 5.0000O0F-02 

5,-0. 


6. 1 . 300000E C3 7, 1.300090F 03 8,-0. 9, 2.000000F 01 10,-0. 


HAT A FOR ALPHA. 

1 1.5000000E C3 2 l.OOOOOOOF 00 3 l .OOOOOOOF 00 A l.OOOOOOOE 00 5 l.OOOOOOOE 00 6 1 . 35000 OOE-06 

7 1 .50C00CCE-09 


DATA FOR ALPHA. 

1 l.BOOOOOOE 03 2 2 . OOOOOOOE 00 3 l.OOOOOOOE 00 A l.OOOOOOOE 00 5 l.OOOOOOOE 00 6 2 . 9AA 7000 E-06 

.7 2 • 5 A OOOOOE-C 9 


OATA FOR ALPHA. 

I l.AOOOOOOE 03 2 3. OOOOOOOE 00 3 2. OOOOOOOE 00 A l.OOOOOOOE 00 5 l.OOOOOOOE 00 6 7. 0200000E-06 

7 B.607 1A29F— l 0 8 3 . 57 1 A290F- 1 3 


CrNOUCTinN FLEMFNT TABIF (TOl) 


1 

i, l.roooF oo 

2, 3. 06 OOF 00 

3, 6 . 0000F-02 

A, 5.0000F-02 

5, 2. OOOOE 

00 

6, 7. 00 OOF 

00 

7, 

0. 



a, o . 

9, l.OOOOF 00 

10 .-0. 

11,-0. 

12,-0. 


13,-0. 





2 

1, 2 . CO OOF 00 

2, 3.0600E 00 

3, 6.0000E-02 

A, 5.00006-02 

5, 3. OOOOE 

00 

6, 8.00006 

00 

7, 

l.OOOOE 

00 


8.-0. 

9, l.OOOOF 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 





3 

1, 3. OOOOE 00 

7 , 3.O6O0E 00 

3, 6.0000F-02 

A, 7.5000F-0? 

5, A. OOOOE 

00 

6, 9. OOOOE 

00 

7, 

2. OOOOE 

00 


8,-0. 

9, 1.00006 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 





A 

1, A.COOCE 00 

2, 3. 06 OOF 00 

3, 6.0000F-0? 

A, 7.5000E-02 

5, 5. OOOOE 

00 

6, l.OOOOE 

01 

7, 

3. OOOOE 

00 


8,-0. 

9, l.OOOOE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 





5 

1, S.COOOE 00 

2 , 3. 06 OOF 00 

3, 6.0OO0F-O2 

A, 7. 50006-02 

5, 6. OOOOF 

00 

6, 1.1000E 

01 

7, 

A. OOOOE 

00 


8,-0. 

9, l.OOOOF 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 





fc 

1, 6.C000F 00 

2. 3.0600E 00 

3, 6.0000F-0? 

A, 7 .50006-0? 

5, 0. 


6, 1.2000E 

01 

7, 

5. OOOOE 

00 


8,-0. 

9, l.OOOCE 00 

10,-0. 

11,-0. 

12, -n. 


13,-0. 





7 

1. 7.C00CE 00 

2, 3.1250E 00 

3, 7.0O00F-0? 

A, 5.0000F-02 

5, 8. OOOOE 

00 

6, l • 30 OOF 

01 

7, 

0. 



a. l . o co of oo 

9, l.OOOOE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 





a 

i, e.coooE oo 

2, 3.1250E 00 

3, 7.0000F-02 

A, 9.0000E-02 

5, 9. OOOOE 

00 

6, 1. AOOOE 

01 

7, 

7. OOOOE 

00 


a. 2.CC00E 00 

9. l.OOOOE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 





9 

1. 9.00 OOF 00 

2, 3.1250E 00 

3, 7 . 00 00 E— 02 

A, 7 . 5000E— 02 

5, l.OOOOF 

01 

6, 1.5000E 

01 

7, 

8. OOOOE 

00 


a. 3.CCOOF 00 

9. l.OOOOE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 





10 

1, l.COOOF 01 

2, 3.1250E 00 

3, 7 . 00 OOF— 07 

A, 7.5000E-02 

5, 1.1000E 

01 

6, 1.6000E 

01 

7, 

9. OOOOE 

00 


8, A.OCOOF 00 

9. l.OOOOE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 





1 1 

1, 1.1000E 01 

2, 3.1250E 00 

3, 7 • OOOOE— 02 

A, 7.5000F-02 

5, 1.2000E 

01 

6, 1.7000E 

01 

7, 

l.OOOOE 

01 


a, 5.CC00E 00 

9, l.OOOOE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 





12 

1, 1.2000F 01 

2, 3.1250E 00 

3, 7 .00006-07 

A, 7.5000F-02 

5, 0. 


6, t . 8000E 

01 

7, 

1.1000E 

01 


8. 6.CC00F 00 

9. l.OOOOE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 





13 

1, 1.3000E 01 

2, 3.1950E 00 

3, 7. 00 OOF— 02 

A, 5.0000F-02 

5, 1. AOOOE 

01 

6, 0. 


7, 

0. 



8, 7.CC00F 00 

9, l.OOOOE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 
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7, 7.6000E 

01 


8, 7.2C00F 01 

9, 2.0000E 00 

10,-0. 

1 1 ,-o. 

12, 0. 


13,-0. 




78 

i, 7.enooE oi 

3.4500E 00 

3, 4 .0000 F— 02 

4, l.OOOOF-Ol 

5, 7. 9000 E 

01 

6, 0. 


7, 0. 



8, 7.3C00E 01 

9* 2* 000 OE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 




79 

1, 7 * 9000E 01 

2, 3.4500E 00 

3, 4 ■ 0000E-02 

4, 1. 0000 E- 01 

5, 8.00 OOF 

01 

6, 0. 


7, 7.8000F 

01 


8, 7.4C00E 01 

9* 2.0000E 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 




80 

l, 8 • 00 OOF 01 

2, 3.45 OOF 00 

3, 4.0000F-02 

4, 8.0000 E-02 

5, 0. 


6, 0. 


7, 7.9000E 

01 


8* 7.5C00E 01 

9. 2.0000E 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 




8 1 

1. 8.1000E 01 

2, 3.5500E 00 

3, 8.0000F-02 

4, 7 .0000 E— 0 2 

5, 8 . 20 OOF 

01 

6, 8.3000E 

01 

7, 0. 



8* 7.6000E 01 

9, 2.0000E 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 




8? 

l, 8 « 2000E 01 

2. 3.5500E 00 

3, 8 . 0000 E-02 

4, 6.5000E-0? 

5, 0. 


6, 8.4000E 

01 

7, 8.1000E 

01 


8, 7.7C00E 01 

5, 2. 000 OE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 




87 

1. 8.3000F 01 

2, 3.6350E 00 

3, 9.0000F-02 

4, 7.0000E-02 

5, 8.4000E 

01 

6, 8.5000E 

01 

7, 0. 



8, 8.1C00F 01 

9* 2.0000E 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 




84 

1, 8.4000E 01 

2. 3.6350E 00 

3, 9.0000F— 02 

4, 6.5000 E-0 2 

5, 0. 


6, 8.6000E 

01 

7, 8.3000E 

01 


8, 8.2000E 01 

5* 2.0C00E 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 




85 

l, 9.5000E 01 

2, 3.7250F 00 

3, 9. OOOOE-02 

4, 7 .0000 E-02 

5, 8.6000E 

01 

6, 8.7000E 

01 

7, 0. 



8, 8.3C00F 01 

9, 2.0000E 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 




86 

l. 8.6000F 01 

2 » 3.7250E 00 

3, 9.0000 F— 02 

4, 6.5000 E-02 

5, 0. 


6, 8.8000E 

01 

7, 8.5000E 

01 


8. 8. 4 CO OF 01 

9, 2. 000 OE 00 

10,-0. 

11 ,-0. 

12,-0. 


13,-0. 




87 

l, 8 * 7000E 01 

2, 3.8150F 00 

3, 9.0000F-02 

4, 7. 0000E-02 

5, 8.8000E 

01 

6, 0. 


7, 0. 



8, 8.5COOF 01 

9, 2.C000E 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 




88 

1, 8.EOOOE 01 

2. 3.8150E 00 

3. 9 . 00 OOF— 02 

4, 6 .5000 E-0 2 

5, 0. 


6, 0. 


7, 8.7000E 

01 


8. 8.6CC0F 01 

9, 2.000CE 00 

10,-0. 

11,-0. 

12,-0. 


13,-0. 




89 

1, 8.S000E 01 

2 1 3.24 OOF 00 

3, l.OOOOF-Ol 

4, 1 . 2500 F-0 1 

5, 9 .0000 F 

01 

6, 5.90 00E 

01 

7, 0. 



8, 0. 

9, 3. 000 OE 00 

10,-0. 

11, 2.7900E 02 

12,-0. 


13,-0. 




90 

1. 9.C000F 01 

2 , 3.2400E 00 

3, l.OOOOF-Ol 

4, 1 .0000 E-0 l 

5, 0. 


6, 6. 00 OOF 

01 

7, 8.9000E i 

01 


8. 0. 

9* 3. 000 OE 00 

10,-0. 

11, 5 ■ 4438 E 05 

12,-0. 


13,-0. 




9 1 

1, 9.1000F 01 

2, 2 • 96 5 OF 00 

3, 1.3000F-P1 

4, 8 .0000 E— 02 

5, 9.2000E 

01 

6, 9.3000E 

01 

7, 0. 



8* 0, 

5* 2*000 CE 00 

0 

1 

o 

11, o. 

12,-0. 


13,-0. 






I 


9? 

It 9.2000E 01 
8 * 0 . 

2 « 2.9650F CO 
9, 2.0000F 00 

3, 1.3000E-01 
1 0 ,- 0 . 

4, 8.0000E-02 
1 1 ,- 0 . 

5, 0. 
12 ,- 0 . 


6 , 9.4000E 
L 3 ,- 0 . 

01 

7, 

9.1000E 

01 

97 

It 9.70 00E 01 
8 , 5.1CP0F 01 

2, 3. 06 OOE 00 
9, 2.0000E 00 

3 , 6 . 000 nF -02 

10 . - 0 . 

4, 8 .00 00F-02 
11 ,- 0 . 

5, 9.40 00 E 
12 ,- 0 . 

01 

6 , 9.5000E 
13,-0. 

01 

7, 

0 . 


04 

It 9.400OF 01 
8 . 9.2CO0E 01 

2, 3.0600E 00 
9, 2.0000E 00 

3, 6.0000E-O2 
10 , - 0 . 

4, 8 .0000 F— 02 
11 ,- 0 . 

5. 0. 
12 .- 0 . 


6 , 9.6000E 
13,-0. 

01 

7, 

9.3000E 

01 

9 * 

t. 9.5000F 01 
A. 9.3C00F 01 

2, 3.1250E 00 
9, 2. OOOOE 00 

3, 7.0O00E-02 
10 ,- 0 . 

4, B .00 00 E— 0 2 
11 ,- 0 . 

5, 9. 60 00 E 
12 ,- 0 . 

01 

6 , 9.7000E 
13,-0. 

01 

7, 

0 . 


96 

1, 9 • 60 OOF 01 

a, 9.4C00F 01 

?, 3.12 50E CO 
9, 2.0000F 00 

3, 7.0000E-02 
10 ,- 0 . 

4, 8 .OOOOE— 02 
11 ,- 0 . 

5, 0. 
12 ,- 0 . 


6 , 9.8000E 
13,-0. 

01 

7, 

9.5000E 

01 

97 

It 9.7000H 01 

a, 9 .srooF oi 

2 « 3.19S0F 00 
9* 2 • 000 OF 00 

3, 7.0000F-0? 
10 ,- 0 . 

4, 8.0000E-02 
11 ,- 0 . 

5* 9.80 OOF 
12 ,- 0 . 

01 

6 , 9.9000E 
13,-0. 

01 

7, 

0 . 


98 

It 9.POOOE 01 
8 . 9-6CQ0F 01 

2, 3.1950E 00 
9, 2. OOOOE 00 

3, 7 • 0000 F— 02 
10 ,- 0 . 

4, 8 .0000 F— 02 
1 1 . -0 . 

5, 0. 
12 ,- 0 . 


6 , 1. OOOOE 
13,-0. 

02 

7, 

9.7000E 

01 

99 

It 9 . 9000F 01 
A, 9.7C00F 01 

2 , 3.2550E 00 
9, 2. 000 OE 00 

3, 5.0000F-0? 
10 ,- 0 . 

4, 8.0000E-02 
11 ,- 0 . 

5, 1.00 OOF 

12 ,- 0 . 

02 

6 , I.OIOOE 
13,-0. 

02 

7, 

0 . 


100 

l. 1 . CO OOF 02 
8 , S.8C00F 01 

?. 3.2550E 00 
9; 2.0000F 00 

3, 5 . OOOOE— 02 
10 ,- 0 . 

4, 8.0000E-02 
11 ,- 0 . 

5, 0. 
12 ,- 0 . 


6 , 1.0200E 
13,-0. 

02 

7, 

9 . 90 OOE 

01 

101 

1. 1 . Cl OOE 02 

8 , 9.9C00E 01 

2, 3.3050E 00 
9, 2.0000F 00 

3, 5.0000F-0? 
10 ,- 0 . 

4. 8.0000E-02 
11 ,- 0 . 

5, 1.02 00 F 

12 ,- 0 . 

02 

6 , 1.03 OOF 

13,-0. 

02 

7, 

0 . 


10 ? 

It 1.C200F 02 
8 , l.OCOOF 02 

2, 3.3050F 00 
9, 2. 000 OE 00 

3, 5.0000F-02 
10 ,- 0 . 

4, 8.0000E-02 
11 ,- 0 . 

5, 0. 
12 ,- 0 . 


6 , 1.0400E 
1 3,-0. 

0 ? 

7, 

I.OIOOE 

02 

101 

It 1.C30CF 02 
8 , l.OIOOF 02 

2, 3.3550E 00 
9, 2.000 OE 00 

3, 5 • OOOOE— 02 

10 ,- 0 . 

4, 8.0000F-02 
11 ,- 0 . 

5, 1.0400E 

12 ,- 0 . 

0 ? 

6 , 1.0500E 
13,-0. 

02 

7, 

0 . 


104 

t. 1.C4006 02 
8 , 1 . 0200 E 0 ? 

2, 3.35 50E 00 
9, 2. 00 OOF 00 

3, 5.0000F-02 
10 ,- 0 . 

4, 8 .000OE-O2 

11 ,- 0 . 

5, 0. 
12 ,- 0 . 


6 , 1.0600E 

13,-0. 

02 

7, 

1.0300F 

02 

10 5 

1. 1.C500F 02 

A. l.OCOOF 02 

2, 3.4150E 00 
9, 2.0000F 00 

3, 7.0000F-0? 
1 0 ,- 0 . 

4, 8 .0000 E— 0 2 
LI ,-0 . 

5, 1 . 06 OOF 

12 ,- 0 . 

0 ? 

6, 1.0700F 

13,-0. 

02 

7, 

0 . 


106 

1, 1.C600F 0? 

A, 1.C400E 0? 

2* 3.4150F 00 
9, 2. 000 OE 00 

3, 7 ■ OOOOF-O? 
10,-0. 

4, 8.0000E-02 
11,-0. 

5, 0. 
12,-0. 


6, 1.0800E 

13,-0. 

02 

7, 

1 . 05 OOE 

02 

107 

It 1 . C700E 02 
R. 1.0500E 02 

2* 3.4800E 00 
9, 2. 000 OF 00 

3, 6.0000E-02 
10,-0. 

4, 8 .OOOOE— 0 2 
1 l , - 0 • 

5, 1 .0800 E 

12,-0. 

02 

6, 1.0900E 

13,-0. 

02 

7, 

0. 


10 8 

It l.CROOE 0? 
A. l.CCOOF 02 

2, 3 . 4 8 OOF 00 
9, 7.0000E 00 

3, 6.0000E-0? 

10 ,- 0 . 

4, 8 .OOOOE— 02 

11 t -0 . 

5, 0. 
12,-0. 


6, 1.10 OOE 

13,-0. 

02 

7, 

1 -0700E 

02 

109 

l, 1.C900F 02 
fl. 1.C700F 02 

2. 3.5500E 00 
9* 2. 000 OE 00 

3, 8.0000E-02 
10,-0. 

4, 8.0000E-02 

11,-0. 

5, 1.10 OOF 

12,-0. 

02 

6, 0. 
13,-0. 


7, 

0. 


1 10 

1 , 1 . 1000F 02 

8, l.CfOOE 02 

2, 3.55 OOF 00 
9, 2. 000 OE 00 

3, 8.0000E— 02 
10,-0. 

4, 0.OOOOE-O2 

11,-0. 

5, 0. 
12,-0. 


6, 0. 
13,-0. 


7, 

1.0900E 

02 


ELFMFNT STARTING TEMPERATURES 


1 

1097.0 

2 

1 110.0 

3 

1 123.0 

4 

1163.0 

5 

1 205.0 

6 

1251.0 

7 

1042.0 

8 

1157.0 

9 

1080.0 

10 

1127.0 

1 1 

11 72.0 

l? 

1 223.0 

13 

1 016.0 

14 

1130.0 

1 5 

1052.0 

16 

1100.0 

17 

1141 .0 

18 

1174.0 

19 

1040.0 

20 

1053.0 

21 

1090.0 

? 2 

1 129.0 

23 

1 158.0 

24 

94 3.0 

25 

947.0 

26 

1054.0 

27 

1085 .0 

28 

1120.0 

29 

1146.0 

30 

941.0 

31 

946.0 

32 

1 T53.0 

33 

1085.0 

34 

1113.0 

35 

1138.0 

36 

1058.0 

37 

1065.0 

38 

1086. 0 

39 

1108.0 

40 

1132.0 

41 

106 2.0 

4? 

1 C73.0 

43 

l C 86. 0 

44 

1101.0 

45 

1 120.0 

46 

1005.0 

47 

1007.0 

4 8 

1 CO 9 .0 

49 

935.0 

50 

987.0 

51 

9 89.0 

52 

l 181.0 

53 

1 179.0 

54 

1 L 80. 0 

55 

1177.0 

56 

1177.0 

57 

1179.0 

58 

1 L 75.0 

59 

l 169.0 

60 

1165.0 

61 

1163.0 

6? 

1 162.0 

63 

1 174.0 

64 

1 175.0 

65 

l 176.0 

66 

1174.0 

67 

1169.0 

68 

1165.0 

69 

t 163.0 

70 

1163.0 

71 

1172.0 

7? 

1 173.0 

73 

1 164.0 

74 

1163.0 

75 

1163.0 

76 

1 170.0 

77 

117 1.0 

78 

1164.0 

79 

1163.0 

80 

1162.0 

8 1 

1169.0 

82 

1 169.0 

83 

1167.0 

84 

1167.0 

85 

1166.0 

86 

1166.0 

87 

1165.0 

88 

1165.0 

89 

1107.0 

50 

1060.0 

91 

8 71.0 

9? 

F81 .0 

93 

869. 0 

94 

877.0 

95 

866.0 

96 

872.0 

97 

864.0 

98 

868.0 

99 

861.0 

100 

865.0 

101 

R 60. 0 

102 

E64.0 

103 

858. 0 

104 

860. 0 

105 

356.0 

106 

858.0 

107 

855.0 

108 

856.0 

109 

853.0 

1 10 

854.0 


85 



Ft FMFNT* l. TYPF*7. SUB T YPF*4. NFXT=1. ITEMS READ- 9. ARFACOOF=4. SPECIAL A R F A *- 0 • SO. FEET 

rrNr.FNTPir. cylincfrs 

T-B=nrO. K-Fs 0 . 1 0 1 8F- 0 4 C-P= 0. 2679F 00 Mlj* 0.26B3F-04 NU= 0.5593F-04 R-0* 0.2525E 00 R-l = 0.24I7E 00 OMEGA* 1570. 

V- AX= G.3134F 00 df= 0.1502E 06 pP**I/3= 0.8903F 00 H COFFF* 0.9706356F 00 COEFF+AREA* 0 . 64 1 632 2E-02 

FI FMFNT= 1 . TYPF = 6. SUB T YP F = 2. NFXT = 0. ITEMS RF AD=t 4. AREA C0CF=3. SPECIAL A R F A =— C • SQ. FEET 

V AR.T FMP .RACI AL SFAI GAP 

T-B*1 0 62 . K— F = 0.°e c 4F-C5 C-P* 0.2647E 00 MU* 0.2562E-04 NU = 0.4990F-04 R-AV= 0.2550E 00 GAP* 0.4333E-04 OMEGA* 1570. 
R-0* 0.2992F 00 FI AG-C=0. FLAG-I*?. ANN. NBR. = 94. FLO.NRR •= 0. MASS FLOW* 0.2382E-02 

RF* 0.2046F C7 PR**l/3= 0.8RSCF 00 RE-C* 0.6754F 07 H COEFF= 0.3961579E 00 COEFF*ARFA= 0. 3173 645E- 02 

HT .CFN .FACTOR* 0.10CP0F 01 

FI FMFNT* ?. TYPF=7. SUB TYPF=4. NFXT*0. ITEMS REAn= q. AR F A C0DE=4. SPECIAL ARFA=~0 » SO.FFFT 

CONCENTRIC CYLINCFRS 

T— R = J 3 00 . K-F* p .1 0? IE— 04 C-P= 0.2680F 00 MIJ* 0.2690F-04 NU= 0.5632F-04 R-0 = 0.2525E 00 R-I* 0.2417F 00 OMEGA* 1570. 

V - AX = C.3134F 00 PF* 0.1492F 06 PR**l/3= 0.8904F 00 H COFFF* 0.9672406F 00 CO€FF*ARFA= 0.63 938806-02 

FI FMFNT* 3. T YPF =2 . SUBTYPF = 4. NFXT=0. TTFMS P E AD* 9. AREA C0DF=4. SPECIAL AREA*- 0 . SO.FFFT 

CONCENTRIC CYL INTERS 

T-R = 130O. K— F = 0.1024F-04 C-P= 0.2682F 00 MU* 0.7696F-04 NU* 0.5664E-04 R-0* 0.2525F 00 R-I= 0.2417F 00 OMEGA* 1570. 

V-AX* T.3134F 00 RF= 0.1483E 06 PR**l/3= 0.8905E 00 H COEFF* 0.9643760E 00 COEFF+ARE A= 0 . 95 6 24 1 6F-0 2 

FLFMFNT* 4. TYPE*?. SUB7YPF=4. NFXT=0. ITEMS PFA0= o. ARF A C0DE=4. SPFCIAL ARFA=— 0. SO.FFFT 

ttncf ntp i r. cyl inters 

T— R = 1 3 TO . K-F= 0.1 0 33F- 04 C-P= 0.2687E 00 MIJ= 0.2716E-04 NU= 0.5774F-04 R-0* 0.2525E 00 R-t = 0.2417E 00 OMEGA* 1570. 

V- AX* T.3134F 00 RF= 0.I455F 06 P(?**l/3= 0.8907E 00 H COEFF* 0.9549769E 00 COEFF*AR£ A= 0. 946921 8E-02 

FI FMFNT* 5. TYPF = 2. SURTYPF=4. NEXT=0. ITFMS P E AO* 9. AREA C0DE = 4. SPECIAL A R F A =- 0 . SO. FEET 

C0NCFNTPIC CYl INTERS 

T— R= 1 3 CO . K-F* 0.1042E-04 C-P* 0.2692F 00 MU= 0.7736F-04 NU* 0.5890E-04 R-0= 0.2525F 00 R-I= 0.2417E 00 OMEGA* 1570. 

V-AX* C.3134F 00 RF= 0. 1427F 06 PR**l/3= 0.8910F 00 H COEFF* 0.9452609F 00 COEFF*ARFA = 0. 93 7 287 7E-02 

FI FMFNT* 6. T Y PF =7 . SUPTYPF=4. NEXT=l. ITEMS READ* 9. AREA COOF= 1. SPECIAL ARE A=-0 • SO. FEET 

CONCENTRIC CYL INCURS 

T-R=1300. K— F = C . I 05 IE - 04 C-P= 0.2698F 00 MIJ= 0.2759E-04 NU= 0.6018F-04 R-0* 0.2550E 00 P-I* 0.2417E 00 OMEGA* 1570. 

V-AX* C. RF= 0.1727F 06 PR**l/3* 0.8912E 00 H COEFF* 0.9318566E 00 COEFF*ARFA= 0 . 74 6 5 160F-02 

Ft FMPNT* 6. T YPF *2 . SLR T YPF =4 • NFXT=0. ITFMS READ* 9. AREA C0DE = 4. SPFCIAL ARFA=-0. SQ.FEFT 

CONTENTS I C CYLINCFRS 

T-R=t30n. K-F* 0 .1 05 IE- 04 C- P = 0.2698E 00 MU* 0.7759E-04 NU= 0.6018E-04 R-0* 0.2525E 00 R-I* 0.24I7E 00 OMEGA* 1570. 

V-AX* 0.3134F 00 PF* 0.1396F 06 PR++1/3* O.0912F 00 H COEFF* 0.9348986F 00 COEFF*ARFA= 0. 92 701 2 9E -02 

FLFMFNT* 12. TYPE*?. SURTYPF=4. NFXT=l. ITFMS READ* 9. AREA CODE* 1 . SPECIAL A REA=-0. SO. FEET 

CrNCFN'TR IC CYL INTERS 

T-R=1300. K-F* 0 .1 045F-04 C-P* 0.2694F 00 MU= 0.2745E-04 NU* 0.5940E-04 R-0= 0.2608E OO R-I* 0.2417E 00 OMEGA* 1570. 

V-AX* 0. PF* 0.2543E 06 PR**l/3= 0.8911F 00 H COEFF* 0.9354691F 00 COEFF+AREA* 0. 892 8836E-02 

FLFMFNT* 12. TYPP=6. SUBTYPF=1. NFXT=0. ITFMS RpAD=12. AREA C0DE=5. SPECIAL AREA* 0.23562E-02 SO. FEET 

VAR. TEMP. JUCL FLOW, LIwUlO 

T— B = 12 85 . K-F* 0 . 1 042E-04 Z~P= 0.2b93t 00 MU= 3.273?. -34 MASS FLOW* 0.56856 01 D-HYDR = 0. 3333E-0Z L-=LOrf* 3.3125E-32 
ANN.NBR.* j. FLOW NbR* 0. ELAG-C=0. FLAG~I=1. X-SfcC= 0.3142E-03 

RE* 0.6920E 01 PI* 0.7074E 00 U CuEFF* 3 . 3122549F-01 COEFF^ARF A= 0 .7357350F- 04 

El FMFNT* 57. TYPE*?. SUP TYP F=4 . NFXT=1. ITEMS R EAO= 9. AREA C0DE=1. SPECIAL A R F A *— 0 • SQ.FEET 

CONCENTRIC CY1 INTERS 

T-R*I3m. K-F* T.1 036F-C4 C-P* 0.2689F 00 MU= 0.2725F-04 NU= 0.5824E-04 R-0* 0.26F2F 00 R-I* 0.2417E 00 OMEGA* 1570. 

V-AX* C.4726F 00 RF = 0.3366F 06 PR**l/3= 0.890BE 00 H COEFF* 0.9456614E 00 COEF F *ARE A* 0. 92 2 83 04E-02 

FI FMFfrT* 52. T Y PF = 7 . SUBTYPF = 4. NFXT=0. ITFMS READ* 9. AREA C00F=4. SPECIAL ARE A=-0. SQ.FEET 

CONT.PNTRTC CYl IN^prs 

T— R* 1 3 CO . K-F* 0*103 6F -04 C-P* 0.2689F 00 MU* 0.2725F-04 NU* 0.5824F-04 R-0= 0.2633F 00 R-I* 0.2417E 00 OMEGA* 1570. 

V-AX* C.4726T 00 RE* 0.2949F 06 PR**t/3= 0.8908F 00 H COEFF* 0.9450770F 00 COEF F*AR EA* 0 . 1 1 0761 9E-0 1 

FLFMFNT* 54. TYPF=2. SURTYPF=4. NFXT=l. ITEMS READ* 9. AREA C0DE=1. SPECIAL AREA=-0. SQ.FEET 

CONCENTRIC CYL INOFRS 

T— R = 1 3 TO. K-F* 0.1036F-04 C-P* 0.2689E 00 MU* 0.2725E-04 NU* 0.5824E-04 R-0* 0.2712F 00 R-I= 0.2417E 00 OMEGA* 1570. 

V-AX* 0.4776F 00 PF* 0.4090E 06 PR**l/3= O.8908E 00 H COEFF* 0.9476610E 00 COHFF*ARFA= 0. 672963 3E-02 

FI FMFNT* 54. TYPF*7. SURTYPE=4. NFXT=0. ITEMS READ* 9. AREA C0DE=5. SPECIAL AREA* 0.I7617E-02 SQ.FEET 

CONCENTRIC CYLINCFRS 

T— B = 1 3 00 • K— F = 0.10366-04 C-P* 0.2689F 00 MU* 0.2725F-04 NU* 0.5824E-04 R-0* 0.2692F 00 R-I* 0.2417E 00 OMEGA* 1570. 

V-AX* 0.4726E 00 RF* 0.37R7E 06 PR**l/3= 0.8908E 00 H COEFF* 0.9467012E 00 COEF F*AR E A* 0 . 1 6 678 04E-02 

FLFMFNT* 57. TYPF = 7. SURTYPF*4. NFXT=0. ITFMS R F AD* 9. AR F A C0DE = 4. SPECIAL A R E A =— 0 . SO.FFFT 

CONCENTRIC CYL INTERS 

T-R* 1 3C0 . K-F* 0 .1 036F-04 C-P= 0. 2689E 00 MU* 0.2724F-04 NU* 0.5818E-04 R-0* 0.2742E 00 R-I* 0.24L7F 00 OMEGA* 1570. 

V-AX* 0.4776F 00 RF* 0.4524F 06 PR**l/3= 0.8908E 00 H COEFF* 0.9497571F 00 COEFF+ARFA* 0 . 1363408E-01 

FLFMFNT* 58. TYPF=7. SURTYPF=4. NEXT=0. ITEMS RFAD* 9. ARF A C0DE=4. SPECIAL ARFA=-0. SQ.FEET 

CONCENTRIC CYL INTERS 

T— 8=13 CO. K-F* 0.1035E-04 C~P= 0. 2688E 00 MU* 0.2722F-04 NU* 0.5807E-04 R-D* 0.2742F 00 R-I* 0.2283E 00 OMEGA* 1570. 

V-AX* 0.937 1 F 00 RF* 0.6228E 06 PR**1 /3= 0.8908E 00 H COFFF* 0.9274974E 00 COFFF+AREA* 0. I 464599E -01 
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EL FMF NT = 30. TYPF=3. SUBTYPE*!. NFX7*1. ITEMS RF AO* A. AREA C0DE*2. SPECIAL AREA*-0. SQ.FEET 

L lOUrC FTLM COOL ING 

T-B* S?5. W/L = C.I000F-01 C-P* O. 5030E 00 L-FLOW* 0.88A9E 00 H COFFF= 0 . 5 6 0 A 258E-02 COEFF* ARE A* 0. Al 91 5 8 LE-OA 

FLFMFNT* 30. TYPE = 6 . SUBTYPE = 2. NEXT*0. ITFMS PEA0*1A. AREA fOOE=3. SPECIAL AR FA=— 0 • SO.FFFT 

VAR. TEMP. RADIAL SFAL GAP 

T-B* 955. K-F= C.90A8E-05 C-P= 0. 2611F 00 MU* 0.2A23E-0A NJ = 0.1571E-03 R- AV* 0.2796E 00 GAP* 0.A333E-0A OMEGA* 1570. 
R-D* 0.2997E 00 FLftG-C=0. FLAG- 1=2. ANN. NBR. =1 OA . FLP.NflR.* 2A . HASS FLOW* 0.595AF-02 

RF= 0.781 IE 06 PR**l/3= C.R876E 00 RE-C* 0.67SAE 07 H COEFF= 0.3706976E 00 CO£FF*ARFA= 0. 271 32 75E- 02 

HT .GEN .FACTflR* 0.10000E 01 

FI EMFNT = 36. TYPF*?. SU8TYPE = 3. NEX7=0. ITEMS RF A0= A. AREA C00E=3. SPECIAL AR E A=-0 • SO.FFET 

l TOUT C FTI M COOL ING 

T-B* 575. W/L = C.lOOOF-01 C-P* 0. 5030E 00 L-FLOW* 0.2339E 00 H COEFF= 0 .2 15CA92E-0 1 COEFF* ARE A* 0. 22A3075E-03 

FLFMFNT* AT. TYPE*!. SUBTYPE*!. NFXT=l. ITEMS READ* A. AREA C0DE=2. SPECIAL AREA=-0. SO. FEET 


L I01JID FILM Com ING 

T-B* 525. W/L = r.lOOOE-Ol C-P= 0. 5030F 00 L-FLOW* 0.91B9E 00 H COEFF* 0.5A73936E-0? COEFF*ARE A= 0. A1 91 739E-0A 

FL FMF NT * Al. TYPE*!. SURTYPE = 3. NEXT=0. ITEMS RFAO= A. AREA C0DE=3. SPFCIAl AREA=-0. SO. FEET 

LIOUIC FILM COOLING 

T-B* 525. W/L = C.IOOOE-Ol C-P= 0. 5030E 00 L-FLOW* 0.2339E 00 H COEFF* 0 .2 1 50A92E-01 COEFF*AREA= 0. 1959231E-03 

FLFMFNT* A6. TYPF=3. SUBTYPF*!. NEXT=0. ITFMS READ* A. AREA CnOE=3. SPECIAL AR FA=-0. SO. FEET 


L IOUTT FILM CHCL ING 

T-R= 525. W/L = C.7000E-01 C-P* 0. 5030F 00 L-FLOW* 0.2A50E 00 H COEFF* 0 . A 106 122E-0 1 COEFF*AREA= 0. 50882AOE-03 

FLFMFNT* A9. TYPF=3. SUBTYPF=3. NEXT=l. ITEMS RFAP* A. AREA C0DE=2. SPFCIAL ARFA*— 0 . SO. FEET 


L IOUID FILM COOL ING 

T-B= 525. W/L* C .1 500E-01 C-P* 0.5030E 00 L-FLCW* 0.9608E 00 H COEFF* 0 . 7852831 E-02 COEFF*AREA= 0. 9A31 286E-0A 

FL FMF NT = A9. TYPF*3. SUBTYPF=3. NFXT=0. ITFMS REAP* A. AREA C0DF = 3. SPECIAL ARFA=-0. SQ.FEET 

L IODIC FILM COOL ING 

T-B* 525. W/L* C.2000E-01 C-P* 0. 5030E 00 L-FLOW* 0.2A50F OP H COFFF* 0 .A 106 122E-01 COEFF*AREA* 0. 5202905E-03 

FLFMFNT* 50. TYPE = 3. SUBTYPF = 3. NEXT = 0. ITFMS RF AO* A. AREA C0DF=2. SPECIAL AREA*— 0 • SO. FEET 

l I QUID FILM COOL ING 

T-B* 575. W/L* C.1500F-01 C-P* 0.5030F 00 L-FLOW* 0.9608F 00 H COEFF* 0 .7 85 2 831 E-02 COFFF+ARE A* 0. 9A31 2 86E-0A 

FLFMFNT* 51. TYPF=3. SUBTYPF = 3. NFXT = 1. ITFMS RE AO* A. ARFA CODF=l. SPECIAL ARE A=- 0 . SO.FFET 

L IOUTT FTLM COOL TNG 

T-B* 525. W/L* C.5000F-02 C-P* 0. 5030E 00 L-FLOW* 0.2A50E 00 H COEFF* 0 . 10265 31 F-0 1 COEFF*AREA= 0. 1 300 726E-03 

FLFMFNT* 51. TYPE*!. SUBTYPF=3. NFXT=0. ITEMS READ* A. AREA C0DE=2. SPECIAL ARF A=-0 . SO.FEFT 

L IOUIC FILM COOL ING 

T-B* 525. W/L* C.1 500E-01 C-P= 0. 5030E 00 L-FIOW* 0.9608F 00 H COEFF* 0 .7 85 2831 E-0 2 COEFF *ARE A* 0. 9A31 286E-0A 

FLFMFNT* A8. TYPF=3. SUBTYPE*!. NEXT=0. ITFMS RFAP* A. AREA CODE=l. SPFCIAL AREA=-*0. SQ.FFFT 

L IOUIC F !LM COOt ING 

T-B* 575. W/L* 0 .50 OOF— 02 C-P* 0. 5030E 00 L-FLOW* 0.2A50F 00 H CnEFF* 0 . 102653 1 E-0 1 COFFF* ARF A = 0. 1 2 72 060E-03 

FLFMFNT* A5. TYPP*^. SUBTYPF = 3. NEXT = 0. ITFMS RF AD= A. ARFA CO0E=2. SPECIAL AR F A=- 0 . SO.FEFT 

1 IOUIC FILM COOL ING 

T-B* 5?5. W/L* C.7000F-01 C-P= 0. 5030F 00 L-FLOW* 0.9189E 00 H COFFF* 0 . 109A 787E-01 COEFF* ARE A* 0. 1 2 57522E— 03 

FLFMFNT* 81. TYPF = 3. SUB T YPF =3. NEXT=0. ITFMS RF AD* A. AREA C00E=3. SPECIAL ARE A=-0 . SO.FEFT 

L IOUIC FTLM COOL ING 

T-B* 525. W/L* C.1503F-0L C-P* 0. 5030F 00 L-FLOW* 0.236IE 00 H COFFF* 0 . 3 1 95680E-01 COFFF* ARE A* 0. 396003AE-03 

FI FMFNT* 83. TYPF = 3. SUBTYPE*!. NEXT=0. ITFMS RF AO* A. AREA C00E = 3. SPFCIAL AR E A =- 0 . SQ.FEET 

L IOUIC FILM COOL ING 

T-R* 525. W/L* 0.1 50 OF — 01 C-P= 0. 5030C 00 L-FLOW* 0.2361E 00 H COEFF* 0 .3 195680E-0 1 COEFF *ARE A* 0. A561 709E-0 3 

FLFMFNT* 85. TYPF=3. SUBTYPF*3. NFXT=0. ITFMS READ* A. ARFA C0DE=3. SPECIAL ARFA*— 0. SO.FFET 

LIOUIC FILM COOLING 

T-R* 525. W/L* 0.1 500F— 0 1 C-P* 0. 5030F 00 L-FLOW* 0.2361E 00 H COEFF* 0 .3 19568 0E-01 COEFF*AREA= 0. A67A653E -03 

FLFMFNT* 87. TYPE*!. SURTYPF = 3. NEXT=1. ITFMS RE AO* A. AREA C0DE=2. SPECIAL ARFA=-0. SO.FEFT 

L lOllir F II M COOL ING 

T-R* 525. W/L* 0 .1 500E-0 1 C-P= 0. 5030E 00 L-FLOW* 0.1010E 01 H COEFF* 0 .7A6660 1 E-0 2 COEFF*ARE A* 0. B 802897E-0A 

FLFMFNT* B7. TYPF=3. SUBTYPE*!. NEXT=0. ITEMS READ* A. AREA C0DE = 3. SPECIAL AREA*- 0 * SQ.FEET 

i tain r film cool ing 

T-B* 525. W/L* 0 .1 500E-0 1 C-P* 0. 5030E 00 L-FLOW* 0.2631E 00 H COEFF* 0 .2 867731 F-01 COFFF*AREA= 0. A2962 82E-03 

FLFMFNT* 80. TYPF=3. SUBTYPE*!. NFXT=1. ITEMS RE AO* A. AREA C0DF=1. SPECIAL ARE A*— 0. SQ.FEET 

t IOUIC FILM COOL ING 

T-R* 575. W/L* C.2500F-01 C-P* 0. 5030F 00 L-FLOW* 0.2669E 00 H CHFFF* 0 .A7 I 1502F-0 l COEFF*ARE A* 0. 7058523E-03 

FLFMFNT* 88. TYPF = 3. SUBTYPF = 3. NEXT=0. ITFMS RF AD* A. AREA C0DE = 2. SPFCIAL ARE A*- 0 . SO.FFFT 
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H COEFF- 0.7466601E-02 COEFF*AREA= 0. 81741 19E-04 


1 f CtJ I r film cnoi IMG 

T— B= 575. W/l = 0.1500F-01 C-P = 0. 503 OF 00 L-FLOW- 0.1010F Oi 
FI.FMFNT= 36, TYPF-3. SUPTYPF-3. NEXT-0. ITEMS READ- 4. AREA COOE-1. SPFCIAL AREA— -0. SO. FEET 

L TOUT r FILM COOL IMG 

T— R= 575. W/l = 0.7 500E- C I C-P= 0. 5030E 00 L-FLOW- 0. 2669F 00 H COEFF- 0.47115026-01 COEFF*AREA- 0.6892005E-03 

FLFMPMT= 34. TYPF-3. SUBTYPF-3. NFXT-O. ITEMS BEAD- 4. AREA CODE-l. SPECIAL ARFA— 0. SO. FEET 

i tout r film cnoi TNG 

T-R= 575. W/L = 0.2500E-01 C-P= 0. 5030E 00 L-FLOW- 0.2669F 00 H COEFF- 0.47J1502E-01 COEFF*AREA- 0. 6725486E-03 

FL FMFNT- 87. TYPF-3. SURTYPE-3. NEXT-0. ITEMS RFAO- 4. AREA CODF-1. SPECIAL AREA— 0. SO. FEET 

I lOlHT FILM COOL IMG 

T-B- 5 25. W/L = 0.2500E-01 C-P = 0. 5030F 00 L-FLOW- 0.7669E 00 H COEFF- 0 .471 1 502E-01 COEFF*AREA= 0. 58304176-03 

FL FMFNT- 77. T YPF- * . SURTYPF-3. NFXT-O. ITEMS READ- 4. AREA CODE-l. SPECIAL ARFA— 0. SO.FEFT 

L TOUT C FILM CHOI TNG 

T-R= 525. W/L = C. 25006-01 C-P- 0.5030F 00 L-FLOW- 0.2669C 00 H COEFF- 0 .471 1 502E-0 1 CO€FF*AREA= 0.4292470E-03 

Fi FMFNT- 72. TYPF-3. SUBTYPF-3. NEXT-0. ITFMS RFAO- 4. AREA CODF-1. SPECIAL ARFA— -0. SO. FEET 

l fouir FILM cnoi IMG 

T- R= 525. W/L- C.7500F-01 C-P- C. 5030F 00 l -FLOW- 0. 2669F 00 H COFFF- 0.4711502F-01 C06FF*ARE A= 0. 45 25801E-03 

FI FMFNT- 65. TYPF-3. SUPTYP6-3. NFXT-O. ITEMS RFAO- 4. AREA COOF-?. SPFf. IAL AREA — 0. SO.FFFT 

L IQUI C FTI.M COOL 1*0 

T-B- 5?5. W/L- C.120PF-01 C-P= 0. 5030F 00 L-FLOW- 0.8849E 00 H COEFF- 0.6a?ltl0E-0? COEFF*AREA= 0. 1 005979E-03 

FI FMFNT- 66. TYPF-3. SURTYPF-3. NEXT-0. ITFMS REAP- 4. AREA COOE-7. SPECIAL AREA— 0. SO. FEET 

I TOUIT FILM Cnoi IMG 

T- R= 525. W/L = O.UOOF-OI C-P = 0. 5030F 00 l-FLOW- 0.8R49E 00 H COFFF- 0 .625 2684 E-0 2 COEFF *AREA= 0. 1 01436 3E-03 

FL FMF NT = 67. TYPF=3. SUBTYPF-3. NFXT-O. ITFMS o^AO- 4. AREA CODE-2. SPECIAL AREA — 0. SQ.FFFT 

l IOUTC FILM COOL ING 

T-R = 575. W/L- G.12nOF-01 C-P- 0. 5030F 00 L-FLOW- 0.8849F 00 H COEFF= 0 .6 8 2 l 1 l OF-O 2 COEFF*AREA= 0. 1257474E-03 

FI FMFNT- 73. TYPF-3. SURTYPE-3. NEXT-0. ITEMS READ- 4. AREA CODE-3. SPECIAL AREA — 0. SQ.FEFT 

I T OH 1 r FILM COOL ING 

T— R = 525. W/L - 0.1 000F-01 C-P= 0. 5030F 00 L-FLOW- 0.2298F 30 H COEFF- 0 .2 1 88 860F-0 1 COFFF*ARFA= 0. 1454094E-03 

FL FMF NT = 73. TYPF-3. SUBTYPF-3. NEXT-l. ITFMS RF AO- 4. AREA CO0E=2. SPECIAL ARF A— 0. SO. FEET 

l IOUI T FI! M COOL ING 

T-B- 575. W/l.= Q.1 500F-01 C-P= 0. 5030F 00 L-FLCW- 0.90846 00 H COEFF- 0 . 8 305 81 2 E-02 COEFF*AREA- 0. 1257561E-03 

FL FMFNT- 78. TYPE = 3. SUBTYPE =3. NEXT-0. ITEMS READ- 4. AREA COOE-3. SPECIAL AREA — 0. SO. FEET 

LIOUIT FILM COOLING 

T-R = 525. W/L = O.lOOne-Ol C-P- 0.5030E 00 l-FLOW- 0.2298E 00 H COEFF- 0.2188860E-01 COEFF*AREA= 0. 1317997E-03 

ELFMFNT- 79. TYPE-3. SUBTYPE-3. NFXT=0. ITEMS RF AD- 4. AREA CODE-?. SPECIAL AREA— 0. SO. FEET 

I 10(1 T T FI! M COOL ING 

T-R= 525. W/f - 0.I5O0F-0I C-p= C.5030F 00 L-FLCW- 0.90R4F 00 H COEFF- 0 . 8 305 81 2F-0 2 COEFF*ARE A- 0. 1 25756 IE-0 3 

FI.FMFNT- 80. TYPF-3. SUBTYPE-3. NEXT-1. ITEMS RFAO- 4. AREA CODE-1. SPECIAL AREA— 0. SO. FEET 

i Tom c f i! m com ing 

T— R— 5 75. W/L = C . t 5 OOF- 01 C -P* 0. 5030F 00 L-FLOW- 0.2325E 00 H COFFF= 0.3245161F-01 COEFF* ARE A- 0. 1 954037E-03 

FI EMFNT- 30. TYPF-3. SUPTYPE-3. NFXT-O. ITEMS READ- 4. AREA COOF-2. SPECIAL AREA— 0. SO. FEET 

l IOUTC FT! M cnoi TNG 

T-R= 525. W/L- C.L500F-01 C-P= 0. 5030F 00 L-FLCW- 0.9084F 00 H COEFF- 0. 8305812E-02 C06FF*AREA= 0. 1006 049E-03 

FI FMF NT- 75. TYPF-3. SUBTYPE-3. NEXT-0. ITFMS RF AO- 4. AREA COOE-1. SPECIAL AREA— 0. SO. FEET 

I I QUIT FI! M COOL TNG 

T— B= 525. W/L- C .1 500F-0 1 C-P= 0. 5030E 00 L-FLCW- 0.2325F 00 H COEFF- 0.3245161E-01 COEFF*AREA- 0. 21 558 12E-03 

FI FMF NT - 70. TYPF-3. SUBTYPF-3. NCXT-O. ITEMS READ- 4. AREA COOE-1. SPECIAL AREA— 0. SQ.FEET 

i louir f ii m com ing 

T— R= 575. W/L- C.1500F-01 C-P= 0. 5030F 00 L-FLOW- 0.2325F 00 H COFFF- 0 .3245 16 1F-01 C0EFF*AREA- 0. 2375287E-03 

FL FMF NT = 62. T Y PF = 3 . SUBTYPF-3. NEXT-l. ITEMS RFAD- 4. AREA COOE-1. S°FCIAL AREA — 0. SO. FEET 

I I QUIT FTLM CHOI IMG 

T-R- 525. W/L- 0.1503F-01 C-P= 0.5030F 00 L-FLOW- 0.73756 00 H COFFF- 0 .3 245 16 IE~0 1 COEFF*ARE A- 0. 1 874742E-03 

FL FMF NT = 67. TYPF-3. SUB TYPE- 3. NEXT-0. ITEMS R E An- 4. AREA CODE-4. SPETIAL AREA — 0. SO. FEET 

i tout r f n m cnn? ing 

T-R= 575. w/L= 0.1009F-01 C-P= 0. 5030F 00 L-FLOW- 0.8613F 00 H COEFF- 0 .5840C09E-02 COEFF*ARE A= 0. 67O6822E-04 

FI FMFNT- 61. TYPE-’. SURTYPF-3. NEXT-0. ITFMS READ- 4. AREA C00E-4. SPECIAL AREA— 0. SQ.FEFT 

L IOUTC FT! M COOL ING 

T— 8= 575. W/L- C.5000E-07 C-P= 0.5030F 00 L-FLOW- 0.8613E 00 H COFFF- 0.2920005E-02 COEFF*AREA- 0. 41 91 764E-04 

FI FMFNT- 11. TYPF=6. SUBTYPE-1. NFXT-O. ITEMS READ-12. AREA COOE-5. SPECIAL AREA- 0.16115E-02 SQ.FEET 
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VAR. TEMP. OUCT FLUW, LIWU1U 

T — B= 12 55 . K-F = 0.1035E-U4 S-P* 0.2o38E 00 MU= 3.2721E-34 MASS FLOW* 0.56S5E 01 D-HYDR* 0. 3333E-02 L-=L0W* 0.733SE-02 

ANN .NBR . = 0. FL3W NEJR* 12. FLAG-C=0. FLAJ-l=l. X-SCC* 0.3142E-33 

RE* 0.6965E 03 PR* 0. 70o9E 00 H C0EFF= 0. 2280963E-0L COE FF#AR EA= 0 .3675772E-04 

a LMEV T = L 7. TYPE = 6. SUBTYPE*!. NEXT=U. ITEMS READ* 12 . AREA C3DE=5. SPECIAL AREA* 0. 161 155-02 SO. FEET 

VAK.TEMP. OUCT FLOW, LIQUID 

T— B = 1240 • K-F = 0.1028E-04 Z-P= 0.2684b 00 MU* 0.2705E-34 MASS FLOW* 0.5685E 01 D-HYDR* 0.3333E-02 L-=LOW* 3.U33E-01 

ANN.N8R.* 0. F w JW N3K= It. FLAQ-C=0. FLAG~l=l. X-S EC= 0.3142E-J3 

RE= 0.7006E Oi PR* 0. 7064E 00 H COcFF* 0 . 1982752E-01 COEFF*AREA= 0 .3195205E-04 

EL fcM£MT= 16. TYPE=6. SUBTYPE=l. NEXT = 0. HEMS READ* 12 • AREA CJDE=5. SPECIAL AREA= 0.23562E-02 S3. FEET 

VAR. TEMP. OUCT FLOW, LIQUID 

T— 8= 12 10 . K-F = 0 . 10 19E-04 C-P= U.2679E 00 MU* 3.2685E-04 MASS FLOW* 0.5685E 01 D-HYDR* 0. 3333E-02 L-FLOW* 0.1643E-01 

ANN .NBR • = 0. FLJW NdR* 17. =LA3-C*0. FLAG-I*!. X-SEC* 0.3142E-03 

RE* 0.7057E Oi PR* 0.7059E 00 H CUEF F= 0 . I 765946E-0 l COEFF*AREA= 0 .4 160923E- 04 

EL EMENT * 15. TYPE=6. SUBTYPE=1. NEXT=0. ITEMS READ* 12 . AREAC00E = 5. SPECIAL AREA 35 0.15087E-02 SQ. FEET 

VAR. TEMP. OUCT FLOW, LIQUID 

T-B=U80. K-F = 0. 1009E-04 C-P* 0.2673E 00 MU* 3.2661E-04 MASS FLOW* 0.5685E 01 D-H YDR= 0. 3333E-02 L-FLOW= D.2133E-31 

ANN. NBR. = 0. FLOW NBR= 16. FLAG-C=0. FLAG-I=l. X-S EC= 0.3L42E-03 

RE* 0.7121E 03 PR* 0. 7053E 00 H COEFF= 0 . 16 143 LO E-0 l COEFF*AREA= 0 .2435057E- 04 

a EMEN T = 20. TYP E = 6. SU8TYPE*1. NEXT=0. ITEMS REA0=12. AREA CHOP-s. sdccial arfa* 0.15087E-02 SQ.F55T 

VAR. TEMP. OUCT FLOW, LIQUID 

T— B= llo5 . k-F = 0. 1 009E-04 C-P= 0.2673E 00 MU= 3.2661E-04 MASS FLOW* 0.5685E 01 D-HYDR= 0.3333E-02 L-FLOW* 0.2497E-01 

ANN .N8R . * 0. Fl 3W N8R= 15. =LAG-C*0. FLAG-I = l. X-SEC* 0.3142E-03 

RE* 0.7121E 03 PR* 0. 7053E 00 H COEFF= 0 . 152* 155E-0 1 C0EFF4AREA* 0 .229949 2E- 04 

a EMEN T* 19. TYPE =6. SUBTYPE*!. NtXT=0. ITEMS READ= 12 • AREA C3DE=5. SPECIAL AREA= 0.15708E-02 SQ. FEET 

VAK.TEMP. OUCT FLOW, LIQUID 

T — B = 1 150. K-F = 0.1 00t>E-04 S-P* 0.2671E 00 MU* 3.2655E-04 MASS FLOW* 0.5685E 01 D-HYDR= 0.3333E-02 L-FLOW* 0.2343E-DL 

ANN .NBiR . = 0. FL 3 W NBR* 20. FLAG-C=0. FLAG-I*!. X-SEC* 0.3142E-33 

R E- 0.7137E 03 PR* 0. 7051E 00 H COEFF* 0 . 1 4 5 638 8E -0 l COE FF *AR EA= 0 .2287223E- 04 

FLFMFNT= 7. T YPF = 6 . SUBTYPF=2- NFXT=0. ITFMS RFA0=14. AREA C0DF=3. SPECIAL AR EA=-0« SQ. FEET 

VAR.TFMP .RACIAL <FAl GAP 

T- R= 976. K-F = C^PPF-O? C-P* C.2628E 00 MU- 0.2488E-04 NU= 0.4696E-04 R-AV* 0.2604E OC GAP= 0.4333F-04 OMFGA* 1570. 

R-n= 0.P992F 00 FLAG-C = 0. FLAG-I=2. ANN. NBP.= 96. FLO.NBR.* 1. MASS FLOW* 0.2382F-02 

RF* 0.2267F 07 PR**l/3= C.8883E 00 RE-C= 0.6754F 07 H COEFF* 0.3824877F 00 COEFF*AREA= O . 3650757E-02 

HT. GEN .FACTOR* 0. 10000E 01 

El FMFNT* 13. TYPF=fc. SUBTYPE*?. NFXT=0. ITFMS REA0=14. AREA COOE=5. SPECIAL AREA* 0.16805E-01 SQ.FEFT 

VAR.TFMP.RAOIAL SEAL GAP 

T-B* 946. K-F - 0.9431F-05 C- P - C.2634F 00 MU* 0.2511F-04 NU= 0.4919F-04 R-AV* 0.2675F OC GAP* 0.9152E-03 OMEGA* 1570. 

R-P= 0.299PF 00 FLAG-C=0. F LAG-I -2. ANN . NRR . = 98. FL O.NBR.= 7. MASS FLOW* 0.2382F-02 

RE= 0. 22P4F 07 PR**l/3= C.8885F 00 RE-C* 0.2278E 06 H COEFF= 0.2927B13E 00 COEFF*AR E A= 0. 492 0307E-02 

HT.C-FN .FACTOR* 0.19M5E 01 

Fl FMFNT* ?4. TYPF=6. SUBTYPE=2. NEXT=0. ITFMS REA0=1 4. ARFA COOF=5. SPECIAL ARFA= C.14366E-01 SQ.FEET 

VAR.TFMP.R ADI AL SEAL GAP 

T-B=10C5. K— F = 0.9170E-05 C-P= C.2618F 00 MU* 0.2451F-04 NU = 0.2928E-04 R-AV* 0.2744E 00 GAP* 0.1078E-02 OMFG A= 1570. 

R-0= 0.29S7F 00 FL AG-C = 0 . FLAG-I = 2. ANN- NBR. =102. Fl C-NBR.* 100 . MASS FLOW* 0.5954E-0? 

RF= 0.4037F 07 PR**1/3= 0.8879E 00 RF-C* 0.I900F 06 H COFFF* 0.4274941E 00 CnEFF ♦ARE A* 0. 6141509E-0? 

HT .GFN .FACTOR* 0.18512F 01 

FL FMFNT = 9?. TYPF=5. SUBTYPE = C. NEXT=0. ITEMS READ* 2. AREA CODE* l . SPECIAL AREA=-0 • SQ.FEFT 

SPECIFIED COEFFICIENT 

T— R = 1 3 00 • H COFF F = C .1 380662E-0? C OFF F ♦AREA* 0 .2 32 2056 F-04 

FL FMF NT* 1 10. TYPF = 5. SUBTYPE = 0. NEXT=I. ITEMS READ* 2. AREA CODE* 2- SPECIAL AREA=-0. SQ.FEET 

SPFCTFIFD COFFFICIFNT 

T-B* 676. H COEFF* C .3 7331 OOF-04 COEFF ♦AREA* 0 .4678 1 2 1 E-06 

Fl FMFNT* 1 10. TYPF* 6 . SUBTYPE=2. NEXT=0. ITEMS RE AD* 14. AREA CODE=I. SPECIAL AREA=-0. SQ.FEFT 

VAR. TEMP. RADI AL ^FAL GAP 

T-B* 9 CO. K— F = C.8707E-05 C-P* 0.2590E 00 MU* 0.2345F-04 NU= 0.8278E-03 R-AV* O.2950E 00 GAP* 0.8375E-02 OMEGA* 1570. 

R— 0= 0.2997E 00 FLAG-C=0. FLAG-I*2. ANN. NBR. = 0. FLO.NBR.* 108. MASS FLOW* 0.5954E-02 

RF* 0.1659F 06 Pp**l/3= C.8869E 00 RE~C= 0.1947F 05 H COEFF* 0 .2 106 143E-01 COEFF *AR E A* 0. 2 6C9898E-03 

HT .GFN .FACTOR* 0.16334E 01 

Fl FMFNT* 94. TYPF=6. SUBTYPE*?. NEXT=0. ITEMS READ=1 4. AREA C00E*1. SPECIAL AR EA=— 0. SQ.FEET 

VAR.TFMP. RADIAL SFAL GAP 

T— B = 10 £?. K-F* 0.91 47E-0 5 C-P* C.2617E 00 HU* 0.2446E-04 NU= 0.4423E-04 R-AV* 0-2550E 00 GAP* 0.4333E-04 OMEGA* 1570. 

R-D= D.2992E 00 FL AG-C *0 . FL AG- T = 2. ANN. NBR. = 1. FLO.NBR.* 0. MASS FLOW* 0.2382E-0? 

RF* 0. 2308E 07 PR^+l/3* 0.8878F 00 RE-C* 0.6754E 07 H COEFF* 0.3748462E 00 COEFF*AR E A* 0. 3002915E-02 

HT .GFN .FACTOR* O.IOOCOF 01 

ELFMFNT* 96. TYPF=6. SUBTYPE*?. NEXT=0. ITEMS READ* 14. AREA CODE=l. SPECIAL AR EA=-0 . SQ.FEFT 

VAR.TFMP.R ADI Al SFAL GAP 

T-B* 975. K-F= 0.8929E-05 C-P* 0.2604E 00 MU* 0.2396E-04 NU= 0.4261F-04 R-AV* 0. 2604E 00 GAP* 0.4333E-04 OMEGA* 1570. 

P-n= 0.2992F 00 FL A G-C =0 • FLAG- 1 = 2. ANN. NBR. = 7. FLO.NBR.* 94. MASS FLOW* 0.2382E-02 

RE* 0.249RF 07 PR+*l/3= 0.8874F 00 RE-C* 0.6754E 07 H COEFF* 0.3657236E 00 COEFF*AREA= 0 . 349074 7E-02 

HT .GFN .FACTOR* 0.10CCCE 01 

FL FMFNT* 98. TYPE = 6. SUBTYPE*?. NEXT=0. ITEMS REA0=14. ARFA CODE=l. SPECIAL AREA=-0 • SQ.FEET 
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VAR.TFMP.RATI AL SFAL GAP 

T— B= 945. K-F= 0.8847F-05 C-P= 0.2599F 00 MU= 0.2377F-04 NU 3 0.4245F-04 R-AV= 0.2663E 00 GAP = 0.4333E-04 OMEGA 3 1570. 

R-R= 0.2992E 00 RAG-00. FLAG~I 3 2. ANN. NBR. = 13. FI O.NBR. 3 96. HASS FLOW 3 0.2382F-02 

RF= 0.26?3F 07 PR**l/3= 0.8F72F 00 RE-C 3 0.6754F 07 H C0EFF= 0.3623054F 00 C0EFF*AREA= 0. 353 5583E-02 

HT.FFN .FACtnR= O.IOOCCF 01 

EL FMFNT 3 100 . TYPF=6. SUBTYP£=2. NEXT=0. ITEMS RE AD=1 4. AREA CODE=l. SPECIAL AREA=-0. SO. FEET 

VAP.T FMP. RADIAL SFAL GAP 

T — R=10 c 0. K-F = C.9090F-05 C-P= 0.2614F 00 MU 3 0.2433E-04 NU= 0.4552F-04 R-AV= 0.2713E 00 GAP 3 0.42C8E-02 OMEGA 3 1570. 

R-n= 0 .2997F 00 FLAG-C = 1. FLAG~I=0. ANN. NBR. = 0. FLO. NBR. 3 98. MASS FLOW 3 0.5954E-0? 

R F= 0. 2539F 07 P P**l/3= G. 8877F 00 RE-C 3 0.4183E 05 H COEFF= 0.2351904E 00 COEFF*AREA= 0. 1670160E-02 

HT .GFN »F ACTOR = 0.IPP66E 01 

COMBINED FLOW ELEMENT* CONDUCTION NO. = 100 FLOW FPCM CCND. ELEMENTS 98 AND 19 

FI FMFNT 3 102. TYPF=6. SURTYPF = 7. NEXT=0. ITEMS PEAD=l4. AREA C0DE=1. SPECIAL AREA=-0 • S <3. FEET 

VAO. TFMP. RACIAL <FAI GAP 

T— R= 10 C5 . K-F = 0.R981E-05 C-P= 0.26C7F 00 ML 3 0. 240SF-04 NU 3 0.5563E-04 R-AV = 0. 2754E 00 GAP 3 0.4333E-04 OJ* GA= 1570. 

P-0= 0.2997F 00 FI AG-C=0. FLAG-I 3 2. ANN. NBR. 3 24. FLO. NBR. =100. MASS FL0W= 0.5954E-02 

RF= 0.2I41F 07 PR** 1/3= 0.8875E 00 RE-C = 0.fc754F 07 H COEFF= 0.3679046E 00 COEFF*AREA= 0. 2652736E-02 

HT .FEN .FACTOR 3 C.100C0F Cl 

FL FMFNT=I 04. TYPF = 6. SUBTYPF=2. NEXT=0. T TFMS RE AD= 14. AREA CODE 3 1. SPECIAL AREA=-0. SO. FEET 

VAfl.TFMP.RADIAL sfal GAP 

T-B= 9 * 5 . K_c- o.-.n?2F-C5 C-P= 0. 2599F 00 MU= 0.2378F-04 NU= 0.1497E-03 R-AV= 0.2796F 00 GAP= 0.43336-04 OMEGA 3 1570. 

R _ n _ 00 P L AG— C = 0. FLAG-I=2. ANN. NBR. = 30. FLO. NBR. = 102 . MASS FLOW 3 0.5954F-02 

PF= 0.P197F 06 PR** 1/3= 0-8872F 00 RE-C= 0.6754E 07 H COFFF= 0.3625071F 00 C06FF*AREA= 0 . 26 53361 E-02 

HT .CFN .FACTOR 3 0.10000E 01 

FL FMF NT= 106. TYPE =6 . SUBTYPF = 2. NFXT=0. I TFMS RFA0=14. AREA C0DE=1. SPECIAL ARE A=- 0. SQ.FEET 

VAP. T FMP. RADIAL SEAL GAP 

T— R= 9^5. K-F= 0.8799F-05 C-P= 0.2596F 00 MU= 0.2366F-04 NU 3 0.8472E-03 R-AV 3 0.2P46E 00 GAP 3 0.42C8E-02 OMEGA 3 1570. 

P-n= 0.7992F 00 FL AG-C = 0. FLAG-I=2. ANN. NBR. = 0. F L O.NBR. 3 1 04. MASS FLOW= 0.5954F-02 

RF= 0. 1 50 IF 06 PR**l/3= C.P871F 00 RF-C= 0.4183F 05 H COEFF= 0 .2275467E-0 1 CQEFF*AREA= 0. 2373432E- 03 

HT .GFN .FACTOR 3 0.17420E 01 

FL FMF NT = 108. TYPF = 6. SUBTYPF=?. NEXT = 0. I TFMS RF AD= 14* AREA CODE=t. SPECIAL ARFA=-0. SQ.FEET 

VAR. T FMP. RADIAL SFAL GAP 

T— B= 915. K— F= C.8746F-05 C-P= C.2593F 00 MU= 0.2354F-04 NU 3 0.8360F-03 R-AV= 0.2900E 00 GAP= 0.42086-02 OMEGA* 1570. 

R— n=0.?99?F QO FL AG-C=0. FLAG-I=2. ANN. NBR. = 0. FL O.NBR. 3 106 . MASS FLOW 3 0.5954F-02 

PF= 0.1579F 06 PR**l/3= 0. 8870F 00 RF-C= 0.4183F 05 H COEFF 3 0 . 23 l 89 10E-0 l COEFF*ARFA= 0.21 12670E-03 

HT.GFN .FACTOP 3 0. 16882F 01 


BOUNDARY COEFFICIENT MATRIX - ( H* A) RTU/F.SEC AND (H*A*T-BI OR W*C-P BTU/SEC 


BNDRY 

C ON C . 

FIXFD FLUX 

FORCFD CONVECTION 


FREE CONVECTION 


RAOIATION 

VAR. TEMP. 

CONVECTION 

FI FM. 

FLFM. 

W 

H* A 

H*A*T-B 


H*A 

H*A*T-B 


H*A H*A*T-E 

H*A 

W*C-P 

1 

1 

0. 

0.64163F-02 

0 ■ 8341 2 E 

01 

0. 

0. 

0. 

0. 

0.31 736F-02 

0.630616-03 

2 

2 

0. 

0.63939F-02 

0.83120F 

01 

0. 

0. 

0. 

0. 

0. 

0. 

3 

3 

0. 

0.95624F-02 

0. 1243 1 E 

02 

0. 

0. 

0. 

0. 

0. 

0. 

4 

4 

0. 

0.946S2F-02 

0. 12310E 

02 

0. 

0. 

0. 

0. 

0. 

0. 

5 

5 

0. 

0.93729F-02 

0.121B5E 

0? 

0. 

0. 

0. 

0. 

0. 

0. 

6 

6 

0. 

0.16735F-01 

0. 21 756 E 

0? 

0. 

0. 

0. 

0. 

0. 

0. 

7 

1? 

0. 

0. 89? BQF— 02 

0.11607E 

0? 

0. 

0. 

0. 

0. 

0.735746-04 

0.15308E 01 

R 

52 

0. 

0. 203 04E-01 

0. 26 396 E 

02 

0. 

0. 

0. 

0. 

0. 

0. 

9 

54 

0. 

0. 83974F-02 

0.1 09 1 7 E 

0? 

0. 

0. 

0. 

0. 

0. 

0. 

10 

57 

0. 

0.13634E-01 

0.1 7724 E 

02 

0. 

0. 

0. 

0. 

0. 

0. 

1 1 

5P 

0. 

0. 14646F-01 

0.19040E 

0? 

0. 

0. 

0. 

0. 

0. 

0. 

12 

?0 

0. 

0. 

0. 


0.41 916E-04 

0.22006E— 01 

0. 

0. 

0.27133E-02 

0.155466-02 

1 3 

’6 

0. 

0. 

0. 


0.22431 E-03 

0.11776F 00 

0. 

0. 

0. 

0. 

14 

41 

0. 

0. 

0. 


0. 23784E— 03 

0.12487E 00 

0. 

0. 

0. 

0. 

I 5 

46 

0. 

C. 

0. 


0.50882E-03 

0. 267 13 E 00 

0. 

0. 

0. 

0. 

i 6 

49 

0. 

0. 

0. 


0.61460E— 03 

0.32267E 00 

0. 

0. 

0. 

0. 

17 

*0 

0. 

0. 

0. 


0.94313E-04 

0*495 14E-01 

0. 

0. 

0. 

0. 

18 

51 

0. 

0. 

0. 


0. 22439E-03 

0. 1 1780 F 00 

0. 

0. 

0. 

0. 

1 9 

4B 

0. 

0. 

0. 


0. 12721E— 03 

0.66783 E-01 

0. 

0. 

0. 

0. 

20 

45 

0. 

0. 

0. 


0.12575E-03 

0.66020 E-01 

0 . 

0 . 

0 . 

0 . 

21 

PI 

0 . 

0 . 

0 . 


0.39600E-03 

0. 20790 E 00 

0 . 

0 . 

0 . 

0 . 

22 

F3 

0 . 

0 . 

0 . 


0.45617F-03 

0.23949F 00 

0 . 

0 . 

0 . 

0 . 

23 

f 5 

0 . 

0 . 

0 . 


0 • 46 74 7 E-03 

0. 24542 F 00 

0 . 

0 . 

0 . 

0 . 

24 

P7 

0 . 

0 . 

0 . 


0.51 766E-03 

0.27177E 00 

0 . 

0 . 

0 . 

0 . 

?5 

EB 

0 . 

0 . 

0 . 


0.78759E— 03 

0.41349E 00 

0 . 

0 . 

0 . 

0 . 

26 

P6 

0 . 

0 . 

0 . 


0.68920E-03 

0. 36I83E 00 

0 . 

0 . 

0 . 

0 . 

27 

P4 

0 . 

0 . 

0 . 


0.67 255E-03 

0.35309 F 00 

0 . 

0 . 

0 . 

0 . 

?fl 

P? 

0 . 

0 . 

0 . 


0.58384E-03 

0.30652E 00 

0 . 

0 . 

0 . 

0 . 

29 

77 

0 . 

0 . 

0 . 


0.42925 E-03 

0.22535E 00 

0 . 

0 . 

0 . 

0 . 

30 

7? 

0 . 

0 . 

0 . 


0 • 45 258E-03 

0. 2 3760 E 00 

0 . 

0 . 

0 . 

0 . 

31 

65 

0 . 

0 . 

0 . 


0. 10060E-03 

0 .528 14E-0 1 

0 . 

0 . 

0 . 

0 . 

3? 

66 

0 . 

0 . 

0 . 


0. 10 144E-03 

0. 532 54 F- 01 

0 . 

c. 

0 . 

0 . 

33 

67 

0 . 

0 . 

0 . 


0.12575E-03 

0.66017F-01 

0 . 

0 . 

0 . 

0 . 

34 

73 

0 . 

0 . 

0 . 


0. 14541E-03 

0.763 40 F-Ol 

0 . 

0 . 

0 . 

0 . 

35 

7B 

0 . 

0 . 

0 . 


0 ■ 25756E-03 

0 . 13522 F 00 

0 . 

0 . 

0 . 

0 . 

36 

79 

0 . 

0 . 

0 . 


0.12576E-03 

0.66022E-01 

0 . 

0 . 

0 . 

0 . 

37 

PO 

0 . 

0 . 

0 . 


0.29601E-03 

0.15540E 00 

0 . 

0 . 

0 . 

0 . 

38 

75 

0 . 

0 . 

0 . 


0.2155BE-03 

0. 1 13 18 F 00 

0 . 

0 . 

0 . 

0 . 

39 

70 

0 . 

0 . 

0 . 


0. 2375 3E-03 

0.124706 00 

0 . 

0 . 

0 . 

0 . 

40 

6? 

0 . 

0 . 

0 . 


0.25454E-03 

0.13363E 00 

0 . 

0 . 

0 . 

0 . 

41 

61 

0 . 

0 . 

0 . 


0 • 4191 8E—0* 

0. 22007E— 01 

0 . 

0 . 

0 . 

0 . 
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42 

1 1 

0. 

u. 

0. 

0. 

0. 

0. 

o. 

0.35758E-04 

3.433 11 E-3 3 

43 

17 

0. 

0. 

0. 

3. 

0. 

0. 

0- 

0. 3L952E-04 

3. 47943E-03 

44 

16 

0. 

0. 

0. 

0. 

0. 

0. 

0- 

0. 41609E-04 

3. 47849E-03 

45 

15 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 2V351E-04 

0.47743E-03 

46 

20 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0.22995E-04 

3.47740E-03 

47 

19 

0. 

0. 

0. 

0. 

0. 

0. 

o. 

0. 22 872E-04 

3.4771 3E-33 

4d 

7 

0. 

0. 

0. 

0 . 

0. 

0. 

0. 

0.3S508E-02 

3 • 62599E-3 3 

49 

13 

0. 

0. 

0. 

0. 

0. 

0. 

0 . 

0. V92 03E-02 

3.627 44E-33 

50 

24 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 614 15E-02 

3. 15591E-32 

51 

92 

-0. 10299E GO 

0.232216-04 

0. 30L87E-01 

0. 

0. 

0. 

o. 

0. 

0. 

52 

1 10 

0. 

O.46701E-O6 

0.25028E-03 

0. 

0. 

0. 

o. 

0.2S099E-03 

3. 15424E-32 

53 

94 

0. 

0. 

0. 

3. 

0. 

0. 

o. 

0. 30029E-32 

3 . 62337E-33 

54 

96 

0. 

0. 

0. 

0 . 

0. 

0. 

0. 

0.3V907E-32 

3.62325E-33 

55 

98 

0. 

0. 

0. 

0. 

u . 

0. 

o. 

0. 3535 SE-02 

3.6190BE-33 

56 

1U0 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 1S732E-02 

3. 15 55 IE-32 

57 

102 

0. 

0. 

0. 

3. 

0. 

0. 

0. 

0.2S527E-32 

3.I5523E-32 

50 

104 

0. 

0. 

0. 

0. 

0. 

0. 

o. 

0.2S53VE-02 

3.15476E-02 

59 

106 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0.2373VE-03 

3. 15457E-32 

60 

100 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 21 1 27E-03 

3.15433E-32 

61 

91 

0. 

0. 

0. 

0. 

0. 

0. 

o. 

0. 

3. 

62 

93 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3. 

63 

95 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3. 

64 

97 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0 i 

65 

99 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3. 

66 

101 

0. 

0. 

0. 

3. 

0. 

0. 

0. 

0. 

3. 

67 

103 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3 • 

60 

105 

0. 

0. 

D. 

0. 

0. 

0. 

0. 

0. 

3. 

69 

107 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3. 

70 

109 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3. 

71 

09 

0. 

0. 

0. 

3. 

0. 

0. 

0. 

0. 

3. 

72 

90 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3. 


VAK. TEMP.CJNVECT ION 0 I MENS 1 ONL ESS COEFFICIENTS ( H*A/W*C-P ) 

I 1, 5.032 7E GO 2 12, 1.5299E-01 3 33, 1.7453E 00 4 II, 7.6561E-02 5 17, 6.6650E-02 6 IS * 8.S960E-32 

7 15, 5 • LOG 76- 02 0 20, 4.8157E-32 9 19, 4.7937E-32 10 7, 5.8320E 00 LI 13, 7.8418E 00 12 2V, 3.9393E 33 

13 110, l • 6S2 IE-0 1 14 94, 4.8172E 00 15 96, 5.6280E 00 16 90, 5.71LOE 00 17 100, 1.0733E 00 IB 102, 1.7389E 00 

19 104, 1.7145b 00 20 106, 1.5355E-01 21 : 1.3695E-0L 


FLUID SHFAR HEAT GFNFRATIQN, 0TU/SEC 

1 1, 3.0329E-C2 ? 12. O. 3 30, 8.1501F-02 4 11, 0. 5 17, 0. 6 16, 0. 

7 15. 0. 6 20. 0. 9 19, O. 10 7, 3.6591E-02 11 13, 3.6094E-03 12 24, 2.2674E-03 

13 110. 4.8B3 7F- 04 14 94, 2.8952F-02 15 96, 3.5241F-02 16 98, 3.7368E-02 17 100, 5.6112E-04 18 102, 2.9926E-0 2 

19 104. ■* .C920F— C2 20 106, 8.1461E-04 21 108, 7.1239F-04 

DATA T I MF= 0. 1 3f€ MINUTFS 

CONDUCT .FLFMFNT NO. 1 VAR . TE MP . F LF MF NT NO. 1 FLOW COOE= 2 

C.PNOUfT .El. FMFNT ND . 12 VA P . TFMP . E LE ME NT NO. 2 FLOW COOF= 1 

C.nNOU CT . FL FMFN T NO. 30 VAR. TFMP. ELF MF NT NO. 3 FLPW COOE= 2 

CPNOUCT .FI FMFNT NO. 11 VA R . TEMP . FLE PE NT NO. 4 FLOW CODE= 1 

CONDUCT. FLEMFNT NO. 17 VAR. TEMP. FLFPENT NO. 5 FLOW COOE= 1 

CONDUCT. FLEMFNT NO. 16 VAR . TFMP. p| F ME NT NO. 4 FI PW CODE= 1 

CONOUCT.EI FMENT NO. 15 VAR . TEMP. ELF ME NT NC. 7 FLCW COOE= 1 

CONDUCT .FLFMFNT NO. 20 VA R . TFMP . E LE MF NT NO. 8 FLOW CnOE= 1 

CONTUCT. FLFMFNT NO. 19 VAR . TEMP. ELF MFNT NO. 9 FLOW CODE= 1 

CONDUCT. FLFMENT NO. 7 VAR. TEMP. ELEMENT NO. 10 FLOW COOE= 2 

CONDUCT. FLCMFNT NO. 13 VA R . TF MP . E 1 E ME NT NO. 11 FLCW CODE= 2 

CPNOUCT. FL FMFNT NO. 24 VAR . TFMP . 6 LF ME NT NO. 12 FLOW CODE= 2 

CONDUCT. FLFMFNT Ml . 110 VAR . TEMP . ELE ME NT NO. 13 FLOW CODE= 2 

CONDUCT. FLEMFNT NO. 94 VAR . TFMP . E LE ME NT NO. 14 FLOW CODE* 2 

CONDUCT. FLFMFNT NO. 96 VAR . TEMP. ELE MF NT NO. 15 FLOW COOE= 2 

CONDUCT. El FMFNT NO. 98 VAR . TFMP . ELE MENT NO. 16 FLOW CODE= 2 

COMBINED FinwS FLFMFNT, COMBINATION NO. I CCNDUC T. NO .= 100 BNDRY • NO . = 17 PREVIOUS ELFMENT VAR . TEMP . NOS- =* 16 ANO 9 

CONDUCT. ELEMENT NO. 102 VAR . TEMP .E IE ME NT NO. 17 FLOW C0DE = 2 

CONDUCT. FI. FMFNT NO. 104 VAR . TFMP . ELE ME NT NO. 18 FLOW CODE= 2 

CONDUCT. FI FMFNT NO. 106 VAR . TEMP . E LE ME NT NO. 19 FLOW CODE= 2 

CONDUCT. FL FMENT NO. 108 VAR . TEM P. ELE ME NT NO. 20 FLOW CODE= 2 
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UNVtKSENCt CRITERIA, I T ER AT I UN= 100 


COND. NO. 

NbH TEMP. 

OLD TEMP. 

DIFFERENCE 

MAX. 3 IFF. 

1 

0. 12071 08 E 04 

0. 1207189E 04 

-0 * 9 00 2686 E- 03 

0.5000000E-01 

2 

0.1262034E 04 

3 . 12620666 34 

-3. 1525879E-02 

0.50000006-01 

3 

0.1279193E 04 

0.1279195E 34 

-0 • 2 0446 78 E- 02 

0. 5000000E-01 

4 

O.1205909t 04 

O.L205991E 34 

-3 .215L489E-02 

0. 5003000E - 01 

5 

O.129O909E 04 

0. 1290990E 34 

-3. 1678467E-02 

0. 5000 OOOE-Ol 

6 

0* l 2983505 04 

3 . 1 2980 50E 34 

-3. 3 66 2109 E- 03 

0. 5000 OOOE-Ol 

7 

0.1047630E 04 

0 • 1347635 6 34 

-3 * 5 1574 7 IE- 02 

0.5000000E-01 

0 

0* 1 1 2 007 9 E 04 

3.1120338E 34 

-3.91857916-02 

0.50000005-01 

9 

0.1173756E 04 

0.11737686 34 

-3.1 26 9531 E- 01 

0. 50000006-01 

10 

0. L 2095 3 OE 04 

0. 1239544E 04 

-0.1 3626106-01 

0. 5003000E- 01 

IL 

0.123949 IE 04 

0 ■ 12395 02 E 04 

-3 .10986336-01 

0. 5000000E-01 

12 

0.1283117E 04 

3.12831206 34 

-3.3 25 01226- 02 

0. 5000000E -01 

13 

0. 984 1284 E 03 

0.9841373E 03 

-3. 8911133 E- 02 

0. 5000000E-01 

14 

0* 1 0431 3 3 E 04 

3 ■ 13 4315 IE 04 

-0. 1814270E-01 

0. 50000005-01 

15 

0.1095047E 04 

3.13950 72E 34 

-3 , 2456 665E-0 1 

0. 50000005-01 

15 

0.U38706E 04 

0.11 3 873 3E 34 

-3. 265960 7E-01 

0. 5000000E-01 

17 

0 . 1 1 75a5 4E 04 

0 • 1 L 75b77E 34 

-3 * 2363586E- 0 1 

0. 50000005-01 

18 

0. 1211240E 04 

0. 1211257E 04 

-0.1673889E-01 

0. 5000000E-01 

19 

0.10134736 04 

0.13 135 04E 04 

-0 • 3 1 54755E- 0 1 

0. 50000005-01 

20 

0. 1 0447 90 E 04 

3 . 1J44827E 34 

-0.3697205E-01 

0. 50000005-01 

21 

0.1084893E 04 

3.1384932E 34 

-3 • 3944397E-0 1 

0. 50000005-01 

22 

0.1124L59E 04 

0. 1124195E 34 

-3 . 354 7668E-0 1 

0. 50000006-01 

23 

0. 1 1 63354E 04 

0.116338 IE 34 

-3 .26245126-01 

0. 50000006-01 

24 

0.9142733E 03 

0.91 42952 E 03 

-0.21888736-01 

0. 50000005-01 

25 

0.95306596 03 

0.9631024E 33 

-3.36514286-01 

0. 50000006-01 

26 

0. 1 0026 95 b 04 

0.13327426 34 

-0.4 72 3 3 5 8 E- 01 

0.50000005-01 

27 

0.1043344E 04 

3 . 13 4339 5 E 34 

-0*506 59 18 E- 01 

0.50000005-01 

28 

0.10B4J42E 04 

0. 10843386 34 

-3 .4553223E-01 

0.50000006-01 

29 

0. 1127134E 04 

0. 1 1 2 7 1 68 E 34 

-3 . 335 3882E- 0 1 

0.50000005-01 

30 

0.909L553E 03 

0.90919206 03 

-3 . 257 5684E-0 1 

0. 5000000E-01 

31 

0.Q3J02&4E 03 

0.9330827E 03 

-3 • 4633 33 IE-0 1 

0. 50000005-01 

32 

0.9659154E 03 

0.9&S9751C 03 

-0 • 596 1 6Q9E-0 1 

0. 5000000E-01 

33 

0.10041186 04 

0.130419 IE 34 

-3 .63323976-01 

0. 5000000E-01 

34 

0.1045186E 04 

0.13452 42E 34 

-3 . 5645 752 E- 0 1 

0. 5000000E- 01 

35 

0.1092378E 04 

0.L092418E 34 

-3 .4060364E-01 

0. 5000000E - 01 

36 

0.89631746 03 

0 • 896 3848 E 03 

-3 *673599 2E-0 1 

0. 50000005-01 

37 

0.9249237E 03 

0.925031 7 E 33 

-0.78033456-01 

0. 5000000E-01 

38 

0. 95 845 5 4 E 03 

0 . 95 85373E 03 

-0.8 09 4788 6- OL 

0. 50000005-01 

39 

0.99843946 03 

3.99851086 33 

-0,7141 L13E-01 

0.5000000E-01 

40 

0.10517776 04 

0.13518256 34 

-0 .4855347E-01 

0. 5000000E -01 

41 

0. 85 763 35 £ 03 

0 . 8> 77 148 E 33 

-0.81291 20 E- 01 

0. 50000005-01 

42 

O.B03589OE 03 

0.88369096 33 

-0.10 19058E 00 

0. 5000000E-0I 

43 

0.90705086 03 

0.9371563E 33 

-0 ■ 1055069E 00 

0. 5000 OOOE - 01 

44 

0. 941 13 9 9E 03 

0.941 2338E 33 

-0.9381104E-01 

0. 50000005-01 

45 

0.10115816 04 

3.1311636E 34 

-0.5464172F-01 

0. 5 0 0 0 0 0 05 - 01 

46 

0.814 92 816 03 

3.81508966 33 

-0*16149 14E 00 

0. 5000000E-01 

47 

0.81762136 03 

0.8177812E 33 

-0 • 1 599 426E 00 

0. 5000000E-01 

48 

0.81869056 03 

0.81904736 03 

-0 . 1567459E 00 

0. 50000006-01 

49 

0.812 762 3 E 03 

0.8 1292 1 8E 03 

-0 . 1 5949256 00 

0.5000000E-01 

50 

0.8151225E 03 

0.81528046 33 

-0 • 1 57936 IE 00 

0. 50000005-01 

51 

0, 81 597 3 l E 03 

3.81612806 33 

-0*1 548 69 IE 00 

0. 50000005-01 

52 

0.123 983 76 04 

0.12398456 34 

-0.79193126-02 

0.5000000E-01 

53 

0.121 00 87E 04 

0.1210098E 34 

-3. 1 136 7806-01 

0. 5000 OOOE-Ol 

54 

0.12108266 04 

0.1210B37E 34 

-0.1 10 7788 E- 01 

0. 5000000E-01 

55 

0.11851086 04 

3.1185122E 34 

-0. 13717656-01 

0.5000000E-01 

56 

0. 1 L 8994BE 04 

3 . 11899616 34 

-0.13259896-01 

0. 5000000E -01 

57 

0.12 12945E 04 

0.12129586 34 

-3. 12725836-01 

0. 5000000E- 01 

58 

0.12129106 04 

0.1212927E 34 

-3 .1690674E-01 

0. 5000 OOOE- 01 

59 

0.11630286 34 

3.11 6305 9E 34 

-0.30273446-01 

0. 5000000E-01 

60 

0.L131052E 04 

0.11311066 34 

-0.43960576-01 

0. 5000 OOOE-Ol 

61 

0.11154386 04 

3.11 15537E 34 

-0. 48675546-01 

0. 5000000E-01 

62 

0. 1 l 05840E 04 

0.110 59 39E 34 

-3.49L7908E-01 

0. 5000000E-01 

63 

0.11 59076E 04 

0. 1159092E 34 

-3. 15396126-01 

0.5000000E-01 

64 

0.116 698 4E 04 

0. 1166998E 34 

-3 .14633186-01 

0. 5000 OOOE-Ol 

65 

0. 1 202627E 04 

0 • 120 264 l E 04 

-3.13748176-01 

0. 50000005-01 

66 

0. 12 0444 3 E 04 

0 . 123446 IE 34 

-0. 18066416-01 

0. 5000000E-01 

67 

0. 1 163022 E 04 

0.1 1 63053E 04 

-3.30303966-01 

0.50000005-01 

68 

0. 1 1 2 8745 E 04 

3 • 1 1 2879 0 E 34 

-3.43914796-01 

0. 5000000E-01 

69 

0.1114913E 04 

0.11149616 34 

-0.48004156-01 

0. 5000 OOOE - 01 

70 

0. 1 1O5042E 04 

0. 11058906 34 

-3.4837036E-01 

0. 5000000E-01 

71 

0. L12 32 32E 04 

3.11232496 04 

-3 .16906746-01 

0. 50000005-01 

72 

0.11235576 04 

0.U23583E 04 

-3.15975956-01 

0. 5000000F -01 

73 

0.11216886 04 

0.LL21733E 04 

-3.45181276-01 

0. 5000 OOOE-Ol 

74 

0.1113244E 34 

3.11132926 34 

-3. 4 74 3958 6- 01 

0. 50000005-01 

75 

0.11051526 04 

0.11351996 34 

-3.47515876-01 

0. 50000005-01 

76 

0.1087554E 04 

0. 1087571E 34 

-3 . 1757 81 3E -01 

0. 5000000E -01 

77 

0.10863316 04 

0.1386048E 04 

-3.1 666 260E- 0 1 

0. 5000000E-01 

78 

0.11183456 04 

0,11 1 8390 E 34 

-3.44952396-01 

0. 50000005-01 

79 

0.11L1B526 04 

3.1L11898E 34 

-0 . 4653 93 IE- 0 l 

0.50000005-01 

80 

0. 1 1 04255E 04 

0.11343026 34 

-3.46463016-01 

0. 50000005-01 

81 

0.10553706 04 

0.1355 38 7E 34 

-3 . 1 7227 1 7E-0 1 

0.50000006-01 

82 

0.10539326 04 

0 . 10 53999 E 34 

-3 . 167003 7E -01 

0. 5000000E -01 

83 

0. 102604 7E 04 

0.1326064E 34 

-3.1692200E-01 

0. 50000005-01 

04 

0.10247456 04 

3.1324763E 34 

-3.16494756-01 

0. 50000005-01 

85 

0.1005677b 04 

0. 1305693E 04 

-3.16578676-01 

0. 5000000E-01 

86 

0.10044076 04 

0. 13044236 04 

-3 .16181956-01 

0. 50000005-01 

87 

0.99531526 03 

0.99533136 33 

-3.1615906E-01 

0. 5000000E-01 

88 

0.99393586 03 

0 • 99 395 1 6 E 03 

-3.1 577759E-0 1 

0. 5000000E -01 

89 

0. 1 1 07300E 04 

0.1107000E 34 

3. 

0.5000000E-01 

90 

0.10600006 04 

0. lO&OOOOE 04 

3 . 

0.5000000E-01 

91 

0.8710030E 03 

0.87 10000E 33 

3. 

0. 5000000E-01 

92 

0. 87262 2 86 03 

3 .8 7 2622 8E 33 

-3.7 629 3 95 E- 05 

0. 5000000E-01 

93 

0. 86 9 0000 E 03 

0.85 90300 E 33 

3. 

0. 5000000E-01 

94 

0. 8791839E 03 

0.879L839E 33 

3. 

0. 5000000E-01 

95 

0. 8560030b 03 

0 ■ 85 60000E 33 

3. 

0. 50000005-01 


96 

0. 0 7602 2 BE Oi 

0.3760228E 03 

-0.7629395E-05 

0. 50000006-01 

9 7 

0*864 DO OOF 03 

D.85400CQE 03 

0. 

0. 5000000E-01 

90 

0.07121506 03 

0.8712L50E 03 

-0. 1525879 E- 04 

0. 5000000E-01 

99 

0. 06 l 00 JOE 03 

0.85 l 00 00 E 03 

0. 

0. 5000000E - 01 

100 

0.8661229E 03 

0.8561230E 03 

-0.22888 18E-04 

0. 5000000E-01 

101 

O.ObOuODOfc 03 

0 ■ 85 OOOOOE 03 

0. 

0. 50000006-01 

102 

0* 867 3022 fc 03 

0.8673022E 03 

-0.15258 7 9E-04 

0. 5000000E -01 

103 

0.8580030E J3 

3.8580000E 03 

0. 

0. 5000000E-01 

104 

0.8648023E 03 

0.8648024E 03 

-0.7629 39 5E-05 

0.5000000E-01 

lot 

0. 8560000E 03 

0*8560000 E 03 

0. 

0.50000006-01 

10b 

0.85999A2F 03 

0 • 8599932E 03 

0 . 

0.5000000E-01 

107 

0. 86500006 03 

0 . 85 50000E 03 

0. 

0. 5000000E - 01 

100 

0.8572584E 03 

0.85725846 03 

0. 

0. 50000006-01 

109 

0. 85 300 UOE 03 

0.85 30000E 03 

0. 

0. 5000000E-01 

no 

0. 85546 8 7 E 03 

0.8 554687 E 03 

0. 

0. 50000005-01 


CONVERGENCE CRITERION 4CHEIVEJ, ITERATIONS- 131 


ELEMENT ( AT CENTER] TFMPE MATURES (OFOPEFS F ) > ut:RATIONS = 130 


1, 

I.2073900F 

C3 

2. 

1.26219966 

03 

3, 

1 .2792890F 

03 

4, 

1.2860528E 

03 

5* 

1.291025RE 

03 

6 , 

1.29B0559E 

03 

7, 

1.048521 IF 

03 

a. 

l.l 20R364E 

03 

9, 

1. 17434R3F 

03 

10, 

1. 209935 IE 

03 

11, 

1. 239736 3E 

03 

12, 

l • 28 3 162 36 

03 

13. 

5.F575R97F 

C2 

14. 

l • 04^-433 5F 

03 

15, 

1.0960944F 

03 

16, 

1. 1394177C 

03 

17, 

1 . 1759955E 

03 

18, 

1.21129946 

03 

19, 

1 . 0149709F 

C3 

20, 

l . 045R352E 

03 

21, 

1.0B536R66 

03 

22, 

1. 1243037F 

03 

23, 

1. 1633123E 

03 

24, 

9.16540686 

02 

25. 

5.644 207 1 F 

C2 

26, 

1.003393 8E 

03 

27 , 

1 . 04354 96E 

03 

28 ♦ 

1.0839858E 

03 

29, 

l . 12697716 

03 

30, 

9.10708286 

0? 

31. 

9.33R5793F 

C2 

32. 

9.661 765 86 

02 

33, 

1.0039RP6F 

03 

34, 

1.044BRB1F 

03 

75, 

1 • C920B69E 

03 

36, 

8.96238366 

02 

37, 

9. 245743 ?F 

C2 

38, 

5. 57P6542F 

02 

39, 

9 .9780642 p 

0? 

40, 

1.05 1 3206F 

03 

41, 

R.6706945F 

02 

42, 

8.825 9656F 

02 

4 3, 

5.C5RR249F 

C? 

44. 

9-4004122E 

02 

45 , 

1.0109R61F 

03 

46, 

8. 1260R96E 

02 

47, 

8 . 153 2953E 

0? 

4 8, 

8. 16651906 

02 

49. 

P. 1 0471 8 RF 

C2 

50, 

B. 128611 7E 

0? 

51 , 

8.1 376455 F 

02 

52, 

1.239740BE 

03 

53, 

l -2099419F 

03 

54, 

1 .21067886 

03 

55. 

1.1 849299F 

C3 

56. 

1. 1 35766 9F 

03 

57, 

1.2127559F 

03 

58, 

1.2126390F 

03 

59, 

1- 1625276E 

03 

60, 

1. 13032726 

03 

61. 

1.1 146724F 

C3 

62, 

1. 1P50158F 

03 

63 , 

1. 158R73RF 

03 

64, 

l. 1667836F 

03 

65, 

1 . 2024 2? 8F 

03 

66, 

l • 204 154 l E 

03 

67. 

1 .162521 5F 

C3 

6R » 

1. 12P012CE 

03 

69, 

1 .1 141OR0F 

03 

70, 

1. 1050314F 

03 

71, 

1 . 12300R0F 

03 

72, 

1.12335256 

03 

73. 

1.1 2P9323F 

C3 

74. 

1. 1124499E 

03 

75, 

1.1 04356 2 F 

03 

76, 

1.08731BBE 

03 

77, 

1 .C858072E 

03 

78, 

1 .11759326 

03 

79. 

1.1 11072RF 

03 

BO, 

1. 1034779F 

03 

81 , 

1.055 1 364F 

03 

82, 

1.053756 OF 

03 

83, 

l .0258155E 

03 

84, 

1 . 0245209E 

03 

85. 

1.00544R7E 

C3 

R6 » 

1.004] R52E 

03 

R7 , 

9.9 509 299 F 

02 

R8 » 

9.937L971F 

0? 

P9, 

1. 10 700006 

03 

90, 

1.06 OOOOOE 

03 

91. 

F. 7 09 9 99 9 F 

r? 

92, 

8. 728C04 7E 

02 

93, 

R.6R99999F 

02 

94, 

8. 796955PF 

02 

5 5, 

8.6600000E 

02 

96, 

8.77112916 

02 

9 7. 

P.6400000E 

C2 

98, 

R.73P3481E 

02 

99, 

B . 60999 99 F 

02 

100, 

R.7193322E 

02 

101, 

8. 5999999F 

02 

10?, 

8.7032722F 

02 

103. 

P.57 99999F 

C2 

104. 

R. 663988 8E 

02 

105, 

R . 5 5999 99 F 

02 

106, 

8.6073503F 

02 

107, 

8.55000006 

0? 

l C 8 , 

fl • 57 6 1 78 IF 

02 

109. 

P . 5 300000F 

C2 

1 10, 

8. 556591 7F 

0? 














E LE ME NT ( A T CENTER) TEMPERATURES (DEGREES FI, I T 6* AT I 3N S = 131 


1, 

1.20 7 172 5E 

03 

2, 

1. 262 03736 

03 

3, 

1.2791568E 

03 

4, 

1.2859512E 

03 

5, 

1 .2909592E 

03 

5, 

1 • 2 980 4 3 3E 

03 

7, 

1.04753906 

03 

b , 

1. 11991716 

03 

9 , 

l. 1 7353 14E 

33 

13, 

1.2092B99E 

03 

11, 

1.23929 73; 

03 

12, 

1.283 05 93E 

03 

13, 

9. 83970 >36 

02 

14, 

1 . 042 0 1 2 9E 

03 

15, 

1 .0946141E 

03 

16, 

1. 1382373E 

03 

17, 

1 .17523746 

03 

18, 

1 . 2 L 09470E 

03 

19, 

1.01 2 915 06 

03 

20, 

1.04413716 

03 

21 , 

1.08419716 

03 

22, 

1. 1235345E 

03 

23, 

1 .16289386 

03 

24, 

9.139855 5E 

02 

25, 

9.624 209 56 

02 

2 6, 

l. 00106066 

03 

27 , 

1 .0424498E 

03 

20 , 

1.08 32400E 

03 

29, 

1 .12654526 

03 

30, 

9.09711 1 IE 

02 

31, 

9. 322 17366 

02 

32, 

9. 64 862 2 BE 

02 

33 , 

1. J030300E 

03 

34, 

1 . 04419 05 E 

03 

35, 

1.0916644; 

03 

36, 

8. 95 1 2 6 3 7 E 

02 

37, 

9. 23544796 

02 

38, 

9. 57026716 

02 

39 , 

9. 97179 456 

02 

40, 

I.O5O9220E 

03 

41, 

0.6619644E 

02 

42, 

8. 3178763E 

02 

43, 

9.05186256 

02 

44, 

9. 39483406 

02 

45, 

1 .3 136195E 

03 

45, 

8. 12071 72E 

02 

47, 

8,147923 5E 

02 

48, 

3.151131 4E 

0? 

4 9, 

8.099413 3E 

02 

50, 

0. 12329016 

02 

51 , 

8. 1323395F 

02 

52, 

1.2 3970 O0E 

03 

53, 

1 .209B909E 

03 

54, 

1.2105354E 

03 

55, 

1. 18467196 

03 

5b, 

1. 1 89 7 1 91 E 

03 

57, 

1.21272526 

03 

50, 

1.21 26 1756 

03 

59, 

1.1625043; 

03 

50, 

l . 1303005E 

03 

61, 

1. 114644 36 

03 

62, 

1. 1 04 9d 8 36 

03 

63 , 

1 . 15Q8138E 

03 

64, 

1. 16 67 3 1 5E 

03 

65, 

1.20238936 

03 

5 5, 

1.2041305E 

03 

67, 

i. 16249796 

03 

08, 

1. 1 2 79 854 E 

03 

6? , 

1 . 1 143814E 

03 

70, 

l. 10 5004 1 E 

03 

7 1, 

l . 12294 026 

03 

72, 

1 .1232921E 

03 

73, 

1 . 12090566 

03 

74, 

1. 11242276 

03 

75 , 

1 . 13432926 

03 

76, 

1.0872511E 

03 

77, 

1.0857446; 

03 

78, 

l . 1 L 75 6 65 6 

03 

79, 

l. Ill C45 9fc 

03 

80, 

1.10345126 

03 

81 , 

1 .3553741E 

03 

82, 

1 .0536952E 

03 

83. 

1 .0257 570; 

03 

94, 

1. 0244530E 

03 

85, 

1.00 5 392 7E 

03 

80 , 

1.0041295b 

03 

07 , 

9.9503328E 

02 

88 » 

9.93665206 

02 

B9, 

1 .1070000; 

03 

90, 

l. OSOOOOOE 

03 

91, 

8. 7099999fc 

02 

92 , 

8. 72622 756 

02 

93 , 

0.58999996 

02 

94, 

8.79183856 

02 

95, 

8.6600000E 

02 

95 , 

8.7502 273E 

02 

97, 

8.64000006 

02 

9a, 

8. 7121484E 

02 

99 , 

0. 60999996 

02 

103, 

8.66 12252E 

02 

10 1, 
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APPENDIX D 


HEAT TRANSFER COEFFICIENT EQUATIONS 

The equations used to calculate boundary heat transfer coefficients are mostly ob- 
tained from books by McAdams (ref. 9) and by Kreith (ref. 10). 


Radiation 

The radiation equation is based on information from chapter 4 of Kreith and chap- 
ter 5 of McAdams. The basic expression is in terms of the heat flux but can also be ex- 
pressed as 


h = — 0 — 

A AT 

The number of properties required is reduced by substitution of others as explained in the 
references. 


Liquid Film Cooling 

The gravity-flow liquid film cooling equation can be found in chapter 9 of McAdams. 
Strictly, it applies only to flow over horizontal cylinders. However, in the program, it 
is also used for flow down vertical surfaces. 


Other Free Convection Coefficients 

The two other free convection subtype expressions are obtained from a study of chap- 
ter 7 of McAdams and of chapter 7 of Kreith. They are actually for nonenclosed spaces, 
but are considered acceptable for two reasons. First, the boundaries at which free con- 
vection occurs are far enough from the parts of interest (seal plate and nosepiece) that 
they have little influence. Second, the knowledge of coefficients in enclosed spaces ap- 
pears to be limited and to be expressed in terms of overall coefficients. 

Comparison of the two conditions (table I) shows the equations for turbulent flow are 
identical and those for laminar flow differ only in the constant used. Transition from 

Q 

laminar to turbulent flow is considered to occur at GrPr = 10 (ref. 9). 
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Use of Local Coefficients 


The expression for free convection at vertical surfaces is for a local coefficient in- 
stead of an average as is the case for the other cases already mentioned. 

The program deals with small sections of boundary surfaces. There are few, if any, 
places where a fully established velocity profile (e.g. , in duct flow) or boundary film 
(e.g. , in flat plate flow) exists in seal applications. In the problem of appendix C, as an 
example, the use of a local coefficient is generally preferable. 

While free convection around a horizontal cylinder should properly be considered a 
circumferential flow, the assumption of an axisymmetric temperature distribution implies 
a uniform coefficient around the circumference as well as along the cylinder. 

The conversion of average coefficients to local coefficients is based on text and equa- 
tions in chapter 6 of Kreith. 


Flat Plate 

These coefficients and the critical Reynolds number for laminar turbulent transition 
are obtained directly from chapter 6 of Kreith's book. 


Duct Flow, Liquid 

For the most accurate result, the equation used for duct flow depends on the shape of 
the cross section and on the viscosity of the fluid. However, McAdams and also Jakob 
(ref. 11) indicate that the use of pipe flow equations with an equivalent (hydraulic) diam- 
eter is acceptable. 

In the example of appendix C, the only fluid flowing in a duct is a synthetic oil with a 
viscosity in the range covered by McAdams recommended equation for turbulent flow. 

The laminar flow expression as well as that for turbulent flow are taken from chapter 9 
of McAdams. 

McAdams' data indicate that transition from laminar to turbulent flow can be con- 
sidered as occurring over the Reynolds number range between 2500 and 7500. Kreith and 
also Jakob agree with the overall trend during transition. The equation presented in 
table I seems to be an acceptable approximation of the coefficient in this range. 
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Concentric Cylinders 


This expression is obtained from reference 12. The experimental conditions are 
similar to those found in the problem of appendix C. One assumption is made in the 
problem; the axial flow has a negligible effect on the boundary film. This seems reason- 
able for the high rotational speed and the low axial flow velocities. Therefore, the en- 
trance effect term is dropped. 


Rotors 

The equations for the remaining two flow types are derived from the results of ref- 
erence 13. While the completely shrouded disk used does not duplicate the flow in a seal 
gap, no work applicable to laminar flow (the most probable regime in a seal gap) has yet 
been found. Also the derivation involves using the Reynolds analogy so the resulting ex- 
pressions are only rough approximations. 

In applying the results, the assumption is made that one or the other of the following 
situations exists: 

(1) The clearance between the side of a rotating part and the housing or other sta- 
tionary part is so small that the boundary layers are merged at all times (seal gap). 

(2) The clearance is large enough that the boundary layers are separate at all times 
(sides of rotors). This should be a good approximation of the situation in the problem of 
appendix C except for the seal gap. 

The laminar to turbulent transition for the sides of rotors case occurs at a Reynolds 
number of 158 000. For the seal gap case, the critical Reynolds number is a function of 
the ratio of gap clearance to rotor outer radius, that is, 


Re c = 


366.34 
1019 


Derivation of Equations for the Rotor Coefficients 

The Reynolds analogy (ref. 14) 


CJRePr 1 / 3 

Nu = — 

2 


(Dl) 
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requires an expression for the wall shear stress T t because the skin friction coefficient 



(D2) 


The torque coefficients obtained in reference 13 are, in effect, functions of Tj.; 

that is 


C M = 


877 


2 5 
pco a 


f r\ dr 
•'n 


(D3) 


A general form of the equation for shear stress in a pipe (ref. 14, p. 509) 

r t = Ap v u w rW (D4) 

is used to obtain the following expression for the torque coefficient: 

(D5 ) 

M 3 + X , 2-w 2-x 
cu a 

This is matched with an empirical equation from reference 13 to evaluate the constant A 
and the coefficients of the properties. Substitution in equations (D4), (D2), and (Dl) gives 
an expression for the Nusselt number. The heat transfer coefficient h is a factor of the 
Nusselt number; that is, 


k 

The free stream velocity V M must still be evaluated. For the seal gap with merged 
boundary layers, V M = rci>/2. For the case of separated boundary layers, the difference 
between the free stream velocity and the rotor velocity is required for a rotating part. 

For a nonrotating part, of course, the difference is the free stream velocity = /3r. 
Information in reference 13 indicates that /3 is a function of the rotor angular velocity 
and the ratio s/a, or that /3 = Ku where K = F(s/a). Evaluating the data presented 
yields the following egressions for K (the coefficient equation can be expressed in terms 
of K, see table I): 
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K stator = 0 - 50185 -°- 6508 (f) 
K rotor = °* 49815 + 0-6508^ 


K stator = °- 48558 -°* 35778 (- 


K rotor = °* 51442 + °* 35778 ( 


laminar flow 


turbvilent flow 


102 



APPENDIX E 


SAMPLE THERMAL PROBLEM FOR COMPARISON WITH EXACT SOLUTION 


To demonstrate the accuracy of the program, the following sample problem (see 
fig. 15) was solved. 


6 in. 


T„ = 0° F; h - 0. 020 Btu/(sec)(ft 2 )(°F) 
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Figure 15. - Thermal sample problem. Solid cylinder with right end <z = 0) held at 200° F and left end (z = L) held 
at 0° F. Right side boundary elements (i.e., 14, 28, 42, 56, 70, 84, 98, and 112) held at 200° F. Left side ele- 
ments held at 0° F. 


For the finite cylinder with O < r < a and 0 ^ z < L, 

z = 0 held at 200° F (366 K) 
z = L held at 0° F (255 K) 

Forced convection into a medium at = 0° F (255 K) occurs at r = a. 
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For an exact solution, the temperature has to satisfy the following equation for an 
isotropic homogeneous solid: 


®!t + i ar + £r = 0 0 = r sa , 0 = z =L 

ar 2 r Sr az 2 


which is subject to the following boundary conditions: 

T = 0° F (255 K) at z = L and 0 < r < a 
T = 200° F (366 K) at z = 0 and 0 < r < a 


k— + hT = 0 at r = a and 0 < z < L 
ar 


If the following definitions are used: 

J Bessel function of first kind, zero order 
o 

«F derivative with respect to its arguments 
a root of equation aJ^(aa) +h c J Q (aa) = 0 
the exact solution is given by 


OO 



2h c J Q (ra n )sinh[(L - z)oJ 
ah l + oi* J o (aa n )sinh(La- n ) 


(El) 


For the numerical solution, a coarse mesh of 112 elements was chosen. As indica- 
ted in figure 15, the 96 elements between z = 0 and z = L had 0. 5-inch (1. 27-cm) ra- 
dial and axial dimensions. The 16 elements constituting boundary conditions had 0. 25- 
inch (0. 635-cm) axial dimensions. 

The temperatures were calculated by means of the program and of equation (El). 
The following values were required to obtain solutions: 
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a = 4. 0 in. (10 cm) h = 0. 020 Btu/(sec)(ft^)(°F); 0. 0098 cal/(sec)(cm^)(°C) 

L = 6.0 in. (15 cm) k = 0.00361 Btu/(sec)(ft)(°F); 0.0537 cal/(sec)(cm)(°C) 

The results are shown in figure 16 and table III. The maximum error using this 
coarse mesh is 15 percent. A finer mesh would be expected to show a lesser error. 
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Figure 16. - Thermal sample problem showing thermal distribution. 
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TABLE III. - THERMAL SAMPLE PROBLEM COMPARISON WITH EXACT SOLUTION 


Calculated Percent Element 


number temperature, temperature, error 


Element Exact Calculated Percent 

number temperature, temperature, error 
°F °F 


4.925 

4.571 

7. 19 

58 

4.097 

3.920 

4. 32 

14. 943 

12.682 

15.13 

59 

12.448 

10.883 

12.57 

25.465 

22.471 

11.76 

60 

21.270 

19.322 

9.16 

36.830 

33.150 

9.99 

61 

30. 903 

28. 592 

7.48 

49.372 

45.109 

8.63 

62 

41.707 

39.088 

6.28 

63.415 

58. 746 

7.36 

63 

54. 087 

51.250 

5.25 

79.246 

74.446 

6.06 

64 

68.487 

65. 576 

4.25 

97.082 

92.548 

4.67 

65 

85.394 

82.621 

3.25 

117.018 

113.272 

3.20 

66 

105.298 

102.978 

2.20 

138.969 

136.631 

1.68 

67 

128. 585 

127.185 

1.09 

162.620 

162.308 

.19 

68 

155.300 

155.304 

• 13 

187.408 

186.137 

.68 

69 

184.744 

183.395 

! .73 


14.691 

25.044 

36.240 

48.619 

62.515 

78.233 

96.016 

115.990 

138.096 

162.028 

187.216 


14.188 

24.202 

35.056 

47.102 

60.690 

76.162 

93.815 

113.843 

136.253 

160.767 

186.765 


13, 

.437 

22, 

.941 

33. 

.276 

44. 

.802 

57. 

.897 

72. 

,951 

90, 

,343 

110. 

384 

133. 

215 

158. 

647 

185. 

989 


4.532 
12.573 
22.281 
32.876 
44. 749 
58.303 
73.932 
91.985 
112.708 
136.134 
161.968 
186.004 
4.412 
12.242 
21.702 
32.041 
13.654 
56.952 
72.354 
90.252 
110.955 
134.579 
160.892 
185.584 
4.207 
11.678 
20. 716 
30.616 
41.774 
54.618 
69.603 
87.192 
107.809 
131.727 
158.871 
184.781 


3.695 
11.230 
19.210 
27.954 
37. 822 
49.226 
62.661 
78. 724 
98.134 
121.684 
149.987 
182.681 

3.227 
9.811 

16. 797 
24.477 
33.187 
43.330 
55.415 
70.114 
88.353 
111.444 
141.150 
178.863 
2.705 

8.227 
14.091 
20. 551 
27.899 
36.497 
46.814 
59. 503 
75. 544 
96.548 

125.551 

169.610 


3.552 
9.868 
17.536 
25.986 
35.598 
46.819 
60.178 
76.328 
96.063 
120.316 
150.003 
180.937 
3.106 
8.650 
15.391 
22.837 
31.344 
41.338 
53. 359 
68.116 
86. 586 
110.157 
140.806 
175.963 
2.574 
7.260 
12.941 
19.220 
26.411 
34.896 
45.168 
37.913 
74.148 
95.518 
124.981 
163.090 
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